00 



c^ 



Using effective field theory to analyse low-energy Compton 
scattering data from protons and light nuclei 

H. W. GrieBhammer^, J. A. McGovern^, 
D. R. Phillips^ and G. Feldman^ 

^Institute for Nuclear Studies, Department of Physics, 

The George Washington University, Washington DC 20052, USA 

^Theoretical Physics Group, School of Physics and Astronomy, 



(N 

o 

(N 

>% 

cd. The University of Manchester, Manchester, M13 9PL, UK 



^Institute of Nuclear and Particle Physics and Department of 
Physics and Astronomy, Ohio University, Athens OH 45701, USA 



May 7, 2012 

^ . Abstract 

P; Compton scattering from protons and neutrons provides important insight into the structure 

of the nucleon. For photon energies up to about 300 MeV, the process can be parameterised by six 
S^ i dynamical dipole polarisabihties which characterise the response of the nucleon to a monochromatic 

tH- I photon of fixed frequency and multipolar ity. Their zero-energy limit yields the well-known static 

cn ■ electric and magnetic dipole polarisabihties oei and /5a/i, and the four dipole spin polarisabihties. 

The emergence of full lattice QCD results and new experiments at MAMI (Mainz), HI7S at TUNL, 
and MAX-Lab (Lund) makes this an opportune time to review nucleon Compton scattering. Chiral 
f-^ • Effective Field Theory (^EFT) provides an ideal analysis tool, since it encodes the well-established 

CN . low-energy dynamics of QCD while maintaining an appropriately flexible form for the Compton 

amplitudes of the nucleon. The same x^FT also describes deuteron and '^He Compton scattering, 
using consistent nuclear currents, rescattering and wave functions, and respects the low-energy 
k> ■ theorems for photon-nucleus scattering. It can thus also be used to extract useful information 

^ I on the neutron amplitude from Compton scattering on light nuclei. We summarise past work in 

xEFT on all of these reactions and compare with other theoretical approaches. We also discuss 
all proton experiments up to about 400 MeV, as well as the three modern elastic deuteron data 
sets, paying particular attention to the precision and accuracy of each set. Constraining the 
A(1232) parameters from the resonance region, we then perform new fits to the proton data 
up to wiab = 170 MeV, and a new fit to the deuteron data. After checking in each case that 
a two-parameter fit is compatible with the respective Baldin sum rules, we obtain, using the 
sum-rule constraints in a one-parameter fit, a^^ = 10.7 it 0.3(stat) it 0.2(Baldin) it 0.8(theory), 
f^Mi ~ '^■^ T 0.3(stat) ± 0.2(Baldin) ± 0.8(theory), for the proton polarisabihties, and a^[ = 
10.9 ± 0.9(stat) ± 0.2(Baldin) ± 0.8(theory), /S^^^^ = 3.6 T 0.9(stat) ± 0.2(Baldin) ± 0.8(theory), 
for the isoscalar polarisabihties, each in units of 10~^ fm . Finally, we discuss plans for polarised 
Compton scattering on the proton, deuteron, "^He and heavier targets, their promise as tools to 
access spin polarisabihties, and other future avenues for theoretical and experimental investigation. 

Keywords: Compton scattering, proton, neutron and nucleon polarisabilities, spin polarisabilities, 
Chiral Perturbation Theory, Effective Field Theory, A(1232) resonance 
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1 Introduction 

Coinpton scattering has played a major role in the development of modern Physics. In 1871, Lord 
Rayleigh employed the recently discovered nature of light as an electromagnetic wave to demonstrate 
that the cross section for light scattering with frequency u from neutral atoms behaves as a oc w^, 
thereby explaining why the sky is blue [1, 2]. Fifty-six years later, Arthur Holly Compton won the 
Nobel Prize "for his discovery of the effect named after him" : X-rays have a longer wavelength after 
they are scattered from electrons, with the difference precisely predicted by a quantum treatment of the 
electromagnetic radiation and by the relativistic kinematics of the electron [3]. Compton's experiment 
thus provided a unified demonstration of two of the great advances in Physics in the early 20th century: 
relativity and the particle-like nature of light. 

In 1935, Bethe and Peierls performed the first calculation of Compton scattering from a nucleus: 
the deuteron [4]. The field gained momentum in the second half of the 20th century as it was realised 
that the photon provides a clean probe for Nuclear Physics. Its interactions with the target can be 
treated perturbatively, with the fine structure constant Oem as a small parameter.^ The leading term 
in the scattering of radiation from a nucleus of mass Mx and atomic number Z in the long-wavelength 
limit was first calculated by Rayleigh's student, J. J. Thomson: 

da _ Z^al^ 1 + cos^ 9 

dn~ Ml 2 ' ^ ' 

where 9 is the scattering angle. Thirring, and independently Sachs and Austern, showed that Eq. (1.1) 
is not renormalised non-relativistically [5, 6]. It is also recovered as part of a low-energy theorem for a 
spin-| target due to Low, Gell-Mann, and Goldberger [7, 8] which invokes only analyticity and gauge, 
Lorentz, parity and time- reversal invariance, and was generalised to arbitrary spin by Friar [9]. The 
theorem states that the nuclear magnetic moment is the only additional parameter needed to determine 
the leading spin dependence of the cross section. 

1.1 The importance of dipole electric and magnetic polaris abilities 

For a composite system, the first spin- independent piece of the Compton amplitude beyond the Thomson 
limit is parameterised by two structure constants: the electric and magnetic (scalar dipole) polarisabil- 
ities. Polarisabilities arise because the electric and magnetic fields of a real monochromatic photon 
with frequency uj displace the charged constituents of the system and thus induce charge and current 
multipoles, even if the target is overall charge neutral. The dominant contributions are typically an 
induced electric dipole moment rfjnd, often generated by separating positive and negative charges along 
the dipole component of the electric field E, and a magnetic dipole moment mjnd, often generated from 
currents induced by the dipole component of the magnetic field B. In addition, aligning microscopic 
permanent electric and magnetic dipoles in the external fields can generate mesoscopic induced electric 
and magnetic dipoles. The linear response in frequency space demonstrates that the induced dipoles 
are proportional to the incident electric and magnetic fields, i.e., 

rfind(w) = 47rasi(w) E{uj) , mind(w) = 47r/3Mi(c^) B{u) . (1.2) 

Neglecting recoil corrections, these induced dipoles then re-radiate at the same frequency u and with the 
angular dependence characteristic of El and Ml radiation, respectively. The proportionality constants 
in Eq. (1.2) are the electric dipole polaris ability aEi{u}) and the magnetic dipole polaris ability Pmi{^)- 
They characterise the strength of the dipole radiation relative to the intensity of the incoming fields 



^We use the Hcaviside-Lorentz system of electromagnetic units and /i = c = 1, so e = —^/Atto^. The factor Aw in 
Eq. (1.2) is absent in the Gaussian system, but present for SI units; cf. [10]. 



and vanish for objects with no internal structure. Since they lead to different angular dependences, 
they can be disentangled in the differential Compton scattering cross section at fixed frequency cu. 

Polarisabilities encode the temporal response of the target to a real photon of energy u and thus 
provide detailed information on the masses, charges, interactions, etc. of its active internal degrees of 
freedom. As an example, in the long- wavelength approximation, the Lorentz-Drude model assumes 
that the electric field of the photon displaces point-like, charged constituents g„ of mass ?n„, bound in 
classical harmonic oscillators with resonance energies E^ and damping factors r„: 

One can estimate uei (w) in the nucleon semiclassically by representing the constituents of its charged 
pion cloud as objects which are harmonically bound to the nucleon core with an eigenfrequency such that 
the pion cloud has the same root-mean-square radius as the nucleon, namely about 0.7 fm. Eq. (1.3) 
then leads to a value of aEi{^) ~ 10 x 10~^ fm . This model is unrealistic, of course, but the estimate 
is surprisingly close to the experimental value. Throughout this review, we therefore quote values for 
aEi{uj) and (3mi{'^) in the "canonical" units of 10"^ fm^, and for the dipole spin polarisabilities discussed 
below, in units of 10~^ fm'^. 

We also take the term "polarisabilities" to be the "Compton polarisabilities", defined in parallel to 
the above intuitive description by a multipole expansion of the Compton amplitudes, as elaborated in 
Section 2.2. In nonrelativistic Quantum Mechanics, they can be related to the spectrum of nucleon 
excited states \n) with energies En, e.g.: 

aE,{^) = 2^ ^-r;:: + ••• ' (^-4) 
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where D^ is the electric dipole operator. (In Eq. (1.4) we have not explicitly written subtle but im- 
portant corrections beyond non-relativistic second-order perturbation theory, which were emphasised 
by L'vov [11], and Bawin and Coon [12], and summarised by Schumacher [10].) Therefore, a£;i(a;) 
and /3jvfi(a;) are strongly influenced by the lowest state with the quantum numbers of an electric or 
magnetic dipole excitation of the nucleon. For uei-, this is indeed the ttN state. The Ml excitation of 
the A(1232) would appear to provide a sizable paramagnetic contribution to Pmi of order 10, but since 
the experimentally measured value is about an order of magnitude smaller, a diamagnetic contribution 
of similar magnitude but opposite sign exists whose precise nature is not yet determined. 

Such excitations also set the energies at which the polarisabilities are manifest in the single-nucleon 
cross section: cj > 50 MeV. In addition, Eq. (1.4) implies that the dynamical polarisabilities become 
complex once the first inelastic channel opens at the vrN threshold. Finally, polarisabilities are in general 
related to the dielectric function e{(jj) and magnetic permeability function /i(a;) of a macroscopic system. 
Since these, in turn, characterise optical properties, nucleon polarisabilities are related to the index of 
refraction and absorption coefficient of a bulk system of nucleons at a given frequency. 

Though polarisabilities are naturally defined as functions of the photon energy w, historically much 
of the emphasis in the context of the proton and neutron has been on the static polarisabilities uei = 
aEi{uj = 0) and /3mi = I^Miif^ = 0), which are often simply termed "the polarisabilities". For clarity, we 
therefore refer to the functions as energy- dependent or dynamical polarisabilities [13, 14]^. The static 
polarisabilities are formally and uniquely defined via the Compton scattering amplitudes, as detailed in 
Section 2.1. Along with the anomalous magnetic moment, they parameterise the deviation of the proton 



^For completeness, we note that the generalised polarisabihties of vntual Compton scattering are explored by an 
incoming photon of non-zero virtuality and can provide complementary information about the spatial charge and current 
distribution, see e.g. [15] for a review. 



Compton cross section from that of a point-like particle in a low-energy expansion which is valid up 
to photon energies of roughly 80 MeV. Static polarisabilities can be conceptualised, again up to subtle 
corrections [10-12], as the proportionality constants between induced dipoles and external fields which 
would be "measured" were the nucleon placed into a parallel-plate capacitor or a pure N-S magnet. 
They also enter in other processes with sensitivity to nucleon structure; in particular, the relation of 
Pmi to doubly-virtual forward Compton scattering has attracted recent interest in connection with the 
two-photon-exchange contribution to the Lamb shift in muonic hydrogen [16] (see also Ref. [17] and 
references therein), and with the nucleon electromagnetic mass shift, see most recently [18]. 

The Effective Field Theory (EFT) methods which are discussed in this review predict both static 
and dynamical polarisabilities. They can be used to extract the values of static polarisabilities from 
experimental data taken at energies too high for the low-energy expansion to be valid. We present a 
new EFT extraction of polarisabilities from world data in Sections 4 and 5. 

1.2 Compton scattering from nudeons: data, structure and analysis tools 

Therefore, in this review, we examine real Compton scattering from the simplest stable strongly- 
interacting systems, namely protons and light nuclei, in order to obtain information on the photon- 
nucleon scattering amplitude. Compton scattering from larger nuclei is reviewed in Ref. [19]. In 
Section 2.2, we provide the detailed relation of aEi{oj), /3mi{^) and the spin polarisabilities to the 
single-nucleon Compton amplitude. The spin polarisabilities have received much recent attention in 
both theoretical and experimental studies, since they are a low-energy manifestation of the spin struc- 
ture of the nucleon, parameterising its spin-dependent response to external electric and magnetic fields. 

How are polarisabilities explored? Until the last decades of the 20th century, all experiments inves- 
tigating Compton scattering from nucleons and nuclei employed bremsstrahlung beams. This created 
difficulties in accurately measuring the small (nb/sr) cross sections at energies where the polarisabilities 
are particularly relevant. The advent of photon tagging in the 1980s facilitated a clean separation of 
elastic and inelastic processes, enabling measurements of the proton cross section with good energy 
resolution at a; < 200 MeV. At the turn of the Millennium, this led to a wealth of data for the proton 
from Illinois [20], Saskatoon [21, 22] and MAMI [23, 24], and for the deuteron from Illinois [25], Saska- 
toon [26] and Lund [27]. Since most of these experiments were reviewed by Schumacher [10], we only 
summarise the data in Section 3, with particular attention to statistical and systematic errors. 

In parallel with these developments, aEi and Pmi were calculated in various theoretical models of 
nucleon structure [28-38]. Comparing these predictions with Compton-scattering data indicates how 
accurately these models describe electromagnetic excitations of the nucleon. Quark models which do 
not incorporate explicit pionic degrees of freedom tend to underpredict aEi and overpredict /3mi', see 
e.g. Ref. [33]. Computations of aEi and (3 mi in chiral quark models incorporate both long-distance 
(ttN) and short-distance (other excitations) physics [30, 39, 40]. Direct determinations of nucleon 
polarisabilities from lattice simulations of the QCD path integral now appear imminent [41], with 
results reported in quenched [42-44], partially quenched [45, 46, 49, 50] and even full QCD [47, 48]. 

New plans for Compton-scattering experiments on protons and light nuclei at MAMI, the HI7S 
facility at TUNL and MAX-Lab in Lund make this an opportune time to re-examine our knowledge of 
nucleon Compton scattering at energies up to a few hundred MeV. We therefore delineate in this review 
what is known about the Compton amplitudes of the proton and neutron. Equation (1.4) implies 
that the polarisabilities are dominated by the lowest nucleonic states, namely by vrN and A(1232) 
dynamics, while sensitivity to higher excitations is suppressed. This means that Compton scattering 
at low energies, u < 300 MeV, is dominated by long-distance properties of the nucleon. In particular, 
we note that the particles detected in the experiments are photons, nucleons and pions, not quarks 
themselves. Analysing Compton scattering at these energies in terms of quark degrees of freedom, such 
as in lattice QCD or models of nucleon structure, is thus not really profitable. Instead, such calculations 



can be tested against the constraints extracted from data using a theoretical approach that includes the 
pertinent low-energy dynamics, provided that it is sufficiently general to encompass the data without 
undue prejudice. Effective Field Theory (EFT) fits these requirements. 

1.3 The role of Effective Field Theory 

The basic principle of an EFT is that many phenomena can be economically — i.e. effectively — described 
in terms of entities which are not elementary. The fact that the details of nucleon structure are not 
probed at low energies suggests that a low-energy EFT which includes nucleons, pions and photons 
should be a useful tool to extract information on nucleon polarisabilities. 

In general, for an EFT approach to be successful, a separation of scales must exist between the 
energies involved and those required to excite the particular degrees of freedom of the system which 
are not treated dynamically. A famous example of an EFT is the Fermi theory of weak interactions, 
in which /? decay is described by a contact interaction between the neutron and its decay products 
e~, Vc and p. At energy scales of the order of a few MeV, the threshold for production of W and Z 
bosons is far off, and they can be "integrated out" to leave a simple energy-independent four-fermion 
interaction together with a series of further interactions, each of which is suppressed by powers of the 
small quantity ptyp/^w, where ptyp is the typical momentum of a decay product. Similarly, at low 
enough energies — energies below those where pion degrees of freedom become relevant — a theory of 
heavy, point-like nucleons should suffice to describe the interactions of nucleons with one another and 
with photons. This theory has come to be known as the "Pionless Effective Field Theory" (EFT(^)) of 
Nuclear Physics, and for NN scattering, it is equivalent to Bethe's low-energy Effective Range Expansion 
[51-59]. As in all EFTs, "point-like" and "pionless" does not imply the absence of effects such as a 
non-zero anomalous magnetic moment which are due to nucleon structure, pions and heavier particles. 
Instead, they are taken into account through Lagrangian parameters — so-called low-energy constants 
(LECs) — but are not explained within the EFT itself. 

The guiding principle in constructing an EFT Lagrangian is that all terms compatible with the 
symmetries of the underlying theory must be included, each proportional to an a priori unknown LEG. 
This infinite string of terms and couplings is organised according to the power of Ptyp/A that each 
operator contributes to amplitudes, with ptyp the typical momentum of the process, and the breakdown 
scale A set by the mass of the lightest omitted degree of freedom {rriT^, for example, in the case of 
EFT(]f)). Since all terms are included, the counterterms needed to renormalise the divergent loops at 
a given order in Ptyp/A are automatically present. Unless the theory has a low-energy bound state, the 
renormalised loop contributions are suppressed, typically by pi /A? for each loop. Therefore, only a 
finite number of terms in the Lagrangian need to be considered when working to a particular order in the 
small, dimensionless parameter ptyp/A. While the theory is thus not renormalisable in the conventional 
sense, only a finite and usually small number of terms are needed for renormalisability to a given order 
in ptyp/A. In addition, a rigorous assessment of residual theoretical uncertainties can be made for any 
process by estimating the accuracy of its momentum expansion. 

Though the number of possible terms (and hence the number of LECs) in the Lagrangian grows 
rapidly with the order, any given process usually involves only a few of them. Once determined by one 
piece of data, an LEG enters in the prediction of other observables. Some LECs are related to familiar 
properties of the particles involved, like charge, mass, anomalous magnetic moment, decay constant, 
etc., but others are less easy to fix and interpret. Some LECs which govern pion interactions are now 
being computed by direct lattice simulations of QCD [60-62]. But even when LECs can be derived 
from the underlying theory, calculations of more complex processes are often more tractable in the EFT 
framework, as we shall see here for Compton scattering from protons and light nuclei. 

Compton scattering in EFT(^) is discussed in Section 5.4.1. For the proton, it is of limited use, 
as it amounts to an expansion of the amplitude in powers of u /rriT, <^ 1. Since only a small fraction 

6 



of the data is in the regime where this expansion is vahd, the resuhing errors on the polarisabihties 
are inevitably large. In the two-nucleon sector, though, the situation is different. The low-energy NN 
scattering amplitude is given as an expansion in powers of momenta, with the LECs determined by 
the scattering length, effective range, etc. Nuclear binding effects and photon-nuclear interactions are 
then fixed and the leading deuteron Compton amplitudes are predicted with no free parameters, as 
demonstrated in 1935 in a calculation by Bethe and Peierls [4]. 

The EFT(^) expansion breaks down for typical momenta of order ?ti^, and above this point, the 
pion itself must be included. While an EFT including the pion must exist because the next excitation, 
the A(1232), is far enough away to provide a separation of scales, there is no guarantee that it is viable. 
The size of the pion-nucleon coupling constant, g-j^'NTSi ~ 13, suggests that multi-loop diagrams including 
dynamical pions may not be suppressed. However, two crucial aspects make the theory manageable. 
First, the fact that the pion is much lighter than other hadrons is now understood as a consequence 
of the spontaneously broken (hidden) chiral symmetry of QCD. The up and down quarks are nearly 
massless on the scale of typical QCD energies. If they were actually massless, the Lagrangian would be 
invariant under independent isospin rotations of the left- and right-handed quarks, SU(2)l x SU(2)r. 
However, only the vector (isospin) subgroup, SU(2)v, is manifest in the hadron spectrum and the full 
symmetry is hidden in the physical QCD vacuum. The pions are then identified as the three Nambu- 
Goldstone bosons corresponding to the three axial rotations which are symmetry operations on the 
QCD Lagrangian, but not on the QCD vacuum. The small non-zero quark masses lead to pion masses 
which are again much smaller than typical QCD scales, ttt,^ <^ A-;,.. The fact that the pion is a (pseudo- 
)Nambu-Goldstone boson leads to the second key point: in the "chiral" or zero-quark-mass limit, soft 
pions decouple so that their interactions with one another and with other hadrons vanish linearly with 
their momentum. 

This provides a perturbative expansion of amplitudes in powers of -P/A^ = {ptyp,mT,)/A^ ^ 1, 
with the light EFT scales being the pion momenta and mass, and A-^ being the scale associated with 
hadrons which are not explicitly included in the EFT, A^ Ri rrip. This EFT is known as "chiral EFT" 
(xEFT). The version without explicit A(1232) degrees of freedom is often referred to as "Baryon Chiral 
Perturbation Theory" (B^PT), and the purely pionic one as Chiral Perturbation Theory (xPT). Details 
pertaining to Compton scattering are discussed in Section 4.2, including the special role of the A. There 
we also show that photons are included in xEFT partly by invoking minimal substitution and partly by 
including the field strength tensor F^j, as a building block in the Lagrangian. Since the latter generates 
photon couplings which are not constrained by gauge invariance, new LECs enter, including ones which 
can (at successively higher orders) be related to the anomalous magnetic moment and the charge radius 
of the nucleon, and to non-chiral contributions to polarisabihties. 

The resulting EFT framework provides predictions for the low-energy interaction of pions, photons 
and nucleons through calculations consistent with the known pattern of QCD symmetries and their 
breaking. Furthermore, because it is a quantum field theory, xEFT incorporates the requisite conse- 
quences of unitarity and Lorentz invariance at low energies. It is this framework which we use to analyse 
proton Compton scattering in Sections 4.2 and 4.4. In Sections 4.1, 4.3 and 4.5, we also compare xEFT 
to other approaches for analysing low-energy Compton scattering, most notably dispersion relations 
(DRs). (Full details of DRs are given in the review of Drechsel et al. [15].) The main differences be- 
tween xEFT and DRs lie in the careful incorporation of chiral constraints in the former, its stringent 
agnosticism regarding high-energy details, and the fact that the presence of a small parameter allows 
an a priori estimate of residual theoretical uncertainties. 

1.4 The elusive neutron: Compton scattering from light nuclei 

Experiments on the proton reveal only half of the information in Compton scattering on the nucleon. 
The 14.7 minute lifetime of the neutron, coupled with the relative weakness of its electromagnetic 
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interactions, means that the neutron Conipton amphtude must be inferred indirectly. While some 
results exist for scattering neutrons directly in the Coulomb field of heavy nuclei, more accurate data 
are available for Compton scattering in few-nucleon systems, where nuclear effects can be precisely 
calculated and taken into account in the analysis. However, it is clearly advantageous to use the same 
theoretical framework for both the nuclear and photon-nucleon dynamics. The xEFT low-momentum 
expansion again provides a controlled, model-independent framework for subtracting the nuclear binding 
effects and analysing the available elastic deuteron scattering data. Both the number and quality of 
these data, reviewed in Section 3.4, are appreciably inferior to the proton case, since the experiments are 
markedly harder. The substantial progress made to determine the neutron polarisabilities via this route 
is reviewed in Section 5. The key to extracting neutron polarisabilities from few-nucleon targets is that 
the coherent nature of deuteron Compton scattering allows us to observe the photon-neutron amplitude 
through its interference with the proton and meson-exchange amplitudes, while the Thomson limit 
imposes a stringent constraint on the few-nucleon amplitude that can be used to check these non-trivial 
calculations. 

xEFT work on the deuteron originates from Weinberg's seminal papers [63], where the nucleon- 
nucleon potential is computed up to a fixed order in the momentum expansion and then iterated using 
the Schrodinger equation to obtain the scattering amplitude. The resulting wave functions are combined 
with operators for Compton scattering derived in xEFT. The photon-nucleon operators are given by the 
xEFT photon-nucleon amplitudes described above. However, xEFT also expands the nuclear current 
operators in powers of P. This has the crucial advantage that the chiral dynamics which drives low- 
energy Compton scattering is treated consistently in both the single- and few-nucleon operators, so 
that xEFT results for deuteron Compton scattering can be assessed for order-by-order convergence, 
too. Meanwhile, Compton scattering from the deuteron in quasi-free kinematics for the neutron was 
measured in Refs. [64-68], and values for a^l were extracted. These experiments are discussed in 
Section 3.5, but xEFT has not yet been extended to 7d — )■ 7np. Model calculations of elastic and 
inelastic scattering on the deuteron are briefly reviewed in Section 5.5. xEFT is also used to produce 
the first calculations, in any framework, of elastic scattering on ^He [69-71]; see Section 5.6. Since ^He 
is doubly charged, its cross section is larger by a factor of up to 4 compared to the proton or deuteron. 
In addition, polarised ^He is interesting since it serves as an effective polarised neutron target. 

1.5 Results and the future 

In Sections 4.4 and 5.3, we apply the xEFT methodology to the proton and deuteron databases, respec- 
tively, and present new central values and uncertainties for the proton and neutron dipole polarisabilities, 
as well as a detailed comparison of the xEFT predictions with data. Here we preview our best values 
of the proton and isoscalar (i.e. average nucleon) scalar dipole polarisabilities in a two-parameter fit: 



a^J = 10.5 ± 0.5(stat) ± 0.8(theory) , /jJpI = 2.7 ± 0.5(stat) ± 0.8(theory) 
ag = 10.5 ± 2.0(stat) ± 0.8(theory) , /3[^\ = 3.6 ± l.O(stat) ± 0.8(theory) . 

In a one-parameter fit employing the Baldin sum rules for the proton and isoscalar nucleon: 

a^^l = 10.7 ± 0.3(stat) ± 0.2(Baldin) ± 0.8(theory) 
/3[PJ = 3.1 T 0.3(stat) ± 0.2(Baldin) ± 0.8(theory) 
ai^i = 10.9 ± 0.9(stat) ± 0.2(Baldin) ± 0.8(theory) 
/3^\ = 3.6 T 0.9(stat) ± 0.2(Baldin) ± 0.8(theory) . 



:i.5) 



:i.6) 



In Section 6, we close this review with a discussion of future experiments which hold the promise of 
improving the Compton database and the determination of polarisabilities. We also describe upcoming 



experimental efforts to extract the still relatively unexplored spin polarisabilities. Finally, we outline 
the anticipated xEFT developments that will help refine the analysis of these forthcoming data. 

2 Foundations 

2.1 Overview 

We begin by parameterising the T-matrix for Compton scattering of a photon of incoming energy u 
from a nucleon with spin a/2 by six independent invariant amplitudes [72], which read in the operator 
basis first appearing in Ref. [73] and usually employed in xEFT: 

T(w, z) = Ai(a;, z) {e'* ■ e) + ^(a;, z) {e'* ■ k) {e- k') 

+i ^3(0;, z) a ■ (?'* X e ) + i ^4(0;, z) a ■ (k' x k) {e'* ■ e) 



+iA^{u,z) a ■ 
+1^46(0;, 2;) a ■ 



e'* X k\ (e*- k') - (ex k'\ {e'* ■ k) 
(e'* X k'\ {e- k') - (exk) {e'* ■ k) 



(2.1) 



where k {k') is the unit vector in the momentum direction of the incoming (outgoing) photon with 
polarisation e (e*'*), 9 is the scattering angle, and z = cos9. This form holds in the Breit and centre- 
of-mass (cm) frames. The first two amplitudes are spin-independent, while the other four parameterise 
interactions with the nucleon spin. The amplitude is related to the differential cross section by 



dcj 
d?^ 



n.me \T\' , (2.2) 



frame 



where $framo is a frame-dependent flux factor which tends to l/(47r) at low energies; see also Section 2.3. 
For forward scattering, k = k', and only the structures of Ai and A3 survive. The behaviour of these 
amplitudes as w — t- is determined by low-energy theorems (LETs) which rely only on analyticity and 
gauge, Lorentz, parity and time-reversal invariance [7, 8] 

iimTiu, 1) = -^ e'* ■ e- i^a ■ {e'* x e) + 0(c.^) , (2.3) 

where —eQ is the nucleon charge, k is the anomalous magnetic moment (in nuclear magnetons) and 
Mn is the mass of the target. At w = 0, we recover the spin- and energy-independent Thomson term, 
and hence the corresponding cross section of Eq. (1.1). 

Eq. (2.3) may also be obtained in the limit a; — )■ from the pole diagrams in a calculation of 
Compton scattering using photons coupled to a Dirac nucleon via its charge and anomalous magnetic 
moment (with the spinors normalised to uu = 1), as shown in diagram (a) of Fig. 2.1. 

^k k' k' k\ y-k' 



(a) N jf- / , (b) 



TlP 




Figure 2.1: (a) Nucleon Born graphs for Dirac nucleons; (b) pion Born graph; (c) structure contribution. 

As explained in the Introduction, our primary interest is to obtain information on the structure of 
the two-photon response of the nucleon, which requires going beyond the single-photon processes of the 

9 



Born terms. It is therefore useful to separate the amphtudes into "non-structure" and "structure" parts 

Ai{u, z) = Af°™(^, z) + M^. z) ■ (2.4) 

Identifying the former with the pole (Born) contribution mentioned above, while not a unique choice, 
ensures that the LETs are satisfied by these terms alone. It also agrees with usage in dispersion 
relation (DR) calculations (see Section 4.1). There is another known contribution to spin- dependent 
scattering which comes from the t-channel exchange of a neutral pion, as shown in Fig. 2.1(b), which 
first contributes at 0{uj^); we choose to include this term in the Born or "non-structure" part as well, 
but the choice is not universal. The other term, Ai, represented by Fig. 2.1(c), parameterises deviations 
from the "known" part of the amplitude and describes the nucleon as a polarisable object. 

Expanding both nucleon Born and structure amplitudes in the Breit frame to order cu^, we have 

-)2^2 „2 , ,2 



M{^.z) = -^ + ^((Q + «:)2(i + ^) _ Q2j(i _ z)+A7rco'{aE, + z Pmi) + 0{uj' 



v42(w, z) = ^^«(2Q + ^)z - ^nu^f^Mi + 0{u^) 

o2 



AsiuJ, z) = ^ {Q{Q + 2«:) - (g + nfz'^ + Af + Awu^^i - (72 + 274)^) + 0{uj' 



'N 



A^iu, z) = -j^iQ + 't) + 4vra; 72 + Oiu") 

2 
A^ico, ^) = -^iQ + f^f + < + 4vra;'74 + 0{lo') 



Ae{oo, z) = -^]^QiQ + «:) + ^6 + 4vra;S + 0{uj') , (2.5) 

where the tt" pole amplitudes are, with t = luj^iz — 1) and tj, the third Pauli matrix in isospin space: 

.^0 _ e^g^^^ufi{z - 1) A-K^ _ _ a-k"" _ ^ e^gnN^^'^ ,^ p. 

To this order, the only difference in the cm frame is that A2 has an additional Born term e^Q'^u/M^. 
The omitted terms start one order lower in the cm frame (e.g. 0{uj^) in Ai), because the amplitude 
does not have manifest crossing symmetry, whereas in the Breit frame it does. The "structure" parts 
of the amplitude are represented at this order by aEi and (3 mi, namely the static spin-independent 
polarisabilities discussed above, and by the analogous spin polarisabilities 7^ [74]; see Section 2.2. The 
canonical units are 10~^ fm for the former and 10~^ fm for the latter. We denote proton and neutron 
polarisabilities by a superscript (e.g. a^/) and define isoscalar and isovector polarisabilities as averages 

and differences, such that a%\ = |(q;^i + (^ei) ^^"^ ^ei ~ ^ei ~ ^ei- 

For forward scattering, only aEi + /?mi and 70 = 7i — 72 — 274 augment the LETs of Eq. (2.3), while 
the 7r° pole terms do not contribute. Since forward scattering amplitudes can be related to inelastic 
cross sections via the optical theorem, two sum rules can be constructed for these quantities: 

00 00 
1 /• , ariuj') 1 f , 0-1/2(0;') -(^3/2(^0 
«ei + /3mi = ^y do; — ^ , 70 = ^j do; -^ —^ , (2.7) 

where ax is the total cross section for an unpolarised target and ai/2 (o"3/2) is the cross section for a 
photon-nucleon system of total helicity 1/2 (3/2). The integrals are over lab energies, starting at the 
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pion photoproduction threshold. The first of these is known as the Baldin Sum Rule [75, 76]. We adopt 
the evaluation of Olmos de Leonet al. [24] for the proton and that of Levchuk and LVov [77] for the 
neutron: 

«S + /3mi = 13.8 ± 0.4 , «S + /3S = 15.2 ± 0.4 , (2.8) 

which combine to give the isoscalar sum-rule value, with errors added in quadrature: 

ag + /3iJi = 14.5 ± 0.3 . (2.9) 

Using information on pion photoproduction multipoles, together with a parameterisation of data at 
intermediate energies and a Regge form at higher energies, Babusci et al. [78] obtained a^l + /3|'^{ = 
13.69 ± 0.14. However, using different parameterisations for ar, Levchuk and L'vov [77] obtained a 
central value of 14.0. The evaluation quoted above is more recent and more conservative [24]. Given 
the uncertainties associated with extracting neutron multipoles from deuterium data, perhaps it is not 
surprising that evaluations for the neutron are more variable. We regard the direct use of deuterium 
photodisintegration data above pion threshold as ill-advised (c.f. Ref. [78]). The value we quote uses 
neutron photoproduction multipoles obtained from proton ones via isospin considerations in a more 
reliable approach, although significant model dependence is still present [77]. 

Only the combinations aEi — Pmi and 77r = 7i + 72 + 274 enter for backward scattering. Because 
we have separated out the pion-pole contribution, our values for 'y^ differ from those obtained when 
it is included. The difference is — e^(77rNN/(8vr^/^MNm^o) = —46.4, using the proton value gl^^/Air = 
13.64 [79]. 

Various low-energy cross sections have been constructed in the literature. All of these can be 
found by using various approximations to Eq. (2.5) in Eqs. (2.1) and (2.2). The simplest is the Klein- 
Nishina formula of a point-like Dirac particle, with k = [80]. If k 7^ is included but the target 
is still structureless and terms of order u^ are discarded, one obtains the Powell cross section [81]. 
The Petrun'kin cross section additionally allows the inclusion of terms of order u"^ arising from the 
interference of the leading structure contributions oei and (3 mi with the Thomson term in |v4ip and 
hence is complete at C(a;^) [82, 83]. None of these includes spin polarisabilities, the vr" pole, or the 
energy dependence of the polarisabilities. We shall say more about their impact in the next section. 

2.2 Polarisabilities from a multipole expansion 

Many multipoles are induced inside an object that interacts with an electromagnetic field of frequency 
u. Each oscillates with that same frequency and thus emits radiation with a characteristic angular 
distribution. The proportionality constant between each photon field and the corresponding induced 
multipole moment is called a polarisability; each is an energy- dependent function which parameterises 
the stiffness of the internal degrees of freedom with particular quantum numbers against deformations 
of a given electric or magnetic multipolarity and energy. In this section, we generalise the picture 
presented in the Introduction to consider polarisabilities beyond the scalar dipole ones. 

Hildebrandt et al. [13, 14] used the formalism of an energy-dependent multipole analysis established 
by Ritus [84-86] and summarised in Ref. [87] to define energy-dependent polarisabilities. Here we 
proceed differently, constructing the first few multipoles via the most general field-theoretical Lagrangian 
which describes the interactions between a nucleon field A^ with spin (j/2 and two photons of fixed, non- 
zero energy u and definite multipolarities. This includes the structure effects, i.e. the local coupling of 
the two photons to the nucleon. Taking into account gauge and Lorentz invariance, as well as invariance 
under parity and time-reversal, the interactions with the lowest photon multipolarities are 

£poi = 27r N^[aEi{uj) E^ + Pmi{uj) B^ + -fEiEi{oo) a ■ {E x E) (2.10) 

+ lMiMi{uj) a ■ {B X B) - 2'jMiE2{uj) ai Bj Eij + 2-fEiM2{u;) (^i Ej Bij + . . .] N , 

11 



with Tij = \{diTj + djTi), T = E,B. These terms are straightforward extensions to the effective 
Lagrangian of zero-energy scattering in Refs. [88, 89]. The photons couple electrically or magnetically 
{X,Y = E,M) and undergo transitions XI — )■ Yl' of definite multipolarities / and I' = I ± {0,1}. 
The interactions are unique up to field redefinitions using the equations of motion. Dipole couplings 
are proportional to the electric and magnetic field directly, or to their time derivatives. Quadrupole 
interactions couple to the irreducible second-rank tensors Eij and Bij. Eq. (2.10) lists all contributions 
with coupling to at least one dipole field. Polarisabilities of higher multipolarity, e.g. the electric and 
magnetic quadrupole polarisabilities [88, 89], are denoted by ellipses. Thus far, such terms are not 
relevant for Compton scattering below 300 MeV [14, 90-93]. 

The two-photon response of the nucleon in the dipole approximation is therefore characterised by the 
six linearly independent, energy-dependent polarisabilities of Eq. (2.10). The spin-independent terms 
are parameterised by the two scalar functions already encountered in Section 1.1: the electric dipole 
polarisabihty aEi{uj), and the magnetic dipole polarisability f^uii^)- The four spin polarisabilities 
parameterise the response of the nucleon spin to an external field. The two corresponding to dipole- 
dipole transitions, ■jEiEii^) and ■Jmimi{^), are analogous to the classical Faraday effect related to 
birefringence inside the nucleon [94]. They describe how an incoming photon causes a dipole deformation 
in the nucleon spin, which in turn leads to dipole radiation. The two mixed spin polarisabilities, 
7eim2(i^) and 'Jmie2{^), encode scattering where the angular momenta of the incident and outgoing 
photons differ by one unit. In principle, the polarisabilities can be defined in any coordinate system in 
which the initial and final photon energies are identical. In practice, the centre-of-mass frame is usually 
used.^ 

We can now translate the interactions (2.10) into contributions to the structure amplitudes Ai{u, z): 

= 47r [aEi{uj) + z(3mi{uj)] uj^ + ■ . . 

= -4:71 [jElEliuj) + Z-fMlMl{^) + lElAni^^) + Z -fMlE2{uj)] UJ^ + . . . 

= An [-7a/ia/i(c^) + 1mie2{uj)] uj^ + ■ ■ ■ (2.11) 

= 47r7A/iA/i(w)a;^ + ... 
= 47r7£;iM2(w)u;^ + ... , 

where the dots refer to omitted higher multipoles. The relations between the static polarisabilities of 
the multipole expansion and those of Ragusa [74] defined in Eq. (2.5) are 

7l = —lElEl — 1EIM2 , 72 = 'JM1E2 — lAIlMl , 73 = 1EIM2 , 74 = IMIMI , . 

[2.12) 

70 — —1M1E2 — 7M1M1 — 7-BlBl — 7-B1M2 , In — 1M1E2 + IMIMI — 1 ElEl — 1 EIM2 ■ 

The first of the photon multipolarities in the subscripts of the spin polarisabilities is sometimes dropped 
in the literature (so 'Jm2 = 1eim2 etc). 

We reiterate that polarisabilities are identified by a multipole analysis at fixed energy, i.e. only by the 
angular and spin dependence of the amplitudes. Eq. (2.11) emphasises that the complete set of energy- 
dependent polarisabilities does not contain more or less information than the untruncated Compton 
amplitudes Aiiuj^z). However, the information is more accessible, since any hadronic mechanism and 
interaction leaves a characteristic signature in a particular multipole polarisability, as discussed in the 
Introduction. Thus, when the multipole expansion is truncated after the dipole terms, determining the 
six energy-dependent dipole polarisabilities is reduced, in principle, to an energy-dependent multipole 



A(w, 


z) 


Mioj, 


/ 


M^. 


z) 


M^. 


z) 


A(w, 


) 


A(w, 


) 



•^Thc dynamical polarisabilities introduced here via Eq. (2.10) differ from those given in Rcfs. [13, 14] by a factor of 
•^/s/A/n- The polarisabilities are linear combinations of the multipole moments /^y of the Compton amplitudes; the 
details, including the relevant projection formulae, are given in Ref. [14]. 
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analysis of the Compton scattering database. Such proof-of-principle results were reported in [92, 95, 96] 
but suffer at present from rather large error bars. 

At very low energies, the functions encoding the dynamical polarisabilities can be approximated 
by their zero-energy values and therefore some experiments (most recently Ref. [20]) have used the 
Petrun'kin formula to extract the static polarisabilities a^{ and (3^{ directly from their data. Indeed, 
the low-energy expansion of the cross section could conceivably be extended to higher powers in cj^. At 
fourth order, not only do the spin polarisabilities enter, but also the next terms (slope parameters) in 
the expansion of aEi{(^) and /3mi{^)'- 

\im aEi{uj)=aEi + uj^aE^ + 0{uj^) , \im ^mi{uj) = Pmi + oo^Pm^ + 0{uj^) . (2.13) 

Static values of the next multipoles, the scalar quadrupole polarisabilities, also enter at (9(0;^). (Contri- 
butions to T at Lo"^ and co^ were discussed by Babusci et al. [88] and Holstein et al. [89].) However, the 
convergence is governed by the pion-production threshold (the first non-analyticity) and the expansion 
is thus in powers of u/rriT^. With the slope correction to the scalar polarisabilities of size (cj/?72,r)^, this 
leads to a correction of about 10% for photon energies as low as 50 MeV. The expansion is therefore 
useless where most high-accuracy data are taken. 

Consequently, modern extractions of the static polarisabilities choose a different route. It is assumed 
that the energy dependence of the polarisabilities, and hence that of the amplitudes Ai, is adequately 
captured in some framework (e.g. dispersion relations, xEFT). With the long-range part of the interac- 
tions thus fixed, one fits up to six low-energy constants which encode the short-distance dynamics and 
thereby determines the static polarisabilities from data. This is the approach taken in this review. 

2.3 Kinematics 

Consider Compton scattering on a nucleus X with mass Mx, 7(0;, k) X{E,p) — )■ 7(0;', k') X(£", jf), with 
E and p {E' and p^) the kinetic energy and momentum of the target before (after) the reaction. So 
far, uj denoted a generic photon momentum and 9 a generic scattering angle. We now add subscripts 
to differentiate between the centre-of-mass (cm), laboratory (lab) and Breit frames, using relativistic 
kinematics throughout. The coordinate axes are specified as follows: the incident photon beam direction 
defines the z-axis, and the x2;-plane is the scattering plane, with the y-axis perpendicular to it. 

In the centre-of-mass (cm) frame, the total energy is the square-root of the Mandelstam variable s. 



s = uj,r^ + ^Ml + ul^ . (2.14) 

On the other hand, experiments are performed in the lab frame, where the incident and outgoing photon 
energy are related to one another by a recoil correction, and to Wcm; as follows: 

Mx Mx + Wlab ( 1 - COS fc^iab ) 



The scattering angle transforms as 



cos^iab-/? ,^^„. 

cos9c„, = - -— , (2.16) 

1 - /3 cos aljab 



where /3 = ujiah/i^iab+Mx) is the relative velocity between the cm and lab frames. The frame-dependent 
flux (phase-space) factor for cross sections (see Eq. (2.2)) is 

$c. = ^ , $lab = ^ . (2.17) 
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It is also useful to introduce the Mandelstam variables t and u = [s — u)/(4Mx): 

t = 2a;iabWi'ab(cos6'iab - 1) = 2ujl^{cos9cm - 1) = 2a;|^^it(cos6'Brcit - 1) , i^ = |(wiab + ^I'ab) ' (2-18) 
where we have also listed variables in the Breit frame, with 

WBreit = , , COs6'Breit = 1 T^ (1 " COS 6'cm) • (2.19) 

V4M2 - t W^eit 

For uj/Mx — )■ 0, i.e. small photon energy or large target mass, the three coordinate frames coincide. 
The Breit and lab frames coincide for ^ = 0°, and the Breit and cm frames for 6 = 180°. 

In the Breit or "brick-wall" frame, the photon transfers no energy and the target recoils with the 
magnitude of its momentum unchanged but its direction exactly reversed, p^^it — ~PBrcit- This has 
the advantage that the Compton amplitude is manifestly crossing-symmetric, i.e. invariant under the 
interchange of initial and final states: WBrcit ^ — i^Brcit, ^Brcit ^ ~^Brcit ^ ^ ^*- By inspection of 
Eq. (2.1), the spin-independent amplitudes Ai^2 are even in cuBreit, and the spin-dependent ones, A^^q, 
are odd. 

2.4 Observables 

We can now relate the six independent amplitudes of Eq. (2.1) to scattering observables. A complete 
classification of unpolarised, single- and double-polarisation observables for a spin-| target (nucleon or 
■^He) by Babusci et al. [88] demonstrated that experiments in which up to two polarisations are fixed can 
be described by four independent observables below the first threshold and four more observables above 
it. Here, we list some of the combinations which either have been or will soon be explored experimentally 
or theoretically: unpolarised, linearly or circularly polarised beams and unpolarised or polarised targets, 
without the detection of final-state polarisations. We define them as they are measured, and — with the 
exception of the differential cross section — refer to the literature for formulae which relate them to the 
Compton amplitudes Ai. Observables with a vector-polarised nucleus can be understood by replacing 
the spin-polarised nucleon by the polarised nucleus. These will be discussed in Section 6.1. 

In order to take into account nuclear binding for the deuteron and '^He, it is convenient to use a 
helicity basis for these targets: 

A{Mf,Xf;Mi,X,) = {Mj,\j\T\Mi,\.:) , (2.20) 

where Xi/f = ± is the circular polarisation of the initial/final photon, and Mj// is the magnetic quantum 
number of the initial/final target spin, i.e. Mj/j e {0; ±1} for the deuteron and Mi/j G {±|} for ^He. 
A target of spin ^x has [2(photon helicities) x (2S'x + 1) (target spins)]^ (in and out state) amplitudes, 
but parity and time reversal leave only 2(S'x + l)(2S'x + 1) independent components: 4 spin-0 helicity 
amplitudes become 2 independent ones, Ai^2', 16 helicity amplitudes for spin-| reduce to the 6 of 
Eq. (2.1); and 36 amplitudes for spin-1 yield 12 independent structures, constructed e.g. by Chen et 
al. [97]. 

Unpolarised beam and target: In this case, the only observable is the differential cross section. 
For the nucleon, it is obtained from Eq. (2.2) after averaging over the initial target spins and photon 
polarisations and summing over the final states; one finds in the basis of Eq. (2.1) (see e.g. Refs. [98, 
Chapter IV.2] and [92, Chapter 4.1]): 

\Tf = 1 |Ai|2 (1 + ^2) + i lAgp (3 - z^) + (1 - ^2) [4Re[A*A6] + Re[AlA^ + 2A*A, - A*^] z 

+ (1 - z') [i \A2\' (1 -z') + l \A,\' (1 + z') (2.21) 

+ l^sp (1 + 2^2) + 3 |Ag|2 + 2 Re[A; {A^ + 3^)]^ + 2 Re[AlA5]z^ 
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In the helicity basis, the unpolarised differential cross section is built with the ffux factors of Eq. (2.17) 
by summing over all combinations of incident and outgoing quantum numbers (Mj, A/; Mj, Aj) and 
including a symmetry factor from averaging over the initial target and photon polarisations: 



da 

dn' 



u,9) 



frame 



2{2Sy 



frame, X 



J2 \A{Mf,Xf;M,,X, 

Mf,Mt]^f,\i 



(2.22) 



Polarised beam, unpolarised target: The cross section with a circularly polarised photon beam 
of arbitrary helicity is half of the unpolarised one and thus provides no additional information. With 
a linearly polarised beam, two observables can be constructed, see Fig. 2.2, with corresponding exper- 



iments approved at HI7S [99, 100]. The cross section for an incoming photon is denoted by 



da 

'dn 



lin 



^{y) 



d0_ 



lin 1 _^ _^ 



da 



nlin 



Figure 2.2: (Colour online) Observables for linearly polarised photon incident on unpolarised target. 

for an incoming photon polarised parallel (perpendicular) to the scattering plane. Their sum gives the 
unpolarised cross section, and the difference. 



S3(w,l 



"da" 
dVL 


lin 

X 


"da" 
dVL 




(2.23) 



is the beam asymmetry^. Its relation to the deuteron helicity basis is given in Refs. [71, 93, 101]. 

For the forward and backward spin polarisabilities, 70 and 7^, we already saw that the multipole 
expansion is a convenient tool to identify configurations in which a specific polarisability is isolated 
or suppressed. Such configurations were found for the proton by Maximon [102]. The interaction 
in Eq. (2.10) parameterised by aEi{t^) vanishes when the polarisations of the incoming and outgoing 
photons are orthogonal, e.g. when a photon which is linearly polarised in the scattering plane scatters at 
90° as in [^] (6'Breit = 90°), see Fig. 2.3. In that case, since the incoming magnetic field is orthogonal 



c> 



o> 



k' J, 




Figure 2.3: (Colour online) Left: Configuration for which an induced electric dipole cannot radiate an 
El photon to an observer ("eye") at ^Breit = 90°; right: same for radiating an Ml photon. 



to the scattering plane, the induced magnetic dipole radiates most strongly at 90°, providing maximal 
sensitivity to (3mi{uj)- Similarly, [^] (6'Brcit = 90°) is independent of /3mi(i^) but maximally sensitive 
to aEi{uj). When the nucleon is embedded in a nucleus, the relative motion of the 7N system may 
complicate this analysis. 



*The subscript is omitted in Ref. [101]; the symbol used in Ref. [93] is 11 



lin 
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Polarised beam, polarised target: Figure 2.4 depicts double-polarisation observables with circu- 
larly polarised photons, as an example of observables that will be explored in the future; see Section 6.2. 
The target can be polarised in the scattering plane along x or along the beam direction, z. Cross-section 



A"'": ^ 






Figure 2.4: (Colour online) Observables for circularly polarised photon incident on polarised target, 
differences can be defined by flipping the target polarisation: 

The first arrow of the subscript denotes a positive beam helicity, the second the target polarisation. 
A^''''^ compares a target polarised along +x vs. — x, i.e. perpendicular to the beam direction but in the 
scattering plane. Similarly, A^^'''^ is the difference with the target polarised parallel vs. anti-parallel to 
the beam helicity^. Both observables change sign for left-circularly polarised photons (negative beam 
helicity) . 

Polarisation asymmetries are defined as the ratio of the cross-section differences A to their sums^: 

(An_\ _ (Aa_\ (Aa_\ _ ( Aa_\ 

^x /da_\ I ( da_\ ' ^z / da_\ . ( da_\ ' l^Z.ZOJ 

Except for E^'^^'^ for a spin-| target, the denominators are not the unpolarised cross sections [71, 92, 101]. 
Normalising to sums of cross sections removes many experimental systematic uncertainties and the frame 
dependence associated with the flux factors $frame- However, a small spin-averaged cross section in the 
denominator may enhance theoretical uncertainties or hide unfeasibly small count rates: cross-section 
differences A set the scale for the beamtime necessary to perform these experiments. The asymmetries 
^"T vanish in the static limit only for a proton target, but are nonzero for the neutron. 

The relations of the double-polarisation observables (with arbitrary angle between polarisation and 
scattering planes) to the amplitudes Ai are compiled in Ref. [92, Chapter 4.1] for complex amplitudes. 
For real amplitudes, i.e. below threshold, the relations for A™^ and E^^'jf were first reported in Ref. [98, 
Chapter IV. 2]. Those for the helicity amplitudes are compiled for the deuteron in [93, 101] and for 
^He in [70, 71]. Similarly, double-polarisation observables with linearly polarised photons can also be 
defined — see e.g. [88, 92, 93] for definitions and figures analogous to the ones above. 

3 Experimental overview 

In this section, we review the experimental efforts on Compton scattering using proton and deuteron 
targets, spanning the past half-century. For the proton, we have divided the discussion into low-energy 
measurements (below pion threshold) and high-energy measurements (above pion threshold). There is 
a relatively clear distinction between experiments in these two energy regions, although some cases do 
overlap with both regions. There are also some polarised measurements on the proton in the modern era. 



^In Ref. [97], A^'^'^ is denoted by 2Ai^. 

^T,f^^ is denoted by T,2x,2z in Ref. [88] and S^,^ in Refs. [97, 101]. 
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For the deuteron, there are again two categories — elastic Compton-scattering experiments and quasi- 
free measurements in which deuteron breakup is exploited to study the neutron explicitly in quasi-free 
kinematics (with the proton acting as a spectator). 

Special emphasis will be placed on enumerating the statistical and systematic uncertainties of the 
proton and deuteron experiments, since these issues figure prominently in the fitting of the various data 
sets using the xEFT formalism. Somewhat cursory details are given about the experiments: many are 
discussed at greater length in the review by Schumacher [10]. 

3.1 Low-energy proton Compton scattering 

The earliest low-energy Compton-scattering experiments on the proton (up to about cuiab ~ 100 MeV) 
were reported in the mid-to-late 1950s by Pugh et al. [103], Oxley [104], Hyman et al. [105], Bernardini 
et al. [106] and Goldansky et al. [107]. These early experiments were not aimed at measuring the 
electromagnetic polarisabilities of the proton, as we think of them today. In fact, these experiments 
were more motivated by their ability to test recently developed dispersion-theory calculations of Gell- 
Mann and Goldberger and others. Nevertheless, it is noteworthy that in some of these early papers 
(most notably Goldansky et al. [107]), attempts were made to extract the proton polarisability. 

These early experiments were pioneering efforts, given the difficulty of working with continuous 
bremsstrahlung photon beams and detector systems with very poor energy resolution. Large Nal 
photon detectors with good energy resolution were not yet available, and photon-tagging facilities were 
still decades in the future. Normalising the photon flux to obtain an absolute cross section is notoriously 
difficult with bremsstrahlung beams, not to mention that the incident continuous bremsstrahlung beam 
itself has no well-defined energy resolution. For example, the experiment of Oxley [104] had a central 
photon energy of 60 MeV, with a full width of 55 MeV — so it is remarkable that this experiment 
produced results for the cross section between 10.6 and 14.7 nb/sr from 70° to 150°, which is generally 
consistent with modern results. 

Directed efforts were made in the experiment of Goldansky et al. [107] in 1960 to determine the 
proton polarisability. The same group continued these efforts many years later in the mid-1970s in 
the experiment of Baranov et al. [108, 109]. While uncertainties in the extracted polarisabilities were 
reduced by a factor of two in the later experiment, the experimental techniques during that period were 
still relatively crude compared to today. For this reason, it is reasonable to consider the data prior to 
1980 to be exploratory in nature, and to not give very much credence to the absolute scale of the cross 
sections from those experiments. 

Almost 20 years passed before the next proton Compton-scattering experiment was attempted. Two 
major developments in experimental techniques emerged in that period to make the new generation of 
Compton measurements considerably more reliable than their predecessors. First, the method of photon 
tagging was introduced, providing both a mono-energetic beam of photons and a means of normalising 
the photon flux by direct counting of the post-bremsstrahlung electrons. This revolutionised many 
photonuclear experiments. Second, new large-volume high-resolution Nal detectors (25.4 cm diameter 
X 25.4 cm long, and even larger) were available with a resolution AE/E ~ 3% for photons of 50-100 
MeV. These two improvements in the beam and the detectors changed the game signiflcantly, paving 
the way for a new era of Compton-scattering experiments designed speciflcally to pin down the proton 
polarisability. 

The major experiments that have contributed in the modern era to the determination of the proton 
polarisability are those of Federspiel et al. [20], Zieger et al. [23], Hallin et al. [21], MacGibbon et al. [22], 
and Olmos de Leon et al. [24] covering the 10-year period between 1991 and 2001. All but two of them 
used tagged-photon beams of energy wiab < 165 MeV — the exceptions are Hallin, who used a continuous 
bremsstrahlung beam with energies up to 289 MeV, and Zieger, who used a bremsstrahlung beam with 
a proton-recoil detection technique. It is worth examining all of these low-energy experiments in more 
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detail to elucidate their relative merits and potential weaknesses. 

The first of the modern experiments was conducted by Federspiel et al. [20] at the tagged-photon 
facility at Illinois. Tagged photons in the energy range ui^b = 32-72 MeV impinged on a liquid hydro- 
gen target; scattered photons were detected at fixed lab angles of 60° and 135° by two large 25.4 cm x 
25.4 cm Nal detectors. This experimental configuration can be considered a "standard" tagged- photon 
experiment, in the sense that the well-characterised beam and the single-arm detector conditions are 
conceptually extremely simple. The data were divided into 8 bins over this energy range, each 4 MeV 
wide. The resulting cross sections had statistical uncertainties that were roughly ±10%. Systematic 
uncertainties due to photon flux (±1%), target thickness (±1%), detector acceptance (±1.4%) and sim- 
ulation uncertainties ±1%) were correlated among all data points and contributed in quadrature to give 
a total systematic error of ±2.2% for this experiment. The data were predominantly limited by statis- 
tical accuracy, due to the low counting statistics in each tagger energy bin, as well as error propagation 
due to subtraction of random coincidences and empty-target background in the data analysis. 

A follow-up experiment was performed by MacGibbon et al. [22] (including the Illinois group) at 
the tagged-photon facility at the Saskatchewan Accelerator Laboratory (SAL). This experiment was 
essentially identical to the lower-energy version at Illinois and used the same large Nal detectors, which 
had been transported to Saskatoon for this purpose. The tagged-photon energy range was wiab = 70-100 
MeV, so it overlapped with the previous measurement, and the detector lab angles were 90° and 135°. 
One clever twist to this experiment, however, was the ability to extend the usable energy range into the 
untagged region (wiab = 100-148 MeV). The normahsation of the untagged data could be linked directly 
to the well-known normalisation of the data in the tagged region using the bremsstrahlung spectrum 
shape, so a reliable extrapolation of the cross section up to the endpoint energy of 148 MeV could be 
achieved. The tagged data were divided into 4 energy bins, each 8 MeV wide, and the untagged data 
were divided into 5 bins, each 10 MeV wide. Here again, the statistical accuracy of the data was the 
limiting factor — the lowest statistical uncertainty was about ±10% and some points had errors close 
to ±20%. Systematic uncertainties came from similar sources as the previous experiment, but were 
slightly larger. Correlated errors such as photon flux (±1-2%), target thickness (±1.2%) and detector 
acceptance (±1.2%), as well as additional rate-dependent corrections (±1-3%), added up to an overall 
systematic error of ±2.9% for the untagged data and a range of ±3-4% for the tagged data. Only 
a small part of that overall error (±0.4% for the tagged data and ±1.0-1.4% for the untagged data) 
was uncorrelated among the data points. In the overlap region near 70 MeV for the 135° data, the 
MacGibbon results are in excellent agreement with those of Federspiel et al. [20]. 

The experiment of Hallin et al. [21] was also performed at SAL and included participants from the 
two experiments discussed above. The Hallin experiment had three major differences: (1) the photon 
beam was a continuous bremsstrahlung beam, (2) the beam energies were mostly above pion threshold 
(endpoint energies of uy^b = 170-298 MeV), and (3) a single, large- volume Nal detector from Boston 
University was used (the BUNI detector, 50 cm diameter x 56 cm long). Since only a single photon 
detector was used, data for each lab angle in the 25° to 135° range of this experiment had to be 
measured sequentially (as opposed to simultaneously). As a bremsstrahlung experiment, the statistical 
uncertainties of the data can be relatively low, depending on how large a region of the bremsstrahlung 
distribution is utilised to generate the cross section. In this case, different regions were used for the 
angular distributions (15 MeV wide) as compared to the excitation functions (5 MeV wide) — hence, the 
statistical errors in the angular distributions tended to be lower, ranging from ±3.5-5.0% at backward 
angles to ±10-18% at forward angles. The major sources of (correlated) systematic uncertainties arose 
from determining the photon flux (±3%) and the target thickness (±2%), giving ±3.6%. Additional 
contributions from, for example, detector efficiency (±0.3%) and solid angle (±0.3%), as well as other 
factors, led to point-to-point systematic errors in the range of ±1-2%. Taken together, these systematic 
errors (correlated ± uncorrelated) amounted to about ±3.7-4.2% for the total systematic uncertainty 
in the final cross sections, which is quite good considering the inherent difficulties of normalising a 



bremsstrahlung experiment. 

One of the more novel experiments was a 180° scattering measurement by Zieger et al. at MAMI [23]. 
In this case, the forward proton (recoihng at 0°) was detected in a magnetic spectrometer, yielding the 
scattered photon cross section in the backward direction (180°), where the result is exclusively sensitive 
to the difference of the electric and magnetic polarisabilities (a^i " (^mi)- The incident photon energies 
were wiab = 98 and 132 MeV, and each of the two data points constituted an energy bin of width 
16 MeV. The forward proton cross sections were determined by comparing the Compton proton and 
electron yields — this enabled the absolute normalisation to be deduced without detailed knowledge 
of the shape or intensity of the incident bremsstrahlung spectrum. The statistical uncertainties for 
the two data points varied widely (±18.4% and ±5.4%), but the systematic uncertainties due to the 
subtraction of target backgrounds and the normalisation procedure were more consistent between the 
two data points (±6.4% and ±4.3%, respectively). This experiment is the only measurement of the 
180° (backward) Compton-scattering cross section for the proton. 

The most comprehensive experiment to date was the one by Olmos de Leon et al. [24] at MAMI, 
covering the tagged-photon energy range wiab = 59-164 MeV and the lab angular range 59°-155° 
simultaneously. This was possible by using the large TAPS array consisting of six separate assemblies 
of hexagonal BaF2 cells. Statistical errors varied from ±5% at the lowest energies to about ±10% at 
the higher energies. Systematic errors from photon flux (±2%) and target thickness (±2%) were added 
in quadrature to give an overall systematic error of ±3% for the entire data set. In addition, random 
systematic uncertainties (point-to-point) due to the geometry of the individual array elements and the 
detector simulation contributed another ±5% to each data point separately. These data overlapped 
with (and extended higher than) the energy regions of the previous data from Federspiel [20] and 
MacGibbon [22], and all three data sets are in excellent agreement in the overlap region. Regarding the 
extraction of proton polarisabilities, it is worth mentioning that the Olmos de Leon data set tends to 
dominate most global analyses due to the large number of data points and their relatively low statistical 
errors. 

A summary of the low-energy proton experiments, from the 1950s up to the present day, is shown 
in the upper portion of Table 3.1. Plots of sample cross sections are also shown in Fig. 3.1. 

3.2 High-energy proton Compton scattering 

Compton-scattering experiments on the proton have also been explored above pion threshold, dating 
back to the 1960s. In the early experiments of DeWire et al. [110], Baranov et al. [Ill, 112], Gray and 
Hanson [113] and Genzel et al. [114], the main intent was to investigate the region of the A(1232) reso- 
nance. In the 1990s, the emphasis of these experiments evolved to become more speciflc — studying the 
E2/M1 ratio and determining the backward spin polarisability, 7^'' , of the proton. These experiments 
had been interpreted in a dispersion-relation framework, but the recent evolution of EFT has made 
these experiments (below about 400 MeV) accessible to this theoretical treatment. 

Most of the experiments have been conducted at MAMI or LEGS — in the main, they have been 
competing efforts and less complementary in nature. There have been some disagreements in the 
published cross sections, and as a result, there have been diverging interpretations, particularly with 
regard to 'Jt^ . We now review the experiments from each laboratory in turn. 

The MAMI experiments covered the energy region of the A resonance and higher, roughly wiab = 
200-500 MeV, using the photon-tagging facility with a typical tagger-channel energy width of about 
2.0-2.5 MeV. For most of the cases, the scattered photon was detected in coincidence with the recoil 
proton — this kinematic redundancy helps tremendously in the reduction of background due to photons 
from the decay of neutral pions which are copiously produced above threshold. In some cases below 
about 300 MeV, however, the CATS large-volume Nal detector (48 cm diameter x 64 cm long) was 
sufficient to cleanly separate the Compton events from the pion-production events due to its excellent 
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Table 3.1: Compton-scattering experiments on the proton. The column for A^ata shows only the number 
of points for each data set up to a cutoff of 400 MeV. The last column indicates the corresponding 
symbols used in the figures. 

intrinsic energy resolution 

In the paper by Hiinger et al. [115], two independent experiments at forward angles were reported — 
this constituted a summary of work already published. In one case, also reported earlier by Peise et 
al. [116], the CATS detector was located at ^lab = 60° (^cm = 75°) and a multi-element array of small 
Nal detectors was located at the corresponding angle of 6'iab = 47° on the other side of the beamline 
to intercept the forward recoil proton. For additional background suppression, a 27r array of BaF2 
crystals was positioned close to the target, directly opposite CATS, serving as a supplemental veto for 
7r° production events. In the reported cross sections, the statistical uncertainties varied from as large 
as ±40% for the very lowest energies, to ±10-15% below the A resonance, and then closer to ±5% at 
and above the resonance peak. Systematic errors have been estimated to be ±4.4% below 350 MeV, 
±6% for 350-380 MeV, and ±8% above 380 MeV. Within these systematic errors, common sources of 
uncertainty for all data points arose from the photon fiux (±2%) and the target thickness (±2%). 

In the other case, also reported earlier by Molinari et al. [117], sets of 12 photon/proton detector 
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pairs were arranged azimuthally around the target, with all photon detectors at a polar angle of ^lab = 
76° (6'cm = 90°) and the corresponding proton detectors at 6'iab = 44° on the opposite side. The photon 
detectors were Pb-glass blocks and the proton detectors were AE/ E plastic scintillator telescopes. The 
good energy resolution of the proton arm and the small acceptance of each proton telescope enabled a 
clean separation of Compton and pion events. Moreover, the Pb-glass blocks are insensitive to neutrons, 
so the background from charged-pion production was reduced. Statistical uncertainties in the data are 
±4-10% below 400 MeV, ±10-25% up to 450 MeV, and then rather large (±50% or more) up to 490 
MeV. Systematic errors have been estimated to be ±4% below 360 MeV and ±5% above that energy 
due to background subtraction of 7r° events. This included the common uncertainties due to the photon 
flux (±2%), target thickness (±2%), and detector efficiency (±3%). 

A broader survey of energies and angles was presented by Wolf et al. [118], which was also identically 
reported in Galler et al. [119]. This experiment used a large-acceptance detector array covering lab angles 
of 30°-150° and photon energies of cuiab = 250-800 MeV. Much of this data set (the upper half, generally) 
is beyond the applicable scale of xEFT. The photon arm of the detector configuration consisted of 10 
large Pb-glass segments, each containing 15 individual Pb-glass blocks. The proton arm had two wire 
chambers and an array of tall plastic scintillator bars used for proton time-of-fiight measurements. Below 
400 MeV, relatively good separation between Compton and pion events could be achieved. However, 
above that energy, considerable overlap required a separate subtraction procedure based on comparing 
in-plane (Compton-scattering plane) and out-of-plane (pion production with a 27 decay) events to isolate 
the Compton events. This method could add a considerable uncertainty to the extracted cross sections 
for the data above 400 MeV. The random errors associated with each data point (which varied widely 
in magnitude across the extensive data set) included both the purely statistical error due to counting 
statistics as well as individual systematic errors due to detector efficiency, geometrical acceptance and 
background subtraction. Beyond that, there were overall scaling systematic errors related to photon 
flux (±2%) and target thickness (±2%) that affected the entire data set as a whole. 

A separate experiment by Wissmann et al. [120] focused on a single backward angle (6'iab = 131°) in 
the energy range wiab = 200-410 MeV. The motivation was to compare free proton Compton scattering 
with quasi-free scattering, as a precursor to a similar quasi-free Compton experiment on the neutron. 
The experimental setup was simple, consisting of a single photon arm using the large-volume CATS 
detector, which could achieve good separation between Compton and pion-production events below 300 
MeV. Above that point, there is some overlap which can be subtracted using a simulation of the pion 
events. The resulting data set had statistical uncertainties as large as ±20-40% at the lowest energies 
and improved to ±8-13% above about 275 MeV (approaching the A peak and beyond). Systematic 
uncertainties due to photon flux, target thickness and detector acceptance were common to all data 
points and were combined in quadrature to be ±3% overall. 

Finally, in a related experiment, Camen et al. [121] obtained data at a similar backward angle (^lab 
= 136°) at energies of wiab = 200-470 MeV in an effort to extract y^ . This was related to a controversy 
based on an earlier extraction by LEGS [122, 123] which disagreed with standard dispersion theory and 
xPT. As expected, the disagreement can be traced to discrepancies in the measured cross sections at 
the two laboratories near the A(1232) peak. The results from LEGS are discussed below. In the MAMI 
experiment, the CATS detector was used in conjunction with a recoil proton detector array consisting 
of 30 liquid scintillator cells located at ^lab = 18° on the other side of the beamline. This experiment 
was similar to that of Wissmann [120], but the redundancy of detecting the coincident recoil proton 
helped discriminate Compton events from pion-production events without relying on a simulation for 
the latter. The statistical errors were typically in the range of ±7% to ±15% for all data points, but 
there was no mention in the published report of any information regarding systematic uncertainties. It 
is likely, however, that they would be similar to the ones stated by Wissmann [120]. 

Although LEGS covers a more limited energy range than MAMI, it offers the unique capability of 
delivering polarised photon beams. This is accomplished by backscattering polarised ultraviolet laser 
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photons from 2.6 GeV electrons in the Brookhaven National Laboratory synchrotron ring. Energies 
of cuiab = 210-333 MeV can be covered using the photon-tagging facility at LEGS — the photon energy 
resolution is about 5 MeV, which is twice the width of the MAMI tagger channels. At LEGS (as at 
MAMI), these higher-energy experiments were performed by detecting the recoil proton in coincidence 
with the scattered photon, allowing a clean separation of Compton-scattering and pion-production 
events. 

The LEGS polarised experiments were performed in a series of three groups of runs [123-125] covering 
scattering angles of 6'cm = 65°-135°. Scattered photons were detected in a large- volume Nal detector 
(48 cm diameter x 48 cm long); a large array of plastic scintillator bars on the opposite side of the 
beamline was used to detect the recoil proton in coincidence. The proton arm also had wire chambers 
located near the target to reconstruct the proton angle more precisely. The kinematic overdetermination 
for Compton events not only served to isolate them from the larger background due to vr" decay photons, 
but it also enabled detector efficiencies to be determined directly from the data (instead of relying on 
simulations). This helped considerably in reducing systematic uncertainties. 

The statistical precision of the LEGS cross sections is excellent in the vicinity of the A(1232) 
peak, and the systematics also seem to be very much under control. The reported cross sections list 
the combination (in quadrature) of statistical and point-to-point systematic errors. Even with these 
combined errors, the uncertainties are only about ±2% to ±5% for data points between 275 and 325 
MeV. Outside that energy range (the highest energy, 333 MeV, and 213-265 MeV), the error bars are 
in the ±(5-10)% range. The overall systematic error in the cross-section scale was evaluated to be ±2% 
based on uncertainties in the target thickness, photon flux, and geometric solid angles. 

With the exception of a much earlier measurement from Frascati [126] at 318 MeV, these are the 
flrst polarisation asymmetry measurements in this energy region for proton Compton scattering. Near 
the A peak, the asymmetries are not small (0.1-0.5) and the combined statistical plus point-to-point 
systematic errors (arising solely from uncertainties in the photon beam polarisation for the two different 
polarisation states) are typically ±(7-25)% for most data points. Owing to the nature of the asymmetry 
measurement, there are no overall systematic errors that do not cancel out completely. 

Regarding the cross sections, there are clearly discrepancies between the LEGS and the MAMI 
data. This is most apparent in Fig. 11 of Ref. [118] at backward angles close to the A peak. It is this 
discrepancy which gives rise to a disagreement in the values extracted for the E2/M1 ratio and 7^'' . 
At the backward angles, both of these quantities display the greatest sensitivity, especially the spin 
polarisability. The same discrepancy is shown in Fig. 2 of Ref. [121], which is in agreement with [118]. 
Data at the highest SAL energy of 286 MeV [21] appear to be intermediate between the two sets. It 
is worth recalling that both the LEGS experiments [123-125] and the MAMI experiments [118, 121] 
used coincidence methods to detect the scattered photon and the recoil proton, greatly reducing their 
susceptibility to backgrounds from pion production. By contrast, the experiment of Ref. [21] used 
bremsstrahlung (untagged photons) and had a single photon arm using the large- volume BUNI detector. 
It is difficult to know the source of the discrepancy between the MAMI and LEGS data sets. 

A summary of the high-energy proton experiments is shown in the lower portion of Table 3.1. Plots 
of sample cross sections are also shown in Fig. 3.2. 

3.3 A critical look at the proton data 

As a precursor to the EFT fits presented in Section 4.4, we briefly provide a comparison of the various 
data sets that have been used in the proton fltting. Unless otherwise stated, energies and angles will 
always refer to the lab frame. We will use "low energy" to mean the region up to 170 MeV (ending just 
at the upper limit of the Olmos de Leon data), "intermediate energy" to mean 170-250 MeV (ending 
where the Wolf data set starts) and "high energy" to mean 250-400 MeV ( "high" is to be understood 
in a low-energy EFT context). We have attempted to judge the consistency of the world data sets by 
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plotting them as a function of energy for a wide range of angles, as seen in Figs. 3.1 and 3.2 . Initially 
we do not add any theory curves. There is not space here to present all the data in sufficient detail 
to appreciate all the points that can be made; more plots can be found in the supplemental material 
available online. 

In the low-energy region (see Fig. 3.1), we note that there are no data below 45° and none above 155° 
(with the exception of the two Zieger 180° points). The only data below 55 MeV are from Federspiel. 
From 55 MeV to around 120 MeV, multiple data sets contribute; in addition to modern experiments 
(Federspiel, Olmos de Leon and MacGibbon), there are also a number of older data sets which look 
compatible with the modern data. From 120 to 170 MeV, the only data are from MacGibbon, Olmos de 
Leon and Hallin (with the exception of some points from Pugh with large error bars). Figure 3.2 shows 
a selection of data over the whole energy range up to 350 MeV (chosen to highlight kinematics where 
data exist from more than one experiment). From 170-210 MeV, Hallin is essentially the only data set, 
with Bernardini contributing three points and two separate MAMI experiments one point each (both 
with very large error bars). Above 210 MeV, much more data are available, both old and new, but the 
agreement is very problematic. Above 250 MeV, the Wolf data set dominates statistically. 
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Figure 3.1: (Colour online) Data for Compton scattering on the proton below cjiab = 175 MeV. Lab 
cross sections in nb/sr are plotted in bins of 10° lab angle as a function of lab photon energy in MeV. For 
greater comparability of data, a phenomenological parameterisation of the angular dependence has been 
used to shift data points to the nominal bin angle, but the effect is only perceptible for the Oxley and 
Goldansky points at 70° and 120°. The black dot at tuiab = indicates the cross section corresponding 
to the Thomson limit. See Table 3.1 for the key to the data symbols. 

It has been recognised that the extensive data set of Olmos de Leon has larger fluctuations than are 
compatible with the quoted statistical errors, and we follow Wissmann [127] in adding a point-to-point 
systematic error of ±5% in quadrature with the statistical error when fitting this set (not shown in 
plots)^. The Olmos de Leon 133° data tend to be low at energies around 100 MeV compared to other 
data, especially MacGibbon. There is an apparent discrepancy between the Baranov 150° and the 
Olmos de Leon 155° data, with the latter lying high. 

Three points from the modern low-energy data are clearly outliers: Federspiel (135°, 44 MeV), 
Olmos de Leon (133°, 108 MeV) and Hallin (141° (cm), 170 MeV). The Bernardini data set consists 
of 7 points in the range 120-200 MeV (in wide energy bins), and sufficiently many of these are low 
compared with the trend of the rest of the data that we do not include this set in our fits. The Oxley 



^For the other proton data sets, we include all systematic uncertainties in a single floating normalisation for each 
respective set. 
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data set is broadly consistent with other sets, but the statistical errors are so tiny as to give it undue 
weight in any fit. Though we continue to plot all of these, we do not use them in fits. 
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Figure 3.2: (Colour online) Selected data for Compton scattering on the proton below 350 MeV lab 
energy. Centre-of-mass cross sections in nb/sr are plotted in bins of 10° cm angle as a function of lab 
photon energy in MeV. See Table 3.1 for the key to the data symbols. 



In the intermediate-energy region above about 200 MeV, the problems with consistency of the data 
are most significant. The lowest-energy Blanpied data are noticeably high compared to any plausible 
smooth interpolation between the low- and high-energy data; the Hallin data partially reinforce this 
trend. The older Gray and Baranov data lie much lower in this region, while agreeing well a little 
higher, and the few MAMI data points in this region are also lower, though with large error bars. That 
said, the 215 MeV Baranov data points at 91.9° and 107.2° (cm) are too low to be plausible. 

In the high-energy region, the data are dominated by the Wolf data, and the other MAMI data 
look compatible with it (although two of the Wissmann points beyond the peak lie decidedly high). A 
number of points are available from older experiments and most are in good agreement, though as is well 
known, the Genzel, Nagashima and Gray data sets all have one or more points which are implausibly 
low. The major problem, however, is the lack of agreement between the Wolf and Blanpied data for 
angles above 100°, where the Blanpied data sit substantially above the Wolf data. At the worst point, 
the A(1232) peak for 6'cm = 115°, the disagreement is almost 25%. But agreement at more forward 
angles such as 6cm = 90° precludes an explanation in terms of an overall normalisation error. Viewed as 
a function of angle for constant energy, both data sets show more structure than one might expect, and 
it is hard to draw a conclusion about which is more internally consistent based on the data alone. The 
older data sets favour Wolf, while the few Hallin points in this region favour or even exceed Blanpied. 
This issue will be revisited later when discussing the proton EFT fits in Section 4.4. 
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3.4 Low-energy deuteron elastic Compton scattering 

Compton-scattering experiments that have sought to investigate the electromagnetic polarisabihties of 
the neutron are considerably fewer in number than those on the proton, primarily due to the lack of free 
neutron targets. Moreover, for the proton, the electric (a^i) and magnetic {I3mi) polarisabihties enter 
at order u"^ in the cross section (where u is the photon energy) due to an interference with the leading 
Thomson amplitude. For a "free" neutron, there is no Thomson term (the neutron is uncharged), so 
the polarisabihties enter at order u'^ and are much harder to determine. 

It is worthwhile to get a historical perspective on the subject of the neutron polarisabihties before 
proceeding with the Compton discussion. The first experiments designed to investigate this subject 
used a totally different technique. Neutron scattering at very low energies {En < 600 keV) in the 
Coulomb field of a high-Z target (such as ^°^Pb) is sensitive to the electric polarisability of the projectile 
neutron. The initial experiments [128, 129] had very large error bars, but the later experiments of 
Schmiedmayer et al. [130, 131] reported a rather precise extraction of the neutron electric polarisability: 
a^l = 12.0 ± 1.5(stat) ± 2.0(syst). These results were challenged by Koester et al. [132, 133], who 
obtained a^l = 0±5, and by Enik et al. [134], who argued that Schmiedmayer [131] had underreported 
the systematic errors of the experiment and extracted a^l = 7-19 using the Schmiedmayer data set. 
A more recent measurement by Laptev et al. [135] also raises questions about the result of Ref. [131]. 
This debate has not been resolved, and here we can only say that the neutron-scattering technique, 
although successful in determining polarisabihties in more weakly bound systems, appears to have led 
to a weak constraint on a^l, at best. 

The failure of these attempts to measure the neutron electric polarisability in a reliable way moti- 
vated alternate efforts to access the neutron polarisabihties (both a^l and /3j\^i). This led to the use 
of Compton scattering on deuterium, the simplest nuclear target containing a neutron, via either the 
quasi-free d(7,7'n)p reaction or the elastic-scattering d(7,7)d reaction. To date, there have been three 
quasi-free experiments — a pioneering experiment at MAMI [64, 65], a later experiment at SAL [66] and 
the most recent one at MAMI [68]. The last is the most precise, but, even so, the value of a^l was 
obtained to no better than 30% when the statistical, systematic and model errors are combined linearly. 
This emphasises the difficulty of the quasi- free measurements. By contrast, for elastic scattering on deu- 
terium, the Thomson term is recovered, so the polarisability extraction is similar (in principle) to the 
proton case, except for the fact that only the sum of the proton and neutron polarisabihties («^i + 0^1 
and Pmi+Pmi) ^^^ t)e unambiguously deduced from the data. However, a major experimental challenge 
must be confronted due to the fact that the deuteron breakup channel is only separated from the elastic 
channel by 2.225 MeV (plus recoil). As a result of this stringent requirement, only four measurements 
of the d(7,7)d reaction have been performed to date. These experiments are reviewed below, followed 
by the quasi-free experiments in Section 3.5. 

The first elastic-scattering experiment was performed at Illinois by Lucas [25], using an identical 
setup to the experiment of Federspiel [20] on the proton. Two large 25.4 cm x 25.4 cm Nal detectors were 
used to detect scattered photons at four lab angles (50°, 75°, 110°, 140°) at an incident photon energy of 
i^iab = 49 MeV, as well as at two scattering angles (60° and 135°) at a photon energy of wiab = 69 MeV. 
To extract the elastic-scattering cross sections, the inelastic contributions to the measured scattering 
spectra were generated using the Impulse Approximation, and then the sum of elastic and inelastic 
contributions was fitted to the individual spectra. The statistical uncertainties for the deduced elastic 
cross sections were quite good (±4-13%). However, one drawback is that these statistics were obtained 
at the cost of summing all of the tagger channels in the entire tagger focal plane for each energy point, 
so that each data point constitutes an average photon-energy bin width of 6.5 or 7.7 MeV, respectively. 
This is a rather wide energy bin, although if the cross section is linear with energy, there are no serious 
complications when comparing the data to theory. Systematic uncertainties were in the range ±3.6- 
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4.0% for the data, primarily due to correlated factors such as photon flux (±1%), detector acceptance 
(±1-2%), target thickness (±2%), possible contamination of the yield due to inelastic contributions 
(±1-2%), and simulation uncertainties (±1%), leading to a correlated systematic uncertainty of ±3%. 
The remaining systematic uncertainties, at the level of ±2.0-2.5%, were due to subtraction of randoms 
and empty-target contributions, giving an overall total of ±3.6-4.0%. 

The experiment of Lundin et al. [27] was conducted at Lund and was very similar in energy and angle 
to the Lucas experiment [25]. Three large Nal detectors were used at two "nominal" photon energies 
i^iah = 55 and 66 MeV) and at three "nominal" lab angles (45°, 125°, 135°). In this case, however, 
many separate runs were taken over an extended period of time, resulting in multiple independent 
data points which more or less overlap. There were 3 runs at 54.6-55.9 MeV and 3 runs at 65.3-67.0 
MeV, where the actual photon energy varied slightly over those ranges. In each of the 6 runs, the 
detector angles also varied slightly from the nominal values, but the measured angles basically form two 
clusters — a forward cluster near 45° and a backward cluster near 130°. In the end, the total data set 
comprised 18 data points, although the overall kinematic coverage is actually rather limited due to the 
overlap of the various points. For this experiment, the inelastic contribution from deuteron breakup 
was not simulated, but rather it was excluded (to a large extent) from the elastic strength by setting 
a tight summing window in the scattered 7-ray spectra. The statistical uncertainties of the extracted 
cross sections are in the range ±(8-24)%, which are roughly double those of Lucas [25]. In addition, 
the energy width of each point was fairly broad (10 MeV wide) due to summing all tagger channels 
in the entire tagger focal plane to increase the statistics. Systematic uncertainties due to photon flux 
(±5%), detector acceptance (±4%), target thickness (±2%), and possible contamination of the yield 
due to inelastic contributions (±3%) were correlated among all the data points, for a total of ±7.5%. In 
addition, background subtraction gave an uncorrelated systematic error of at most ±4%, leading to an 
overall reported systematic uncertainty of ±7-14% for the measured data. We note that the lower limit 
of this range (±7%) is slightly inconsistent with the overall correlated systematic uncertainty stated 
above (±7.5%); we cannot explain this discrepancy, although there are indeed several data points in 
Table I of Ref. [27] which have systematic errors below ±7.5%. 

Both of the previous experiments spanned an energy range which is only moderately sensitive to the 
polarisabilities. The experiment of Hornidge et al. [26] pushed the energy higher, where the sensitivity 
increases. This experiment was performed at SAL using the large-volume BUNI detector. In many 
respects, this experiment was similar to the Hallin experiment on the proton [21], except that Hornidge 
used tagged photons at energies Wiab = 84-105 MeV and measured 5 angles sequentially in the lab 
angular range 35°-150°. The BUNI detector was crucial for separating the elastic-scattering peak from 
the inelastic contributions which are only 2.2 MeV apart. With a scattered photon energy of ~100 
MeV, a large-volume Nal detector with an energy resolution better than AE/E ~ 2% was required to 
achieve this separation. With the use of a tight summing region in the scattering spectra, the BUNI 
resolution ensured that the elastic-scattering contributions could be cleanly determined. The resulting 
statistical uncertainties for the 5 data points are fairly good (±5-10%). One issue, however, as seen 
before, is that the full focal plane was added together to improve the statistics, and, as such, the data 
points each constitute a broad energy bin with a width of 21 MeV. Systematics due to photon flux 
(±1%), sohd angle (±1.6%), detection efficiency (±3.6%), and target thickness (±2.5%) were correlated 
among the data points for a total of ±5%. In addition, uncertainties due to energy calibration (±1-5%) 
give a range of ±(5-7)% for the overall systematic uncertainty of the data. 

A potential issue in the analysis of the deuteron elastic-scattering data, as mentioned above, is 
that the widths of the energy bins themselves are often non-negligible, e.g. the data of Ref. [26] are 
combined into a single photon-energy bin that is 21 MeV wide. One might be concerned about how the 
energy dependence of experimental acceptances and systematic errors affect these data. However, some 
reassurance is provided by the analysis of Ref. [136], which showed that comparing the average cross 
section over a wide experimental bin with the point theoretical cross section at the centre of the bin, 
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as done in theoretical analyses thus far, is entirely reasonable. In fact, the difference between carefully 
averaging over the energy bin and simply taking the cross section at the central energy is less than 
1%, even for the SAL bin widths of 21 MeV, and nearly an order of magnitude smaller for the other 
experiments; see also [137] for more details. 

As in the proton case, the experimental systematic errors have been separated into two pieces: (1) 
point-to-point contributions which have been added in quadrature to the statistical error, and (2) overall 
scaling contributions for the whole data set which have been subsumed into a floating normalisation, 
see Eq. (4.19) in Section 4.4. We have checked that these separate contributions to the systematic error 
combine to give the overall stated systematic error of the data points. In cases where a discrepancy 
exists, as for some of the Lund data [27], the point-to-point contributions have been taken to be zero. 

A summary of all of the deuteron experiments is included in Table 3.2. Overall, the three published 
experiments [25-27] provide a total of 29 data points in essentially 4 energy bins, centred at lab energies 
between 49 and 95 MeV. The angular coverage is rather limited and overlaps between experiments, 
with only 2 to 5 angles per energy bin. To a large extent, this is attributed to the difficulties in 
clearly differentiating between elastic and inelastic (deuteron breakup) Compton events, which are only 
separated by the small deuteron binding energy. Experimental statistical and systematic errors range 
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Table 3.2: Compton-scattering experiments on the deuteron, including elastic and quasi-free scattering. 
The last column indicates the corresponding symbols used in the figures for elastic scattering. The most 
recent Lund data (2010 and 2012) were not yet available to be included in these figures. 

over ±5 to ±24% and ±4 to ±14%, respectively. At the present time, these are the data that have 
been included in global EFT fits aimed at extracting the neutron polarisabilities. We note that this 
limited deuteron database contains roughly 10% of the volume of data available for the proton, where 
the latter extends far into the A(1232) resonance region and generally has substantially better energy 
resolution and angular coverage, as well as lower statistical and systematic uncertainties. 

There is some promise, however, that the database will be expanded. A more recent experimental 
programme has been conducted at Lund [138, 139] in an energy range higher than the earlier Lundin 
data and comparable to that of Hornidge. The tagged energy range was broader than any previous 
experiment (cuiab = 67-116 MeV) and covered four lab angles (60°, 90°, 120° and 150°). A key improve- 
ment, moreover, is that the measurements were performed using three large- volume high- resolution Nal 
detectors simultaneously in a single configuration — the BUNI detector already described, the CATS 
detector which had been transported from MAMI to Lund for this purpose, and the DIANA detector 
(60 cm diameter x 50 cm long) from the University of Kentucky. This enabled three angles to be mea- 
sured simultaneously in a single run period. It also facilitated a cross-check on systematics by enabling 
the measurement of repeat angles with different Nal detectors over the various run periods. These data 



27 



are still being finalised, and so there are no results to be shown at the time of this review. 

3.5 Deuteron quasi-free Compton scattering 

Compton scattering on a "free" neutron is the objective of the quasi-free reaction d(7, 7'n)p in which a 
kinematic configuration is chosen such that the proton is predominantly a spectator, and the scattered 
photon (backward direction) and recoil neutron (forward direction) are detected in coincidence. Since 
this is a breakup process, it has not yet been treated in Effective Field Theory, although it is now feasible 
to do so and work is beginning along these lines. Nevertheless, a brief review of quasi-free experiments 
to date is appropriate, since it so closely relates to the objective of obtaining precise information about 
the neutron polarisability. 

The pioneering experiment was performed by Rose et al. [64, 65] at MAMI using a 130 MeV brems- 
strahlung photon beam. Two 25.4 cm x 35.6 cm Nal detectors were located at 135° and 90° on one 
side of the beamline and four plastic scintillators were located at 22°, 31°, 39° and 49° on the other 
side. The data were summed over a wide energy range, from 80 MeV up to the endpoint. Since the 
endpoint was below pion threshold, all coincidence events were unambiguously from Compton scattering 
(i.e. no contamination from 7r° production). Due to poor statistics and the low energy range of the 
measurement, only an upper limit could be established for the neutron electric polarisability a^l- Using 
a dispersion theory analysis, this experiment gave results of a^l = 10.7lfo^7 ^^^ I^mi — ^■'^-1%'^ ^^^ ^^e 
neutron, where a sum-rule value of a^l + l3}^j{ = 16 has been assumed. We note that this was the 
first experimental determination (1990) of the neutron polarisability via Compton scattering, since it 
precedes the initial elastic-scattering experiment of Lucas (1994) by 4 years. 

At the low energies of the above experiment, there is a sizable model dependence in the polarisability 
values extracted from the dispersion analysis. Levchuk et al. [140, 141] determined that the model 
dependence could be minimised (and the sensitivity to polarisability maximised) by moving to higher 
energies like 200-300 MeV. At SAL, this was undertaken by Kolb et al. [66] using tagged photons at an 
average energy of 247 MeV (ranging over 236-260 MeV). Scattered photons were measured at 135° in 
the large- volume BUNI detector on one side of the beamline and the recoil neutrons were detected in 
an 85-cell liquid-scintillator array centred at 20° on the other side (note that this is much more forward 
than Rose [65]). An important consistency check was performed in this experiment by simultaneously 
measuring the quasi-free scattering reaction on the proton using the same setup and comparing those 
results to the known case of the free proton. For the extraction of the neutron quasi-free cross section 
at these energies, a simulation of the significant 7r° background is absolutely necessary. This is rendered 
more difficult by the smearing of the Compton and pion events due to Fermi motion of the bound 
nucleon in deuterium. In the end, this experiment yielded an independent lower limit for the neutron 
electric polarisability a^l- Combining this lower limit with the upper limit from Rose, the final values 
obtained from a dispersion theory analysis gave a^l = 13.6lg;o ^^^ /3f^{ = 1.6lo!4 for the neutron, based 
on a sum-rule value of a^l + /3]{^{ = 15.2 [141]. 

The most definitive quasi-free scattering experiment was performed at MAMI by Kossert et al. [67, 
68] using the large- volume CATS detector at 136° and a 30-cell liquid-scintillator array at 21° on the 
opposite side of the beamline. The tagged-photon energy range was 200-400 MeV. This experiment was 
very similar to the SAL experiment [66] in design, but it produced 9 data points over this energy range, 
as compared to a single SAL point at 247 MeV. This data set allowed much tighter constraints to be 
placed on the extracted polarisability values for the neutron. The results from the Kossert experiment 
were a^^l = 12.5 ±1. 8 (stat);j:^(syst) ±1.1 (theory) and l3^\ = 2.7 Tl.8(stat)+°:?(syst)Tl.l(theory) with 

the inclusion of the sum-rule condition a^l +(3^1 = 15.2. The model used to obtain these numbers will 
be discussed in Section 5.5, where a critical assessment of the theory error bar will be provided as well. 
In fact, these results for neutron polarisabihties depended on assuming the model value of 'yi^ obtained 
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by LVov et al. [142]. As pointed out by Levchuk et al. [140], the quasi-free cross section in these 
kinematics is rather sensitive to the backward spin polarisabihty of the neutron, jn . It is therefore 
reassuring that the experiment also afforded an independent determination of this quantity. By allowing 
both a^l — (3^j[ and 7^^^ to vary as free parameters, a fit to the quasi-free cross section yielded a value 
of 77^"^ = 58.6 ±4.0 (including the vr^ pole contribution; cf. Section 2.1). This is entirely consistent with 
the model value employed to obtain the numbers quoted above. The addition of 7Tr as a free parameter 
in fitting the cross section did not alter the extracted values of a^l and l3^{ at all, except to slightly 
increase the statistical error bar. Up to the present, this remains the only experimental determination 
of the backward spin polarisabihty of the neutron. 

4 Compton scattering from the nucleon 

4.1 Dispersion relations 

In this section, we review dispersion-relation (DR) calculations which have been used to extract proton 
polarisabilities from Compton-scattering data. For a thorough review, the reader is referred to Ref. [15]. 
A short description of the DR approach can also be found in Ref. [143]. 

To write down DRs, we must construct a complete set of amplitudes for Compton scattering in 
accordance with the dictates of Lorentz covariance, parity and time-reversal invariance. These ampli- 
tudes must be free of kinematical singularities. The basis of L'vov [142, 144, 145] has been employed in 
Refs. [88, 146]. These amplitudes are linear combinations of the A's defined in Eq. (2.1), with coeffi- 
cients which vary with angle. The relationship can be derived from Appendix A of Ref. [88]. The L'vov 
amplitudes, here denoted by Bi, are written as functions of Mandelstam t and the variable u, defined 
as z/ = (s — m)/(4Mn) = cJiab + t/{4:Mj<f). Under crossing symmetry, they obey Bi{—u,t) = Bi{u,t). 
Within the DR approach, the -Bj's are constructed by adding the nucleon and pion-pole contributions to 
an integral over the spectrum of intermediate excitations. Applying Cauchy's theorem (with a suitable 
contour) to the function Biiy'^t) / {y' — z/ — ie) yields 
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Re B,{u^ t) = Bf-i,., t) + ^P [ du'u' ^'^'^'f^ , (4.1) 
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with i^p™" the contributions from the Born (nucleon- and pion-pole) graphs and P indicating a principal- 
value integral. This is an unsubtracted dispersion relation which applies at the chosen value of t if we 
may neglect the contribution from the semi-circle at infinity. B^ can then be reconstructed from the 
Born terms and its imaginary part, Im^i^j, which itself is l/(2i) times the discontinuity across the 
s-channel cut of the Compton process. This imaginary part can thus be obtained from nucleon-pion 
photoproduction amplitudes, and the integral starts at z/q = m^r + {fn'i + t/2)/{2M^), which is the 
ttN threshold in the lab frame — the frame in which we choose to work throughout this sub-section. 
Pion photoproduction is the lowest channel that contributes to the imaginary part, since we neglect 
the process 7p — )■ e"'"e~7p (Delbriick scattering) and the effect of 7p intermediate states, as these 
contributions to ar are suppressed by powers of Oem compared to the hadronic ones. 

The dispersion relations for t = (i.e. 6 = 0) are of particular interest because the optical theorem 
states that the imaginary part of the forward amplitude is equal to a kinematic factor times the total 
cross section for photoabsorption. As discussed in Section 2.1, only two amplitudes {Ai and As) survive 
for forward Compton scattering, and arioj') does not fall off fast enough as w' — )■ 00 for Ai to obey an 
unsubtracted dispersion relation such as Eq. (4.1). To deal with this, we replace Ai{u}\sh, z = 1) with 
the function y4i(ci;iab, 1) — ^i(0, 1) in the argument that led to Eq. (4.1) and obtain the once-subtracted 
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dispersion relation: 

OD 

Re AiKb, 1) = ^(0, 1) + ^P fdu'^^ , (4.2) 

with Wtt = ■'TiTT + "^^/(2Mn). In the case of ^3, however, Regge arguments allow us to anticipate an 
unsubtracted dispersion relation, 

00 
Re A3(wiab, 1) = P / dw UJ 2 ' (4-3) 

where (T1/2 (o"3/2) is the total cross section for a photon-nucleon system of helicity 1/2 (3/2). 

If the integrals in Eqs. (4.2, 4.3) converge, then each side may be expanded as a Taylor series about 
i^iab = 0. Inserting the low-energy expansions of Eq. (2.5) on the left-hand side then yields a variety 
of sum rules for the structure coefficients in those expansions. For example, from Eq. (4.2) we obtain 
the Baldin sum rule quoted in Eq. (2.7). Meanwhile, equating terms of (9(aJiab) in Eq- (4.3) yields 
the famous sum rule of Gerasimov, Drell and Hearn (GDH) [147, 148] relating the square of k^^^ to 
an integral over cri/2 — ^"3/2. The terms of 0{(jjf^^^ yield the Forward Spin Polarisability (or 70) sum 
rule [8, 149], which was also listed in Eq. (2.7). Recently, Pasquini, Drechsel and Pedroni took such 
arguments one order further, obtaining and evaluating the sum rule for the "higher-order forward spin 

polarisability", 70 [150]: 

00 

7o = ^ydc. -, , (4.4) 

which is defined as the second term in the Taylor expansion of the dynamical forward spin polarisability 
70(0;) about w = 0. These three sum rules demonstrate that higher-order terms in the Taylor expansion 
of Ai and A^ about a; = correspond to integrands with additional powers of the excitation energy in 
the denominator. Consequently, higher-order coefficients are increasingly dominated by near-threshold 
physics, and so their dependence on assumptions regarding high-energy pieces of the integral decreases 
markedly. In particular, the GDH sum rule receives approximately 10% of its total value from con- 
tributions above ^/s — Mp = 2 GeV, while less than 4% of the Forward Spin Polarisability sum rule 
comes from -y/s — Mp > 800 MeV [15]. One corollary is that dispersion-relation results for 70 and 
(even more so) for 70 are quite stable against assumptions about the high-energy behaviour. So, for 
instance, Ref. [15] quoted 7q = —1.01 ± 0.08(stat) ± O.lO(syst), which was updated in Ref. [150] to 
7q = —0.90 ± 0.08(stat) ± O.ll(syst). Ref. [151] makes additional model assumptions and found a 
marginally consistent result: 7q = —0.58 ± 0.20. Another corollary is that dispersion relations for Ai 
and A^, upon which multiple subtractions have been made, predict that the energy dependence of these 
functions is driven entirely by chiral (i.e. near-threshold) physics. 

From such analyses we learn that unsubtracted dispersion relations do not hold for the L'vov am- 
plitudes Bi and B2 and that the unsubtracted dispersion relation for the amplitude B3 converges only 
slowly. It is therefore more appropriate to consider once-subtracted dispersion relations at fixed t [146]: 
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Re B,{u, t) = 5,(0, t) + -i^'P f d^V'^^4^ , (4.5) 
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at least in the cases i = 1, 2, 3. (Here the -Bj's are the structure parts, in analogy to the Aj's of Section 2.) 
In order to obtain predictions from Eq. (4.5), it is necessary to provide the function Bi{0,t) as input. 
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A dispersion relation can also be written for Bi{0,t): 
B,iO, t) = 5,(0, 0) + BJ (0, t) - 57 (0, 0) + - y dt'^L^l^ + _ j at'^L^l^ . (4.6) 

4m2 -00 

The possible intermediate states in the t-channel yield cuts along the positive t-axis, but as long as 
t < 4m|-, this is saturated by vrvr intermediate states. The main contribution which needs to be 
calculated to evaluate the function Bi{0,t) is thus 77 — t- tttt — t- NN. In Ref. [146], the tttt — t- NN 
amplitudes of Ref. [152] were employed, and a model was constructed for 77 — )■ vrvr. 

Meanwhile, the s-channel integral in Eq. (4.5) is largely saturated by ttN intermediate states, which 
mitigates the need for modelling of higher-energy physics in photoproduction. The vrN contribution is 
evaluated using a contemporary parameterisation of the 7p — )■ vrN data. Multi-pion intermediate states 
in the s-channel integral are evaluated from data on the inelastic decay channels of vrN resonances, and 
their helicity structure is assumed to be the same as that of the 7p — )■ vrN amplitudes. This introduces 
some model dependence in the DR evaluation, and there is also model dependence in the computation 
of the second, w-channel, integral in Eq. (4.6), since it requires the evaluation of amplitudes well into 
the unphysical region. However, the overall uncertainty in the evaluation of the DRs is very small since 
"...subtracted dispersion relations are essentially saturated aX u = 0.4 GeV." [15]. Indeed, neglecting 
both the M-channel integral and the two-pion contributions changes the cross sections by only 3-5% [15]. 

The six subtraction constants hi = Bi{0, 0) then remain as the only free parameters in this approach. 
These are equal to linear combinations of the six static dipole polarisabilities. In the results quoted 
below, 64-66 are assumed to be derivable from unsubtracted forward dispersion relations, i.e. to obey: 

00 
^^^2'r^,lu^MtyA ,ith z = 4,5,6. (4.7) 

IT J U' 

where 63 is often fixed from the Baldin sum rule. Therefore, at least three of the six static polarisabilities 
are assumed to be calculable from analyticity arguments, and only aEi, Pmi and one spin polarisability 
(typically chosen to be 7,^) are taken as free parameters in the fit. The formalism does not, however, 
demand that Eq. (4.7) be invoked. All six 6j's could be taken as fit parameters, if sufficient data were 
available to determine them. 

If one does not wish to take these parameters as input, one can continue to work with an unsubtracted 
dispersion relation instead, provided the Cauchy-theorem contour used in their construction is closed at 
a finite radius, rather than at i^' = 00 [142]. This produces a third contribution on the right-hand side of 
Eq. (4.1), Bf^, which is a piece resulting from integrating along a finite semi-circle of radius z/max in the 
complex plane. In Ref. [142] and subsequent works (see Ref. [153] for a review and Refs. [151, 154-156] 
for more recent developments), this contribution is parameterised by the exchange of a single particle in 
the t-channel, e.g. the a meson in the case of Bi. Proceeding in this manner means that one is modelling 
the high-energy behaviour of the Compton amplitudes, and so a significant set of additional dynamical 
assumptions is introduced. The advantage of once-subtracted fixed-t dispersion relations in which the 
subtraction function is fixed via Eq. (4.6) is that they contain minimal dynamical assumptions and are 
sensitive mainly to the well-constrained low-energy dynamics of the vrN system. 

Such dispersion relations do, however, require extrapolation of the imaginary parts into an unphysical 
region. In practice, this extrapolation is difficult to control for large values of t. Hyperbolic DRs provide 
a related approach, which, however, does not necessitate such an extrapolation. In hyperbolic DRs, 
the integration paths are hyperbolae in the Mandelstam plane which correspond to a fixed value of 
cm (or lab) scattering angle. Thus they are often referred to as "fixed-angle dispersion relations". 
The integrals in these calculations have a form similar to that found in subtracted fixed-t DRs, but 
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additional kinematical factors are now present in the integrand. These factors cure the difficuhies with 
the extrapolation of the integrand at large t, but lead to some problems with convergence at small 
angles. The two types of DRs are complementary: fixed-t DRs are good at small angles, while fixed- 
angle dispersion relations work best at ^lab = 180°. We note that such backward DRs have been used 
in their unsubtracted form to directly obtain sum rules for a^l — l3^j[ [157] and 7!^ [158]. Results 
for these quantities from other analogous evaluations at 6'iab = 180° are given in Table 4.2 below. The 
values found for other angles within the domain of validity of the hyperbolic dispersion relations differ 
by ~ 10% for a^l — /3^j[ and 7^'' and by a few percent for other spin polarisabilities. 

Reassuringly, if the same photoproduction data are used as input for both the hyperbolic and fixed-t 
dispersion relations, then very good agreement is obtained for the shape of the low-energy 7p data (see 
Fig. 4.1). (The predictions for hyperbolic dispersion relations are only shown for 6'iab > 107°, as formally 
they do not converge forward of this point.) Within their expected range of validity, it is gratifying that 
all approaches agree at a level which is in accord with the size of the experimental error bars. However, 
we caution that when these curves are continued into the resonance region, more significant differences 
are seen [15, 159]. 
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Figure 4.1: Differential cross section for 7p scattering as a function of cuiab at fixed lab an- 
gles. Experimental data are represented by the symbols, as indicated. Curves are for fixed- 
t subtracted/unsubtracted DRs (solid/dashed) and for fixed-angle subtracted/unsubtracted DRs 
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Figure modified from Fig. 20 in Ref. [15]. Used with permission. 



Regardless, the good agreement in the vicinity of the pion threshold lends strong support to the 
extraction of proton polarisabilities using dispersion relations. Using fixed-t, subtracted dispersion 
relations, the results of a fit to the data of Refs. [20, 22-24] for a)^{ + [5^\ and a^l — j3^\ are shown 
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in Table 4.2 in Section 4.3. The fit also determines 'jn = 10.4 ± 1.8 ± 3.2 and has a x^/d.o.f. = 1.2. 
In this global analysis, the normalisation of each data set was allowed to float (see Eq. (4.19) below) 
within an assigned uncertainty of 3%. The first error bar is obtained by fitting with fixed values 
of the normalisation constants, while the second represents the systematic uncertainty due to these 
normalisation effects. 

4.2 Chiral EFT for 7p and ^n scattering 

4.2.1 Principles of xEFT 

There are many excellent reviews of chiral EFTs, and we will give only a summary here. For those 
interested in learning more, the most comprehensive introduction to the subject are the lectures by 
Scherer and his subsequent book with Schindler [160, 161]. Recent reviews include Bijnens [162] for the 
meson sector and Bernard [163] for the baryon sector; see also references therein. The classic review of 
Bernard et al. [98] is still quite extensively used. Reviews covering the few-nucleon sector include those 
of Beane et al. [164], Bedaque and van Kolck [165], Epelbaum et al. [166], and Machleidt and Entem 
[167]. Phillips' review specifically covers the applications to electromagnetic reactions on nucleons and 
nuclei [168]. 

We start by briefly reviewing xPT, in which the dominant dynamics is directly linked to the chiral 
symmetry of QCD. Subsequently, we will introduce the A(1232) as a dynamical degree of freedom in a 
more general chiral EFT. 

Pionic xPT is a mature subject with a voluminous literature, but most of the details are scarcely 
relevant here. The Lagrangian initially includes all allowed terms with either two derivatives or one 
factor of m"^, reflecting the fact that soft pions are non- interacting in the chiral limit. The resulting 
C)t — see Eq. (4.8) — contains one term of each type and generates the free-pion propagator away from 

(2) 

the chiral limit. Writing Cn iii- a. chirally symmetric way generates non-vanishing amplitudes for 4, 
6, . . . pions interacting at a point, which are predictive at this order. These reproduce Weinberg's 
low-energy theorem for tttt scattering [169] and give lowest-order predictions for the scattering lengths 
which are within 25% of currently accepted values. At next-to-leading order (NLO, (9(P^)) there are 
contributions both from terms in the fourth-order Lagrangian C^^^ (with, e.g., four derivatives) and also 
from one-loop diagrams such as a single rescattering; at next-to-next-to leading order (N^LO) there are 
contributions from C^^\ from two-loop diagrams and also from one- loop diagrams in which one of the 
vertices is taken from C^"^^ [170]. Of course, beyond C^'^^ unknown LECs enter, and various strategies 
exist. One can fit to data and excellent agreement with the scattering amplitudes near threshold is 
obtained [162, 171, 172]. Alternatively, lattice QCD can now constrain the fourth-order LECs [60-62]. 
Whichever approach is taken, tttt scattering is a textbook example of a systematically improvable EFT 
calculation. 

A chiral Lagrangian with nucleon fields coupled to pions and photons was first investigated beyond 
tree level by Gasser et al. [173]. Since the Dirac equation is linear in derivatives, terms with both odd 
and even numbers of derivatives appear. It was immediately recognised that the simplest possible loop 
contribution to the self-energy, from one ttN loop using dimensional regularisation, does not obey the 
power counting found in the pionic case; rather than being suppressed by some specific power of m.,r/A^ 
relative to the leading nucleon mass, it generates an entire series of powers of mT^/M-^, including a 
contribution ~ Mj^. This effectively precludes the possibility of a systematic calculation in which all 
omitted diagrams contribute at a higher order in m^/A^ than the desired order. 

In subsequent years, methods to tame this problem and restore power counting to the theory with 
Dirac nucleons were developed. The first solution which was proposed, and the one which has been 
used in most work on Compton scattering, is called Heavy-Baryon xPT (HB^PT) [98, 174, 175]. Other 
approaches will be mentioned in Section 4.2.7. In the framework of HB^PT, the nucleon mass is 
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recognised to be of the same magnitude as A^, and so the free Dirac Lagrangian contains terms which 
contribute at different chiral orders in low-energy processes: specifically the mass and energy are zeroth 
order while the momentum is first order (i.e. of order mj^). By decomposing the nucleon field into 
"large" and "small" components (upper and lower components in the nucleon rest frame) and by 
integrating out the "small" component, one is left with a theory of Pauli spinors. As guaranteed by the 
relativistic starting point, Lorentz invariance is maintained in the new Lagrangian, but perturbatively, 
with corrections of order (P/Mn)"'*'^ when working to nth order. The mass and its contribution to 
the energy cancel and disappear from the leading-order Lagrangian C^^ — see Eq. (4.9). From C^^ 
onward, terms appear which have been generated by the heavy-baryon reduction, and so terms with a 
fixed coefficient proportional to 1/M^~^ at nth order occur, in addition to those with an undetermined 
LEG. In this theory, corrections to the bare nucleon mass arise as expected at second order, with a 
contribution proportional to m^ (or equivalently mg) with the LEG Ci as its coefficient (this also gives 
the leading contribution to the pion-nucleon sigma term — see Eq. (4.10)). Then, at third order, the 
loop diagram gives a finite contribution proportional to m^ only. (The contribution needs to be finite, 
because terms in the Lagrangian can only involve integer powers of rrig or equivalently, even powers of 
mTT- Thus there can be no LEG to cancel a divergence at m^.) HB^PT does have a transparent power 
counting which matches that in the mesonic theory. The Lagrangian to third order was developed in 
Refs. [176-179], and then to fourth order in Refs. [180, 181]. 

One comment on power counting and electromagnetic processes: since a single Lagrangian term 
with a gauged derivative can contribute either a power of momentum or a factor of |e| = y/AiraEM, and 
since the magnitude of |e| and of m.,r/A^ are comparable, they will be considered interchangeable in the 
perturbative expansion. This is just a matter of bookkeeping provided the only photons in the process 
are external ones: Gompton scattering from the nucleon, for instance, will start with the Thomson term 
at 0{P^), and working to N^LO will involve terms from the fourth-order Lagrangian. Photon loops can 
also be included; their effects will generally be of the same size as strong isospin-splitting effects at the 
appropriate order, and both have so far been ignored in Gompton scattering. 

4.2.2 The Lagrangian for Gompton scattering in xEFT 

Here we discuss the relevant terms in the Lagrangian for the construction of the Gompton scattering 
amplitude to fourth order in HB^PT. The full Lagrangian needs to be written in terms of building 
blocks with appropriate chiral properties, and hence every term can give rise to interactions with 
multiple pions. The usual notation is compact but far from transparent. Below we retain only the 
relevant structures for our purposes: 

£1^) =1 d^cj> ■ d>^(t> + eA^ 63,, (j)^d>^<j), + \ e^A^A^{4>\ + 0^) _ l^l^^ + . . . (4.8) 

6ll=i,^{:wD^g^u-S)i, (4.9) 

+ (c2-^)(t;-«)^ + C3«-«-^[^^^1eF,.((l + «:(^)) + (l + «:W)r3)|^ + ... (4.10) 

(4.11) 

where ip is the nucleon field, 0a are the pion fields, F^'^ is the electromagnetic field tensor and D^ = 
d^ — ieQA^^ is the gauged derivative; v^ = g^^ and S^ = (0, ff/2) in the rest frame of the nucleon. Since 
this is meant to be an expansion around the chiral limit of QGD, the LEGs written here as g^, fn, ^ etc. 
are actually "bare" quantities which will differ from the physical quantities by loop corrections which 
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only vanish in the chiral hmit. This distinction is relevant at 0{P^). The object n^ is given by 

M^ = — ^M^0a + ee'^=^V„</.bA^ + ...) . (4.12) 

Jit 

For the physical values, we use mj = 139.6 MeV, /^ = 92.42 MeV, Mn = Mp = 938.3 MeV, Qa = 1.267, 
K^*^-* = —0.22 and k^^^ = 3.71. The neutral pion mass, m^o = 134.98 MeV, is typically used in the vr*^ 
pole diagram but not elsewhere; other isospin-breaking effects are neglected. The second-order LECs 
are less well known. They have been determined by fits to both vrN and NN scattering (since two-pion 
exchange is a significant part of the NN force); we take the values of Bernard [163] ci = — 0.9to'5, 
C2 = 3.3 ± 0.2, C3 = -4.7tJ:^, all in GeV-^ 

4.2.3 Born terms 

We now discuss the HB^PT nucleon Compton amplitude in xPT. First, we focus on the Born terms. 
The low-energy theorems of Compton scattering provide a simple example of how the same physics plays 
out in the Dirac and heavy-baryon pictures; we refer to Eqs. (2.3) and (2.5). In the Dirac picture, the 
Born (non-structure) terms come from the two diagrams with an intermediate nucleon, with direct and 
crossed photons (Fig. 4.2(a)). If the resulting amplitude is expanded in powers of 1/Mn, the Thomson 
term is reproduced in the amplitude Ai; at one order higher, when the anomalous magnetic moment is 
included in the photon coupling, the Born terms in A^^q appear. 



(a) ^^^ ^ / (b) ! (c)~*~ """" /^ 

Figure 4.2: (Colour online) (a) Nucleon Born diagrams for Dirac nucleons; (b) Pion Born diagram; (c) 
Born diagrams in HB^PTup to 0{P^): solid (red) dots indicate vertices from C^.^ and open sliced dots 
from £^N- ^^ ^^^ counting schemes to be introduced in Sections 4.2.6 and 4.2.7, the first diagram of (c) 
counts as e^ and e^, and (b) and the rest of (c) count as e^ and e^(5^. 

In HBxPT, however, the leading (electric) photon-nucleon coupling is sufficiently simple that the 
direct and crossed diagrams just cancel. Indeed, for real photons it is possible to work in a gauge in 
which this vertex is actually absent (transverse or radiation gauge, corresponding to a purely space-like 
photon polarisation vector k ■ e = in the nucleon's rest frame). The HB^PT diagrams which do 
contribute to the Born terms of Eq. (2.5) up to 0{P^) are shown in Fig. 4.2(c). In transverse gauge, 

(2) 

the leading 7NN vertex comes from C)^-^ and includes a magnetic coupling proportional to the sum of 
the Dirac and anomalous magnetic moments. Thus, the two right-most diagrams in Fig. 4.2(c) give 
a third-order contribution proportional to 1/M^, accounting for the (Q + k)^ Born terms in ^3 to A^, 
and also for the Born term in Aq. However, the missing terms in Ai and ^3 come from photon-nucleon 
seagulls at second and third order, respectively. In particular, the Thomson term comes from gauging 
the leading kinetic term in the Lagrangian; as discussed above, it counts as second order even though 
it has no powers of chiral momenta because it is proportional to Oem- This demonstrates that the 
distinction between Born and non-Born in HB^PT cannot be equated with (apparent) one-particle 
reducibility or the lack thereof. It is always safest to define Born terms via Dirac nucleons. 

Of course, strictly speaking, in all these Born diagrams, it is the chiral-limit magnetic moment, k^^\ 
which appears, not the experimental value. The difference between n and n^^^ is (!?(?7i^), and so this 
affects the amplitude only at fourth order. This issue will, however, recur below. The vr'^-pole diagram 
Fig. 4.2(b) also contributes to amplitudes A3 to Aq at third order; its contribution is denoted by Aj 
in Eq. (2.6). 
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4.2.4 Leading-order structure contributions and polarisabilities 

It is gratifying tliat tlie basic low-energy tlieorems of Eq. (2.3) are reproduced in tliis EFT, but our 
interest is in the predictions made by the theory for the structure-dependent amplitudes, including the 
static polarisabilities a^i, (3mi and the 7's. As just described, the leading-order HB^PT Compton- 
scattering amplitude is simply the Thomson term. At NLO — 0{P^) — there are the spin-dependent 
Born contributions described above, but there are also contributions from pion loops [182], specifically 
the diagrams depicted in Fig. 4.3. Individually these diagrams are divergent and violate the LETs, but 



Figure 4.3: (Colour online) 0{P^) loop diagrams in HB^PT; all orderings of vertices and crossed as 
well as direct photons are implied. Vertices (shown without dots) are all from the LO Lagrangian, that 
is, £^j^ fo^ the nucleonic coupling and Ci for the 77r couplings. These also count as e^ and e^(5^. 



the sum is finite and leaves the Born contributions intact. Thus the sum of the loop diagrams contributes 
only to the structure parts of the six amplitudes and hence vanishes quadratically for Ai and A2 as 
cj — )■ and as the third power of u for A3_q. The coefficients of these terms are the polarisabilities, 
and at this order they are the same for both the proton and neutron. The results, first calculated by 
Bernard et al. [98, 182], are 

Uei = WPmi = j^ 7^ = 12.5 , 'jEiEi = 57M1M1 = -57A/1E2 = -o'jEiAn = - Q„ 2 2 f 2 = ~^-"- 

(4.13) 
It should be stressed that up to third order the full amplitudes, as well as the polarisabilities, are 
entirely predicted in terms of the well-known quantities m.,r, f^^ and Qa', there are no free parameters. Of 
course, the best method to analyse experiments for extracting even aEi and /3mi is the subject of this 
review, but nonetheless, the many attempts made in the past to measure these quantities all come out 
close to these values for both the proton and neutron; in particular, the order-of-magnitude difference 
between aEi and Pmi and their nearly isoscalar nature is not easily understood in most models. This 
has long been lauded as a stunning early success of HB^PT. (As the spin polarisabilities are less well 
known, it is harder to judge these predictions; see Section 4.3.) 

There are a number of caveats, however. Even strictly within HB^PT, one would expect higher- 
order corrections to be of order -P/A^ — around 20% if the scale of the expansion were A^^ ~ nip. There 
is also good reason to expect that for Pmi (as well as •Jmimi), the scale is actually set by the much 
smaller A-nucleon mass difference M^ — M-^. Furthermore, in a relativistic framework, the predictions 
from the diagrams in Fig. 4.3 are substantially smaller: a^l = 6.8, (3^j{ = —1.8 [183, 184]. But, before 
dismissing the success of third-order HB^PT as a fluke, we should step back and remember that the 
calculation gives us full amplitudes as a function of u, not merely the static polarisabilities. As will 
be shown in more detail subsequently, the full third-order cross section extends the region in which 
data can be well described substantially beyond that where the Petrun'kin cross section (Born plus 
static scalar polarisabilities) is valid. In particular, it reproduces the pronounced cusp at the photopion 
threshold which is seen at forward scattering angles (see Fig. 3.1). Beyond that point, the data show a 
huge rise in the cross section which is obviously due to the A(1232) (see Fig. 3.2), and one could not 
expect a theory without the A(1232) to work in that region. 

For completeness, we should mention that a handful of calculations of polarisabilities have been 
done in the framework of SU(3)xSU(3) chiral perturbation theory, involving kaons as well as pions and 
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all the octet baryons. Bernard et al. calculated the spin-independent static polarisabilities in HB^PT 
[185] and showed that for nucleons the effect of kaon loops was small (see also Butler and Savage [186]); 
Vijaya Kumar et al. found a similar result for 70 [187]. Dynamical polarisabilities 0^1(0;) and [3mi{(^) 
have also been calculated at NLO in a covariant framework by Aleksejevs and Barkanova [188]. 

4.2.5 Structure beyond leading order 

Although the ability of third-order HB^PT to qualitatively describe low-energy data is encouraging, 
the lack of any free parameters limits its use as a tool to extract more information from those data. 
This situation changes at fourth order, because at that order we can construct Lagrangian terms like 
ip'^F^'^F^j.il) which are multiplied by new, undetermined LECs. Such terms give rise to photon-nucleon 
seagull diagrams which contribute terms proportional to uP' to the amplitudes Ai and A2 [189]. In the 
enumeration of Ref. [181], there are actually six such terms (numbers 89-94) but in the photon-nucleon 
sector only four independent combinations of LECs enter, which we can call 5a^{, Sa^l, 5f5^\ and 515^\ 
(see C\{^, Eq. (4.11), and Fig. 4.4). These are contributions to the spin- independent polarisabilities 
of the proton and neutron which come from non-chiral physics — for example, quark substructure, or 
resonances, according to perspective, and they obviously encode the leading effects of a A(1232) pole. 
In addition, at fourth order a new set of vrN diagrams has to be included. Finally, all the N^LO terms 
in the expansion of the relativistic Born contributions to Ai and A2 are also generated via fourth-order 

(1') 

seagulls and diagrams like those of Fig. 4.4 with either one vertex taken from C\^ or with an NLO 
nucleon propagator. 



^ ^ ^y./ ,/ ,/ ,/ 



I a L _L 



X^\ 



Figure 4.4: (Colour online) O(P^) diagrams in HB^PT; vertices labelled as in Figs 4.2 and 4.3 with 
the addition of a (magenta) diced dot for the fourth-order counterterms Sa^l etc. of C)^-^- All orderings 
of vertices and crossed as well as direct photons are implied. Omitted are all diagrams obtained from 
those in Fig. 4.3 by substituting an NLO vertex or propagator for an LO one. These also count as e^ 
and e^(5^, though the final diagram is included at one order lower if polarisabilities are fit. 



Of the loop diagrams, many are 1/Mn corrections to the diagrams of Fig. 4.3 (no new LECs enter 
in these). However, there are also two new types of diagrams — those with magnetic- moment couplings 
as well as a pion loop, and those with a pion-nucleon seagull, as shown in Fig. 4.4. In the former, the 
only new LECs are the well-known proton and neutron anomalous magnetic moments. In the latter, 
however, the three vrN LECs Ci, C2 and C3 enter. 

This time, the sum of all loop diagrams does make a 0{uj) contribution to the Born terms. The 
contributions are exactly those which are needed to replace the chiral-limit k^^' with the correction that 
shifts K to its experimental value at this order. In the expansion of the 7N vertex, this shift comes 
from a diagram in which the photon couples to a pion loop, as in the third diagram of Fig. 4.4, but in 
Compton scattering this is not the only diagram which gives 5k corrections to the Born term, nor is such 
a correction the only contribution from this diagram [190-192]. The 0{uj'^) piece of the sum of all fourth- 
order loop diagrams produces a logarithmically divergent result for the spin-independent polarisabilities. 
These divergences are cancelled by the divergent parts of 5a§{ etc. to leave a finite but undetermined 
total fourth-order contribution to the spin-independent polarisabilities [189]. By contrast, the 0{uj^) 
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parts of all these loop diagrams yield a finite shift for the spin polarisabilities [190-194]; in particular 
the contribution of the one-nucleon-reducible diagrams in Fig. 4.4 may not be omitted [194, 195]. 

There are actually a number of subtleties in the fourth-order calculation which are not present at 
third order, one of which is the frame dependence. Since the kinetic energy p'^/[2M^) of a nucleon is 
one chiral order higher than that of a photon (which is the same as its three momentum and assumed 
to be of order ttt,^), at lowest order the incoming and outgoing photon momenta are the same, and 
the same in any frame in which the heavy-baryon reduction is valid. That includes the lab, Breit and 
centre-of-mass (cm) frames. Thus, in the lowest-order {0{P^)) loop calculation, we do not need to 
specify in which frame we are working. At (9(P'^), however, that is no longer the case. The Breit frame 
is the simplest in which to calculate and has the additional merit over the commonly used cm frame 
that the amplitudes are crossing-symmetric. (Already at third order the cm frame violates this: the 
Born terms, but not the loops, are different in the two frames.) Furthermore, in identifying the terms 
a^/cu^ etc. in the amplitudes, it is desirable that this is as close as possible to what is used in other 
determinations of the polarisabilities. In the Petrun'kin cross section, it is wiabi^iab which multiplies the 
polarisabilities; in DR analyses, it is z/^. These differ from one another and from w|j.cit ^Y terms of order 
tu^/M^ (which matters only at 0{P^)), whereas u"^^ differs at 0{w/M^). Having obtained the Breit- 
frame amplitudes, those in other frames can be obtained from a boost, expanded to the appropriate 
order in 1/Mn. Whereas at 0{P^) the cross section differs noticeably depending on which frame is 
chosen, at 0{P'^) this dependence is negligible. (This frame dependence at 0{P^) can be confusing, 
given the lack of frame dependence in the form of the amplitudes. But the variables in the two frames, 
(cUcm, 6*011) and (wBrcit, ^'Brcit), a^G not the same as functions of the lab variables, see Section 2.3.) 

A more physical issue than the choice of frame, passed over in the discussion of third order above, is 
the fact that any amplitude calculated to finite order will put the photoproduction threshold at w = m^ 
(irrespective of the frame in which u is defined). This is because the incoming nucleon kinetic energy is 
included perturbatively and is not present in the nucleon propagator in the vrN loop. Of course, within 
the heavy-baryon amplitude there will be a string of terms which would be generated by the 1/Mn 
expansion of a loop in which the nucleon is allowed to recoil, and, as noted by Bernard et al. [196], the 
solution is to resum these terms to put the threshold at the correct place. In practice, that involves 
replacing u within loop integrals by a new variable Usiuj), which differs from u by terms of 0{u/M-^) 
and which equals ?7i^ when u = Uth, and then writing a Taylor expansion in powers of Ug/M^ and 
retaining only those terms required at the order to which one is working. At third order, no expansion 
is required, and there is simply a substitution of variable; that has always been done before comparison 
with experiment. At fourth order, the third-order loop pieces do need to be expanded, and indeed only 
with such a prescription are the amplitudes finite at threshold [196]. Ug is not uniquely determined 
by this prescription, but it is commonly taken to be Ug = ^/i — M^. A different choice was made in 
Refs. [196-199]. 

In Ref. [197], McGovern took the values of a^l and (3^j[ from the Particle Data Group (which were 
very close to the third-order values) and made a comparison with all extant data below 180 MeV. The 
resulting conclusion was that there was a modest improvement over the third order up to around the 
photoproduction threshold, but no improvement beyond that point, with particularly poor agreement 
at backward angles. Of course, there is still no A(1232), but the tadpole diagrams of Fig. 4.4 would be 
obtained from a theory in which an explicit A(1232) had been integrated out to leave the pion seagull 
terms proportional to the LECs C2 and C3, so it was somewhat surprising that they did not extend 
the fit a little further. Of course, the biggest contribution from an explicit A(1232) would be the pole 
diagram, and it is clear from comparison to the data that by 200 MeV, replacing this with one term 
of a polynomial expansion (6/3^{u'^ etc.) is quite inadequate (see Fig. 4.6). In Ref. [198, 199], Beane 
et al. instead fitted a^l and l3^{ to the data in the same region (with a cut on |t| as well as u). As 
will be discussed in Section 4.4, there are problems with the data set which precluded a fit with an 
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acceptable x^. Nevertheless, a fairly robust conclusion emerged in which a higher value of /3j^^^ than 
that generally accepted was required by the data — the central value quoted was 3.4. If the Baldin sum 
rule was imposed, this value was decreased to 2.8 (see Table 4.2 in Section 4.3 for detailed results). 

4.2.6 Including the A(1232) in the small-scale expansion 

The A(1232) resonance has long been recognised as hugely important in the physics of nucleons. In 
xPT it is not explicitly present, but its influence is felt through LECs such as 02,3 and 5f3^i\. But, as 
the radius of convergence of an EFT is set by the scale of the lowest degree of freedom which has not 
been included, the A(1232) can be expected to severely restrict the applicability of xPT at least in 
those processes in which it contributes, with the convergence governed by the scale Ajv/ = Ma — M^. 
And any glance at Compton-scattering data above 200 MeV, as in Fig. 3.2, confirms that this is such 
a process. 

In fact, since Am ^ 2mn, it could be argued that these two scales are similar and should be included 
on the same footing. This is the basis of the so-called "small-scale" or "e" expansion developed by 
Hemmert et al. [200, 201] (see also the earlier work of Butler et al. [186, 202]). Exphcit A(1232) fields 
are included in the Lagrangian, as first shown by Manohar and Jenkins [175], and Am/A^ is counted 
like ptyp/A;^ and m^/A^ in determining the order of a diagram. The heavy-baryon expansion is used, 
ensuring that only Am and not Ma appears. 

The relevant terms in the Lagrangian are the following: 

4'^ =(A^)t (-i^ . D + Am) A'' (4.14) 

4na = - ^(^t^-^^^A^ + (AJ,)ta>V + • . .) (4.15) 

4na = 1^(^^^P^"'^J - (K^S^F^y) (4.16) 

where A^ is the heavy-baryon reduction of an J = |, S = | Rarita-Schwinger field \l/^, with i and /i 
being the indices on the (iso)spin-l vector coupled to the (iso)spin-i spinor. The coupling constant (^ttNA 
is fit to the A(1232) width and varies significantly depending on whether relativistic or non-relativistic 
kinematics are used. The transition magnetic moment 61 also has no single widely accepted value. Both 
will be specified when we use them later. '^ It is also worth displaying here the alternative form of C^j^/^ 
used in the 6 expansion, see later (here ip is the nucleon Dirac spinor and other notation is that of 
Ref. [205]): 

^ = 2M,(M^+M,) (^"*'^'" - 9-^^^"'')«'.*J - ^n'SAm.F'"- - S.7.F'"')*) • (4.17) 

The leading (magnetic) term in the heavy-baryon reduction of this Lagrangian is equivalent to the one 
above with the identification 

^m = &i(1 + Ma/Mn)/3 (4.18) 

(though it is important to note that if the full vertex is used, there are substantial sub- leading terms), 
but there is also a sub-leading electric coupling g^. The ratio of these two couplings (at the A(1232) 
pole) can be obtained from the E2/M1 ratio to be —0.34 [206, 207]. Being third order, the electric 
contribution to the amplitudes is suppressed by a power of u/M^ relative to the magnetic one. In the 
following, we will sometimes use 62 for the electric coupling, defined via qe = ^2(1 + Ma/Mn)/3. 



^There is some confusion in the literature by what is meant by the magnetic couphng bi. The early papers of Hemmert, 
Holstein et al. [200, 201, 203, 204] were not always consistent with one another and did not always specify which Lagrangian 
was being used, so that signs and factors of 2 come and go. Although it is not stated in the paper, Hildebrandt et al. [92] 
use the Lagrangian displayed above, which gives a contribution to Pmi of 2aEM^i/(9Aj\/Af^). 
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Figure 4.5: (Colour online) A(1232) contributions to Compton scattering at e^ and e^5^ (see text for 
explanation). All orderings of vertices and crossed as well as direct photons are implied. 



In the small-scale expansion, A(1232) propagators {v ■ p — Ajv/ + i?7)~^ ~ e^^ ~ fn"^ have the same 
counting as nucleon ones [v-p + irf)'^, and so the diagrams of Fig. 4.5 are also third order (denoted now 
0{e^) rather than 0{P^)). Those vertices of the theory without external As are unchanged, except that 
the LECs C2 and C3 need to be retuned to exclude the A(1232) contribution; this is done by subtracting 
the A-pole contribution to soft vrN scattering. There are significant complications in including spin-| 
fields, not least of which is the need to ensure that unphysical spin-| degrees of freedom do not propagate, 
but at leading order this is not an issue. The contribution of the A(1232) to the polarisabilities in this 
framework was calculated in Refs. [203, 204] (see also [186, 202]) and to the full Compton amplitudes in 
Refs. [90, 92, 203]. Of course, these are dependent on new and rather poorly known coupling constants, 
namely the vrNA coupling constant and the 7NA transition magnetic moment. 

Long before Hemmert's work, though, it was known that the A-pole diagrams give a huge paramag- 
netic (isoscalar) contribution to /3a/i- The exact value depends on the parameter values, but Mukhopad- 
hyay et al. estimated 5(3^i = 7.0 [208], and Hemmert et al. obtained values of 7.2-12 [203, 204]. In 
the latter work and in Ref. [186], it was also shown that the vrA loop diagrams made a substantial con- 
tribution to oei- With no counterterm to adjust at third order, there was no possibility of agreement 
with data. The pragmatic solution, adopted by Hildebrandt et al. in Refs. [14, 90, 92], was to promote 
the fourth-order isoscalar counterterms 6aEi and 6f3Mi to third order, so that these unwanted large 
contributions to the polarisabilities could be cancelled. The prediction for the rest of the amplitudes 
remained intact, of course. 

In these works, the third-order theory with an explicit A(1232) was fit to the Olmos de Leon and 
Hallin data up to 240 MeV (wcm = 200 MeV). In view of the uncertainty in the A-sector parameters 
in this formulation, the transition magnetic- moment strength bi was also fit along with Oei and Pmi- 
[Am = 271.3 MeV was fit to the real part of the A pole position, and g-j^f^A = 1-12 to the imaginary part 
using the non-relativistic formula for the width.) A good description of the data was obtained. The value 
of fei = 4.7 may seem high compared to the value of 3.7 obtained from a recent fit to photoproduction 
data [206, 207], but since the latter used the Lagrangian of Eq. (4.17), no direct comparison can be 
made. As expected, the problem at backward angles that was evident without the A(1232) was cured 
and no restriction on |t| was required. Perhaps surprisingly, the values of aEi and /3mi obtained were 
close to those of Beane et al., although with less tension with the Baldin sum rule. The value of /3j^{ = 
3.4, in particular, is again higher than the traditional value and appears to be a feature of chiral EFT 
fits. It should be noted that though the theory now includes the A(1232), its reach is still limited to 
i^iab ^ 240 MeV, since the power counting suppresses the loop diagrams which need to be resummed to 
give the A a width, and so the amplitudes diverge as u; — )■ A a/. 

4.2.7 The 5 expansion 

As important as the A(1232) is above the photoproduction threshold, its influence on the cross section 
diminishes rapidly as the energy is reduced. Arguably, with the particular values of Aa/(= Ma — Mn) 
and rriTr obtaining in the real world, counting the two as the same scale gives undue prominence to the 
A(1232) in the region where static polarisabilities are important. An alternative counting was proposed 
by Pascalutsa and Phillips [205] in which tti^/Aa/ and Am/A^ are counted as proportional to the same 
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expansion parameter 6 (the so-called "5 expansion"). In this counting, the A-less theory is an expansion 
in powers of 6"^, and, for low energies, the first contributions from the A (the pole diagrams and the vrA 
loops) intercalate between the third and fourth orders of xPT. The main advantage of this expansion, 
however, is that it allows for two separate energy regions, u ~ ttItt and u ~ Am, and, in the latter 
regime, there is no suppression of vrN loop contributions to the A propagator. Thus, these must be 
resummed, and the A becomes an unstable particle with a width F, with the new propagator going as 
{'^ — MA + iT{p^)/2)~^. For cu ~ Am, the one-A reducible diagram dominates (that is, the direct A-pole 
diagram), with all other contributions being sub-leading. In order to have amplitudes which could 
be used over the entire energy region, Pascal utsa and Phillips added just the resummed, direct A-pole 
diagram to the standard third-order xPT amplitudes. This includes all contributions to 0{e'^6^) in each 
region, at the expense of having an incomplete set of higher-order contributions in each.^ This meant 
that the sub-leading electric 7NA coupling was also included. (Two further technical details relevant to 
Ref. [205] are that: (a) the A Lagrangian used was one developed by Pascalutsa [209, 210] which ensures 
freedom from spin-| degrees of freedom via the imposition of an extra local gauge symmetry; and (b) 
the direct A-pole contribution was included with relativistic kinematics and the full Dirac vertex rather 
than being expanded in powers of 1/Mn.) Importantly, in this approach there is no promotion of the 
SaEi etc. terms of Eq. (4.11), and the polarisabilities are predicted in terms of the two 7NA couplings 
Qm and gE, which were fit to Compton data. Because vrA loops are absent, uei is close to its HBxPT 
value (there is a small negative contribution oc g"^). Since the u-channel A diagram is missing and 
the magnetic coupling is somewhat smaller than previously assumed, the usual huge enhancement of 
(3mi is significantly reduced, and the net value is (3^j[ = 3.9. The fit in the low-energy region is not 
dissimilar to that of Hildebrandt, but a good description of the data is obtained even in the A- resonance 
region. From the point of view of deducing a^i and (3mi from the low-energy data, the main message 
again is that a good description is obtained with a larger value of /3mi than the previously accepted one. 
Unfortunately, the amplitude used by Pascalutsa and Phillips for vrN loops in this work did not have the 
correct analytical continuation above the ttN threshold, thereby rendering both the real and imaginary 
parts of those loops incorrect. The description of data between wiab ~ 150 MeV and the A(1232) peak 
published in Ref. [205] is infiuenced by these loops, and thus cannot be regarded as definitive. This may 
also account for their finding that the best results were obtained with a value of gs ~ —2.3gM which is 
much larger than that usually deduced from the E2/M1 ratio. ^° 

More recently, Lensky and Pascalutsa [184, 211] extended the previous calculation to include all 
contributions that enter at 0{e'^6^), which includes both s- and w-channel A-pole diagrams, as well as 
ttA loops. However, unlike previous work on Compton scattering in EFT, they did not use the heavy- 
baryon framework, but worked in terms of Dirac nucleons. It is beyond the scope of this review to detail 
modern developments in baryon xPT; these avoid expanding the Lagrangian in powers of 1/Mn, while 
ensuring that positive powers of M^ which would spoil the power counting do not enter (see Bernard 
[163] for a summary). The various schemes in use all yield results which have several features in common: 
they aim at a Lorentz-covariant amplitude; when expanded in powers of 1/Mn, they agree with one 
another and with HB^PT up to the order of validity of the calculation; they also agree for any term 
which is non-analytic in ?n^; but they may differ with regard to effects that are higher-order in 1/Mn. 
Ultimately, of course, LECs (which are different in the different schemes) will absorb these differences. 
For low-order calculations, however, the difference between schemes can produce significant effects in 
observables. The calculations of Refs. [184, 211] were performed in one such version of Baryon xPT. 



^In the literature on S counting, contrary to common usage in xPT, a factor of e^ is pulled out before the counting 
starts, so 0{e^S^) - 0{P^). 

^°Note also the following misprints in that paper: in the expression for the A width in Eq. (42), s + M^ — rn^ should 
be replaced by {{yG + M^^)'^ - ml)/2; Eq. (51) should contain O^" = q'^q'al"'' +l'""qW , O^" = q^qaT""" +7^"g'„9"' 
and Og'' = i€^'^°'^qaq' o. Eq. (54) should have — O7 on the left-hand side. The unwritten convention is 69123 = 1 so the 
penultimate equation of footnote (2) should be 7^''" = — ie'^''"^7^75. 
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Although the nucleons obey relativistic kinematics, antinucleons are not — and should not be — included 
in the effective theory. 

While some of the 1/Mn terms which are automatically included in the covariant case produce 
effects relevant to low-energy dynamics, e.g. terms which must be resummed to put the threshold in 
the right place in the HBxPT approach, others act as a regulator on loop integrals. They therefore have 
no meaning when separated from the higher-order LECs, except perhaps as an estimate of the likely 
size of uncertainties due to omitted higher-order physics. Nevertheless, the polarisabilities constitute 
one case in which such higher-order effects make a sizable difference in a low-order calculation. The 
covariant B^PT diagrams to be evaluated at third order include all those in Fig. 4.3, as well as those 
from Fig. 4.4 where the 7NN vertex is generated by the 1/Mn expansion of the relativistic Lagrangian. 
An example is the Dirac magnetic moment of the Lagrangian; diagrams where both photons couple to a 
nucleon line also enter the relativistic calculation of polarisabilities at LO. The integrals differ from the 
HBxPT case in that they depended on the nucleon mass as well as on the pion mass, and the net effect 
of the sub-leading mT,/M-^ terms is to substantially reduce aEi and /3mi- Once the A(1232) is added in 
(with the parameters of Ref. [184, 211], which are taken from photoproduction studies), we return to 
familiar values of uei = 10.8 and /9mi = 4.0 — but the difference from the fourth-order HB^PT case is 
that these are predictions rather than fits. Since this calculation includes all the dynamics relevant for 
i^iab < 170 MeV, it is not surprising that the data description in this domain is good (see Fig. 4.6). 
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Table 4.1: Comparison of the different mechanisms included in different EFT variants. All calculations 
include the nucleon Born and pion-pole term. 0{X^+) indicates some contributions beyond the stated 
order; pNLO indicates some effects which would be NLO in HB^PT. CT denotes the 7N contact 
interactions of Eq. (4.11); at 0{P^) the vrN contact interactions oc q also encode resonance physics. 



4.2.8 EFT fits in relation to tlie data 

At this point, it will be useful to summarise the different counting schemes in use. In HB^PT without 
the A(1232), the counting is the usual one, in powers of P/^x '^ith as P the pion mass, a typical 
momentum or the photon energy. For the multipoles to which the A is known to contribute, though, 
A^ ~ Mp is replaced by Am(= Ma - M^ 



Thomson term which is 0{P ) (since |e| 



The lowest-order contribution to Compton scattering is the 
P/A^). When the A is included in the small-scale expansion 



[200, 201], we treat A^ ~ P and call the expansion parameter e = (Aj/, P)/A^. In the low-energy 
regime, both A and nucleon propagators scale as e~^ and so ttA loops and A-pole diagrams enter at 
0{e^) along with the corresponding nucleon diagrams. Short-range contributions to a^l etc. only enter 
at fourth order, but have often been promoted to third order — partly to provide diamagnetic strength 
to counteract the A; we refer to this as modified 0{e^). In the region of the A resonance, diagrams 
which give the A a finite width must be included, but this has not been done in the e expansion. An 
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alternative counting is given by the 6 expansion [205], in which away from the resonance region A 
propagators are suppressed relative to nucleon ones by one power of 5 ~ 771^^/ Am', here the Thomson 
term is ^(e^), LO N-pole and ttN are 0{e'^6'^), LO ttA and A-pole are 0{e'^6^), and NLO N-pole and 
vrN diagrams are O^e^S"^) (see Figs. 4.2-4.4). The ingredients of several calculations that use different 
XEFT expansions are summarised in Table 4.1. 

An important difference not noted in the table is that, in order to facilitate the treatment of the pole 
region where finite-width s-channel A-pole diagrams are leading order, all work in the 6 expansion has 
used a relativistic form for the A propagator with the 7NA vertex from Eq. (4.17), even for u ^ Am- 
At sufficiently low energies, the two magnetic 7NA coupling constants, 61 from Eq. (4.16) and Qm 
from Eq. (4.17), are related by Eq. (4.18), which translates g^ = 2.9 [206, 207] to bi = 3.8; with this 
identification, the two give exactly the same contribution to /Sa/i- However, the two vertices are not 
equivalent at higher energies. If the full /Ja/i is fit to data and hence does not directly depend on bi 
(qm), significant uj/M-^ effects from Eq. (4.17) mean that a substantially higher value of bi is required 
to give broadly the same cross section. Numerical studies show that values of ^f^j = 2.9 and bi = 4.8-5.0 
give similar cross sections up to 150 MeV, if the parameters and other ingredients of the modified e^ and 
modified e^(5^ are the same. (The closeness of the dynamical polarisabilities can be seen in Fig. 4.7.) 
Hildebrandt et al. [14, 90, 92] found that a slightly lower value of bi ^ 4.7 gave the best results, but 
they use pole-position rather than A-resonance parameters. 

Figure 4.6 shows the results of the various calculations and fits described in the previous section; 
see Table 4.1 for ingredients and Table 4.2 for a summary of results for the static polarisabilities. In 
attempting to extract polarisabilities from the low- and possibly intermediate-energy data, there are 
competing pressures. On the one hand, EFTs are clearly most reliable at low energy. On the other 
hand, the effects of the polarisabilities grow with energy, and the more data that can be included in a 
fit, the better the resulting statistical errors are expected to be. Using HB^PT to (!?(P^) but without 
the A(1232), Beane et al. applied a cut of 180 MeV in both cuiab and a/— t, the latter meaning that 
the maximum energy was in fact only reached for angles up to 60°, dropping to 90 MeV at 180°. The 
resulting extraction was therefore free from any of the issues discussed above concerning the data at 
intermediate energies, but the statistical errors of ±1.1 on a^l and (3j^j[ compared unfavourably with 
those obtained, for instance, by Baranov et al. [212] using all data up to 150 MeV prior to Olmos de 
Leon (±0.8 on a^l and ±0.9 on (3^/i). Hildebrandt et al. did include the A in the small-scale expansion 
and fit to the Olmos de Leon and Hallin data up to 240 MeV. Without a finite width, their A-pole 
amplitude starts to deviate from a Breit-Wigner-like shape which could continue into the resonance 
region above 200 MeV, more strongly at backward angles. Some of the nice agreement with the Hallin 
data in the intermediate-energy region shown in Hildebrandt's thesis may therefore be fortuitous. The 
relatively high errors (of ±1.4 on a^l and /3m{) may be due to the fact that they did not allow the 
normalisation of the data sets to float or that they did not include the older data. 

Pascalutsa and Phillips resum A self-energy diagrams, however, generating a width and a Breit- 
Wigner-like propagator, and the resulting form gives at least a qualitatively good description of the 
data up to the pole region. All the ingredients for an extraction that also incorporates the effects of 
ttN loops, ttA loops and the A-pole diagrams — with a finite width where necessary — now exist. While 
Hildebrandt et al. could not access the resonance region and hence had to fit the crucial 7NA magnetic 
coupling constant to the low- and intermediate-energy data, in the 6 expansion it should be possible to 
fit purely A parameters in the resonance region, leaving only a^l and (3^1 to be fit to low-energy data. 
Such a strategy should also be able to bypass the problems of the intermediate region, although the 
issue of Blanpied (LEGS) versus Wolf (MAMI) cannot be avoided. Details of such a calculation will be 
presented in Section 4.4. 
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Figure 4.6: (Colour online) Previous results. Working from the lowest cross section up in the last 
frame, the lines are as follows: Medium dash (orange): OiP^^ Babusci et al. [78]; short dash (red): 
C(P^) Beane et. al. [199]; solid (purple): partial 0(e^(5^+) Lensky and Pascalutsa [184]; long dash 
(blue): Oie") Hildebrandt et al. [14]; dash-dot (green): C(e^(52+) Pascalutsa and Phillips. [205]. For 
Refs. [14, 199], the fits without imposition of the Baldin sum rule are shown. For Ref. [205], the coding 
error mentioned in the text has been corrected, and photoproduction values are used for h\ and 62 
[206, 207]. The data key is in Table 3.1. The 180° curves are all extremely close to the 155° ones and 
hence are not shown. 



4.3 Comparing x^FT with Dispersion Relations 



We compare the values for various polarisabilities extracted from various DR and xEFT calculations 



(p) 



5(P) 



in Table 4.2. Overall, the agreement is very good. We note that the integral for a^l — (3j^^[ in the 
hyperbolic dispersion relation is somewhat sensitive to assumptions about high-energy physics, but the 
value found there is within the error bar of the various EFT fits. The result for a^l + (3^{ is even more 
sensitive to high-energy physics, which makes it remarkable that the agreement with the 0{P^) xPT 
prediction is so good. The apparent variation of the static values of the spin polarisabilities may be 
settled in the near future by a series of ongoing and planned experiments with polarised targets and 
beams (see Section 6.1). 

Far more informative than the static {u = 0) values of the polarisabilities are their full functional 
forms. These were calculated in (modified) 0{e^) by Hildebrandt et al. [14] and compared with the 
results of the DR analysis of Drechsel et al. [15], as updated in Ref. [14]. In Fig. 4.7 we show results 
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Table 4.2: Previous calculations of proton polarisabilities. f indicates that the Baldin sum rule was 
used. * indicates that the result was obtained by fitting to 7p data. In the (9(P^) column, the second 
errors are theory errors; the errors of the modified 0{e^) calculation are statistical only. In the fixed- 
t DR column, the second errors are from floating normalisations (see Section 4.1) and the numbers 
for spin polarisabilities are obtained by averaging the calculations of Refs. [88, 89] with uncertainties 
reflecting the range. The DR calculations presented here predict all spin polarisabilities, in contrast to 
the fit of Ref. [15] discussed in Sec. 4.1. 
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Figure 4.7: (Colour online) Comparison of DR results [15] (green long-dashed line) and (9(P^), modified 
0{e^) and modified 0{e^5^) EFT results (red dotted, and black dashed and solid lines respectively) for 
the real parts of the six dipole proton polarisabilities, as a function of cm energy. The 0{e^5^) EFT 
results use the parameters uei = 10.7, /3mi = 3.1 and hi = 3.66 as determined in the Baldin-constrained 
fit of Section 4.4, whereas the 0{e^) result uses the same polarisabilities but a non-relativistic A(1232) 
propagator and hence 6i = 5 as discussed in Section 4.2.8. The threshold for pion production and 
the peak of the A(1232) Breit-Wigner cross section are marked by Ut,- and cja respectively. Note the 
difference in scales and the inherent theoretical uncertainties of each approach. 
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for EFT fits with the parameters which will be obtained from the proton Compton data in Section 4.4. 
Keeping in mind the inherent theoretical uncertainties of the DR and EFT approaches, the shapes of 
all six dynamical dipole polarisabilities agree remarkably well up to photon energies of around 200 MeV 
(centre-of-mass). For the modified 0{e'^6^) EFT which includes a A(1232) width, the shape agreement 
continues above 200 MeV for all but the two mixed spin polarisabilities. However, for most energies, 
these are numerically small so that higher-order effects may be more prominent. The similarity em- 
phasises the point that the energy dependence of Biiy^t) is driven by near-threshold dynamics, even if 
the value of Bi{0,t) is sensitive to higher-energy physics. Its origin is the fact that HBxPT captures 
the large Eq^ photoproduction multipole that generates a significant El excitation of the pion cloud 
of the nucleon. This is clearly seen in the cusps at the pion-production threshold, most prominent in 
the polarisabilities of ElEl multipoles. We see in the table that the static 0{P^) value for a^l agrees 
exceptionally well with the DR result, while that for Jmimi i^ i^ reasonable accord with DR. 

However, the Ml excitation is not properly described in xEFT until the A(1232) is included. This 
is reflected in the improved agreement between DRs and the xEFT calculation with dynamical A(1232) 
in the second block of the table, compared to the A-less calculation of the first block, and even more 
so by the prominence of the A(1232) resonance in the polarisabilities containing an Ml multipole in 
Fig. 4.7. The poor agreement around the A peak for 'Jmie2 is rather surprising in this regard, since 
the E2 electric transition vertex is included in the EFT Lagrangian (4.17). The pronounced structure 
in the DR description in the 160-250 MeV region is markedly different from what is seen in the xEFT 
calculation; this discrepancy needs further study. 

Finally, even a 0{e^5'^) xEFT calculation does not have the strong tttt interactions in the t-channel 
which generate a substantial contribution to a]^[ — /3|'^{ that offsets the effect of the A(1232) in that 
quantity. Whether this contribution is modelled as a a or constructed through the t-channel DR (4.6), it 
is present in the numbers given in the last two columns. It is incorporated in current xEFT calculations 
only as a contribution to a a^[ — /3]^j{ 7N contact operator. Therefore, the fact that the static and 
dynamical scalar polarisabilities in xEFT and DR agree suggests that the energy dependence of this 
process is negligible below about 300 MeV, which is consistent with the phenomenology of the correlated 
27r state. 

4.4 A new fit to the proton data 

As outlined in Section 4.2.8, in the 5 expansion, it should be possible to fit purely A(1232) parameters 
in the resonance region, leaving only a]^[ and /5]i^{ to be fit to low-energy data. This section presents 
the results of a fit using this strategy. 

All the experiments have an overall normalisation error. We incorporate this by adding a piece to 
the usual x^ function: 

2 _ 2,2 

A, Xstat. "•" A,syst. 

E.ctB J /A/'j(d0-ij/d(])cxpt - (daij/d(])theoryV , V^' (J^j - iV fA 1 n^ 

where (dcrjj/df2)expt and Ajj are the value and statistical error of the zth observation from the jth 
experimental set, 5j is the fractional systematic error of set j, and Mj is an overall normalisation for set 
j. The additional parameters Mj are to be optimised by minimising the combined x^- This can be done 
analytically, leaving a x^ that is a function of a^l and 13^\^ alone. The number of degrees of freedom in 
the final minimisation is reduced by the number of independent data sets used. For more details on this 
formalism, see Ref. [212]. The best justification we have found in the modern literature for the precise 
form of this expression (specifically, the inclusion oi Nj in the denominators) is from d'Agostini [213]. 
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We perform a fit at modified third order, (9(e^(5^) plus the counterterms of Eq. (4.11), then one 
at fourth order, (9(e^(5^). For reasons which will become clear below, we regard the latter only as a 
consistency check on the former. The basic ingredients of the modified third-order calculation were 
described in Section 4.2 and shown in Figs. 4.2-4.5. These are nucleon Born diagrams and pion- and 
A-pole diagrams^ vrN loops, vrA loops and the photon-nucleon seagull terms from £*^^) which give LEG 
contributions Sa^l and S(3j^j[ to the polarisabilities. The nucleon Born diagrams are calculated to fourth 
order as given by McGovern [197], and they are indistinguishable from the full relativistic form given by 
Babusci et al. [88] for the energies of interest. The pion-pole contributions are given by Bernard et al. 
[98]; this form does not change at fourth order and is indeed relativistically invariant, although at higher 
orders, form factors would still enter at the vertices. The vrN loops are given at third order by Bernard 
et al. [98] and at fourth order by McGovern [197]. At both third and fourth order, we use Us = ^/s — M^ 
to shift the threshold to the correct place. The vrA loops are given by Hildebrandt et al. [14] in his 
Appendix B, except that we use c^Breit throughout in place of Ug and cOu to preserve crossing symmetry. 
The s- and u-channel A-pole diagrams are calculated using the Lagrangian of Pascalutsa and Phillips 
[205]; the expressions given in Appendix A3 of that paper refer to a redundant, covariant set of operators 
which we reduce to our usual six in the Breit frame. Strict consistency in the 6 counting would require 
us to use the third-order nucleon Born contribution at this order. However, since the main contribution 
at fourth order consists of terms of order u"^ in Ai and A2, polarisabilities extracted at different orders 
will be more comparable if these terms are always included. With regard to the A pole, Pascalutsa and 
Phillips have detailed which terms are leading or sub-leading in the low-energy and resonance regions; 
following their scheme, we include all terms in both regions to avoid discontinuities. At third and fourth 
order, the 7NA vertices and A propagator should strictly be expanded in powers of cj/Mn, A/Mn and 
w/A, but again, to avoid artificial shifts in polarisabilities in going from one order to the next, we keep 
the full Dirac structures. The A parameters we use are Am = 293 MeV and q^na = 1.425, which are 
fit to the Breit- Wigner peak position and width, the latter via the relativistic formula. These, and not 
the pole position, must be used if the resonance is to be reproduced. 

First we present our main result, from a modified 0{e'^6^) calculation. Our strategy is to first 
determine the 7NA couplings by considering the resonance region. By eye it is apparent that with 
typical values of a^l = 11 and (3j^j[ = 3.4 and with 62/^1 = —0.34, a reasonably convincing overall 
reproduction of the MAMI data up to the peak around 325 MeV can be obtained with 61 = 3.7 (close 
to the Pascalutsa and Vanderhaeghen value of bi = 3.76 [206]). Above the peak, the cross section falls 
off too slowly at forward angles and too fast at backward angles, but we hardly expect to fit well at 
these high energies. (See Fig. 4.8, solid blue curve.) For that reason, unlike Pascalutsa and Phillips, 
we do not adjust 62- Below the peak, the fit is always better to the MAMI data than the LEGS data 
(where they disagree) even if the peak height is raised. In view of this problem, and considering the 
fact that the former will dominate statistically, we choose to fit only to the MAMI data. 

Thus we fit 61 to the MAMI data from 200 MeV up to 325 MeV, then we fit aJ^J and /3[^i to the 
low-energy data up to 170 MeV, iterating until convergence is reached. If we exclude the Hallin data, we 
obtain a solution with a good x^- Looking at individual experiments, however, we see that the Baranov 
data at 150 MeV (which we treat as an independent data set, following Ref. [212]) have a x^ per degree 
of freedom (d.o.f.) of over 3 and should therefore be excluded from the fits. The Olmos de Leon data, 
which form nearly half the total, are fit with an acceptable x^ of 73 for 65 d.o.f. If we include the 27 
Hallin points in this energy range, though, the x^ of 41 for 27 d.o.f. is hard to accept, and we prefer to 
quote our best results without it. (The specific x^ above is taken from the fit constrained by the Baldin 
sum rule but is hardly changed if that constraint is lifted.) 

The central result for the modified 0{e'^6^) fit without using the Baldin sum rule is then (rounded 
to one decimal place) 

a^P^ = 10.5 ± 0.5(stat) ± 0.8(theory) , /3[^{ = 2.7 ± 0.5(stat) ± 0.8(theory) (4.20) 
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with x^ = 106.1 for 124 d.o.f. 
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axes of the Xmin + 1 elhpses ahgn fairly closely with these axes). These are obtained with bi = 3.66±0.03. 
With the Baldin constraint of a^ ±/3S = 13.8±0.4, the result is a^^l -P^^\ = 7.7±0.6 with x^ = 106.5 



for 125 d.o.f. and the same bi as before, which gives 



a^^l = 10.7± 0.3(stat) ± 0.2(Baldin) ± 0.8(theory) , /^^^i = 3.1 T 0.3(stat) ± 0.2(Baldin) ± 0.8(theory). 

(4.21) 
If the low-energy fits are repeated with 6i varied within its limits, a^l — (3j^j[ changes by less than 0.1. 
Varying the upper cutoff between 300 and 350 MeV scarcely changes a^l — (3^j{- Restoring the two 
omitted data sets (Hallin and Baranov 150°) also results in an upward shift in a^l — (3^j[ of 0.1. We 
therefore conclude that our fit is stable. Visually, the agreement with the bulk of the data continues 
far above the low-energy region, without any obvious systematic problems, as can be seen in Fig. 4.9. 
To arrive at the quoted theory error on our results, we note that we perform an 0{e'^6^) fit in a 
framework in which the polarisabilities first enter at 0{e^5'^). We would expect corrections to be of 
order 5^ ~ 16% of the lowest-order result, or 5 ~ 40% of the shift between the LO and NLO results; 



„(p) 



5{P) 



taking {a^[ ± /3)^^Q/2 ^ 7 to set the scale for the first approach gives 1.1, while taking the shifts in 



„(p) 



5(P) 



the values of a^{ and f3j^^ from third order to fourth order to be ~ 2 gives 0.8 in the second approach. 
In view of the similarity between our third- and fourth-order results (see later), the stability under 
inclusion or exclusion of data sets, and the values obtained in the 0{P^) and 0{e^) fits [14, 198], we 
consider the latter to be already conservative and so we use it. 
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Figure 4.8: (Colour online) Comparison of third-order (blue, solid) and fourth-order (red, dashed) cross 

sections with Compton-scattering data in the intermediate- and high-energy regions. Centre-of-mass 

cm) cross section in nb/sr plotted in bins of 8 MeV as a function of cm photon angle. In both cases, 

2.8; 6i = 3.66 and 3.47, respectively. See Table 3.1 for key. 



a^^l = 10.5 and /^S^^ 



Next we look at the 0{e S ) fit, which differs from the previous one by the inclusion of the NLO 



(p) 



j(p) 



pieces of ttN loops (in the xEFT expansion). The first observation is that, with a^l and /Sl/^ as above 
but with a somewhat reduced bi of 3.48, the general trend of the data in the resonance region and 
above is improved; the exaggerated angular dependence of the third-order cross section is tamed (see 
Fig. 4.8). However, at low energy, the fit is systematically worse, with the fourth-order curve lying 
above the third-order one — and above the data — at all angles. If we repeat the fit procedure detailed 
above, the value bi = 3.48 is confirmed, but the best fit constrained by the Baldin sum rule has a x^ of 



170; relaxing the Baldin constraint drops the x^ to 137 but with a^l + (3^[ = 18.7. In fact, this is not 



3(P) 
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a surprise: McGovern and Hildebrandt et al. independently found that adding either NLO pion loops 
or the A separately raised the cross section around the photoproduction threshold and improved the fit 
compared with third-order pions alone; adding both is too great a correction. 

With both a^l and (3^/^ as fit parameters, the primary contribution of the A(1232) at low energies 
is via the spin polarisability "Jmimi- J^^^ as inclusion of the A at third order requires promotion of the 
strictly fourth-order counterterms for aEi and /3mi, it would appear that, at fourth order, promotion 
of one or more fifth-order counterterms for the spin polarisabilities is required. Various strategies are 
possible: we could promote all four and fit them all, or we could find out which is the most important 
and fit just that one, or we could take 'Jmimi as the one with the largest A contribution. With any 
of these strategies, we can get an acceptable low-energy fit. However, there are very flat directions in 
parameter space if we promote all four, and in the absence of further constraints, it is not clear what we 
learn from the procedure. In fact, the best low-energy results from promoting a single polarisability are 
obtained with 7MIMI5 so we choose this; the Baldin-constrained results are then a^l — l3j^j{ = 7.6 ± 0.8 
and Xmimi = 2-6 ± 0.5 with x^ = 116.7; this is obtained with 61 = 3.59. 

It is reassuring to note that the fitted value of 'Jmimi is "sensible" when compared to DR estimates 
(see Table. 4.2), whereas unfitted it is distinctly large at 6.4. However, in view of the choices we have 
had to make in obtaining 'Jmimi^ ^^ would caution against treating our result as a chiral-EFT extraction 
of this parameter. Only low-energy polarisation measurements will have real power to constrain spin 
polarisabilities. It should also be noted that the fit at higher energies is degraded in the forward 
direction. For this reason, we prefer to consider the third-order results for a^l and (3^[ of Eqs. (4.20) 
and (4.21) as the most reliable, using the fourth-order fit simply for reassurance that these results are 
stable against the inclusion of higher orders. 

Our results agree within errors with the results of Beane et al. and Hildebrandt et al. (see Table. 4.2). 
The errors quoted on the former were taken from the limits of the la curve (Xmin + 2-3) as a function of 
a^l and f3^[ and hence are large compared with those obtained from the more conventional Xmin + 1 
measure used here (equivalent to marginalising over the other parameter) which explains our tighter 
errors. In general, fits which include data above 170 MeV prefer higher values of both a^l and (3^1, 
which may account for their higher central values. Like Baranov et al. [212], we find that most low-energy 
data sets, old and new, are compatible with one another (including, with the caveats of Section 3.3, 
the Olmos de Leon data set which was not available to Baranov at the time), and our results are again 
compatible with his but with a distinctly smaller a^l — f3j^j{. 

Further details of the calculation described in this section are the subject of a forthcoming pa- 
per [214]. 

4.5 Other methods 

We close this section by briefly discussing other approaches which incorporate the vrN and A(1232) 
dynamics that are key to describing 7p data in this energy domain. 

The approach of Gasparyan and Lutz [215, 216] emphasises causality and unitarity. In these works, 
Compton scattering is analysed using the xPT Lagrangian. Partial-wave amplitudes are obtained by 
an analytic continuation of amplitudes computed in (A-less) xPT in the sub-threshold region. An 
integral equation that implements vrN-analyticity and unitarity of the six Compton amplitudes is used 
to extrapolate to the kinematic region of interest. The difference from the work presented in Section 4.1 
is that, in Refs. [215, 216], the input for pion photoproduction that saturates the integral equation 
is from 0{P^) xPT, not from data. Hence, multi-pion states have not, as yet, been included in the 
s-channel integral. Moreover, the absence of an explicit A (1232) field in the Lagrangian is an additional 
dynamical assumption in this approach. 

Nevertheless, this assumption is somewhat vindicated by comparison with the data — at least for the 
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Figure 4.9: (Colour online) Comparison of third-order result with Compton-scattering data. Lab cross 
section in nb/sr plotted in bins of 10° lab angle as a function of lab photon energy in MeV. The insets 
show the fit region. The fits with and without the Hallin data are both plotted, but are indistinguishable. 
The barely visible narrow grey band shows the variation within the statistical error of the one-parameter 
fit. See Table 3.1 for key. 



specific regulator choice and 0{P^) calculation discussed in Ref. [215]. Gasparyan and Lutz can explain 
the 7p data up to energies of roughly wiab ~ 450 MeV with a quality comparable to that obtained 
in the 0{e^6^) computation (c.f. Fig. 4.9). Broadly speaking, this happens because photoproduction 
multipoles in this approach are in good agreement with extant data, and the constraints of gauge 
invariance, analyticity, and unitarity are all implemented. In contrast to the DR results displayed in 
Section 4.1 or the xEFT results discussed in Section 4.4, the 7p cross sections displayed in Ref. [215] 
are not a fit, since there are no additional free parameters in the Compton amplitude at this order. 
Because the calculation of sub-threshold amplitudes is matched to xPT, the 0{P^) results for a^l and 



Ap) 



(p) 



,(p) 



/3]^| of Eq. (4.13) are obtained. The xPT results for 'yEiM2 ^^^ 7a/is2 shown in Eq. (4.13) are also 



„(p) 



(p) 



recovered, but the values found for 'Jeiei ^^^ Tmimi ^^^ closer to those obtained in DRs. 
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Finally, we consider the dressed i^-matrix model of Kondratyuk and Scholten [217, 218], a relativistic 
approach that respects crossing symmetry. Two-body vrN-unitarity is maintained by computing the 
scattering amplitude T in a coupled-channels (ttN, photoproduction and Compton scattering) approach. 
Analyticity constraints are fully respected for one-particle reducible diagrams, but violated by the one- 
particle irreducible diagrams. The degrees of freedom used are the nucleon, the A(1232), and the p and 
a mesons, together with higher nucleon resonances. Most parameters are determined by a combined fit 
to ttN, 7p — )■ ttN and 7p data. For Compton scattering, the two key parameters are in the 770" vertex, 
and they are adjusted to reproduce the backward 7p cross section at moderate energies. In order to 
reproduce the differential cross section around pion-production threshold, and in particular the cusp 
seen there, the one-particle irreducible diagrams shown in Fig. 4.3 are added to the calculation in an ad 
hoc and approximate way. As in Ref. [215], the good description of ttN data and photoproduction up 
to tuiab ~ 450 MeV, together with the implementation of analyticity and unitarity, helps ensure a fairly 
successful result for 7p cross sections. After the fit to the data was performed, polarisabilities were 
extracted from the model, with results similar to the DR and xPT evaluations: a^[ = 12.1, I3^j[ = 2.4. 
The small value of Pmi again results from the a meson cancelling most of the A(1232) contribution. 
The diagrams of Fig. 4.3 are crucial to describing a^l and 'Jeiei^ but are only a small perturbation in 
(3^{- Together with the subsequent results of Ref. [218], this implies that reducibility and analyticity are 
not the best guides for constructing a consistent description of Compton-scattering data for wiab < 350 
MeV. It also reminds us that, especially for structure constants at higher orders in the u expansion, vrN 
and A(1232) physics play a key role and ultimately determine many aspects of the shape of Compton 
observables. 



5 Compton scattering from two- and three-nucleon systems 

Does the neutron have a similar Compton response as the proton? We saw in Section 4.2 that xEFT 
predicts the isoscalar polarisabilities to be one order in -P/A^ ^ ^ larger than the isovector ones, 
indicating that neutron and proton polarisabilities should agree at the 20% level. This is not unexpected, 
since the largest contributions to nucleon polarisabilities of all kinds come from charged-pion dynamics 
and excitation of the A(1232), and both of these mechanisms are predominantly isoscalar. 

As already mentioned, the structure functions Ai (2.11) of the neutron can be probed indirectly 
by embedding it into a light, stable nucleus. In xEFT, elastic Compton scattering has indeed been 
explored for the deuteron in Refs. [91-93, 101, 136, 199, 219-222] and for ^He in Refs. [69-71]. Since the 
deuteron is isoscalar, its Compton amplitude is sensitive only to the isoscalar polarisabilities. Neutron 
polarisabilities are then obtained by combination with the proton polarisabilities found in the preceding 
section, Eqs. (4.20) and (4.21). We first concentrate on the deuteron, since this is the only case for 
which data are presently available, and we compare these data to non-^EFT calculations in Section 5.5. 
The case of ^He is discussed in Sections 5.6 and 6.1 below. In that case, both isoscalar and isovector 
polarisabilities affect the cross section, but '^He behaves approximately as a free-neutron target for 
polarisation observables. 

In nuclei, nuclear effects such as meson-exchange currents and nuclear binding act in concert with 
the single-nucleon Compton amplitude. EFT allows the rigorous and systematic calculation of such 
effects in a model-independent way and with a well-defined theoretical uncertainty. Analysing the 
proton, deuteron and '^He in one common EFT framework can thus lead to high-accuracy determinations 
of the proton and neutron polarisabilities. However, the corresponding coherent Compton-scattering 
experiments are challenging, as described in Section 3.4, and the accuracy claimed in an EFT calculation 
must be checked by comparing to data. To overdetermine polarisabilities by multiple extractions from 
experiments on different systems thus provides important cross-checks on both experiment and theory. 
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5.1 x^FT for few-nudeon systems 

The effective theory of few-nucleon systems is significantly more invilved than that discussed in Sec- 
tion 4.2, as the shallow binding of light nuclei complicates the picture. The S-wave NN scattering 
lengths a(^Si) f» 5 fm and a(^So) ~ —24 fm are much larger than typical large-distance scales of the 
one-nucleon sector, such as the pion Compton wavelength. Concomitantly, the deuteron binding energy 
B^ ^ 2.225 MeV is small compared to the typical QCD energy scale, and the corresponding "binding 
momentum" (inverse size) 7 = \^M^B^ p^ 46 MeV ~ l/a(^Si) is appreciably less than even the typical 
"chiral" xEFT scale m^. Therefore, low-energy S-wave NN rescattering is nonperturbative at low en- 
ergies. Iterates of the NN potential are not suppressed by chiral symmetry and must at least partially 
be resummed. 

In practice this is achieved by defining a leading-order NN potential, which is then iterated via the 
Schrodinger equation, in order to generate the LO wave function of the nuclear bound state. In xEFT, 
at LO, the long-range part of the NN potential is given by one-pion exchange, constructed from the 
chiral Lagrangian, Eqs. (4.8) to (4.10). Weinberg proposed [223, 224] that this long-range potential 
should be supplemented at LO by two S-wave contact interactions, whose strengths are determined 
from the above-mentioned fine-tuned physical scales. This potential yields a Hamiltonian which is 
unbounded from below, and so a cutoff must be placed on the Schrodinger equation in order to obtain 
physically sensible predictions. Refs. [225-229] have pointed out that additional contact interactions 
(e.g. in attractive P waves) must be added to this LO potential in order to make the NN phase shifts 
independent of the cutoff used in solving the Schrodinger equation. A number of proposals for dealing 
with orders beyond LO also exist (see, e.g., Refs. [230-234]). The question of how to build a xEFT for 
NN scattering that is valid over a large range of cutoffs must be regarded as presently unresolved. 

However, Epelbaum and MeiBner argue in Ref. [235] that the xEFT proposed by Weinberg for few- 
nucleon systems is consistent if the cutoff is kept at < 800 MeV, namely roughly the mass of the p meson. 
Weinberg's proposal is that all diagrams contributing to the nuclear potential are classified according 
to their xPT order. For example, the NLO potential includes two-pion exchange (constructed from the 
first-order vertices in Eq. (4.10)), together with additional contact interactions. The entire potential at 
a fixed order is then inserted into the Schrodinger equation to obtain the nuclear wave function. For 
cutoffs in the range 500 to 800 MeV, NN potentials based on this power counting have been constructed 
to N'^LO, with very little cutoff dependence observed in this range [236, 237]. Furthermore, this N^LO 
potential produces a x^ psr degree of freedom with respect to the NN data which is comparable to that 
of "high-precision" NN potential models. Consistent three-nucleon forces have also been obtained to 
N^LO [238-240]. This method has been successfully applied in a number of reactions on deuterium: 
7rd, e~d, . . . ; see e.g. [168] for a recent review. It is this scheme that we shall use to construct xEFT 
deuteron wave functions. Additional details regarding the pros and cons of this procedure can be found 
e.g. in [164-167]. But, as far as Compton scattering is concerned, we will demonstrate in Section 5.2.4 
that any NN potential can be used for deuteron Compton scattering at the energies and level of accuracy 
of the present data — as long as that potential captures the correct long-distance physics of one-pion 
exchange and reproduces NN scattering data reasonably well. The advantage of Weinberg's power 
counting for the NN potential and operators of nuclear Compton scattering is that the consequences of 
chiral-symmetry breaking are included in the calculation in a straightforward manner, and a potential 
of the minimal complexity needed to obtain results at the desired level of accuracy is employed. 

5.2 Compton scattering from the deuteron in x^FT 

5.2.1 Scales and regimes 

We now turn to the presentation of xEFT calculations of 7d scattering. As before, we begin by 
examining the xEFT variant without dynamical A(1232) degrees of freedom and postpone including 
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the A to Section 5.2.3. Thus we first count in powers of the generic scale P and still define the 
leading-order nuclear Thomson term^^ as e^ ~ P^. 

The deuteron binding momentum is not the only new scale in few-nucleon systems. In the one- 
nucleon sector of Section 4.2, recoil effects are suppressed because the nucleon mass is much larger 
than the energy of any pions or photons participating in the reactions of interest. Recoil between two 
nucleons cannot be neglected, however, because they have the same mass. The nonrelativistic kinetic 
energy of a nucleon with momentum p propagating close to its mass shell in the few-nucleon system is 
thus E ^ p^/(2Mn) ~ P^ <^ |p| ~ p. After interaction with a photon, its propagator is given by 

'^" (5.1) 



M^E ± Mnw - {p± kf 

which converts the photon energy into a new NN momentum scale \/M^uj ^ u which is "hard" relative 
to uj. The (+) sign applies to photon absorption by the nucleon and the (— ) sign to photon emission. 
This scale only appears when the photon's energy-momentum flow must be routed through a nucleon 
and not when two photons couple instantaneously to the same nucleon; cf. discussion of Fig. 5.1 below. 

While a fourfold expansion in the low xEFT scales w, a/Mncj, 7 and m^ is possible, it is more 
convenient to approximately identify scales until only one is left. We set 7 ~ 777,7^ and define two different 
regimes of photon energy, relative to energy scales built out of the chiral scale P [199, 220-222]. 

Regime I has comparable chiral and hard photon scales, ^/M^uj ~ 711^^ ~ P, i.e. u; ~ 20 MeV, so 
that the intermediate-state propagator (5.1) scales as P~^. This Regime is discussed in Section 5.2.2. 

Regime II, on the other hand, treats the photon energy as a soft scale which is close to the 
chiral one, uj ~ 140 MeV ~ P, so that the hard scale is a/Mnw ~ P2 and the intermediate-state 
propagator (5.1) counts as P~^, i.e. one order less. Numerically, this implies ^JM^UJ ~ 360 MeV, which 
is still small compared to the breakdown scale of xEFT when the A(1232) is included dynamically. 
But convergence for these energies is in powers of yjM^uo/ K^ and hence not as rapid as in Regime I. 
This will be discussed in Section 5.2.3, where we also include a dynamical A(1232) as we did in the 
single- nucleon sector of Section 4.2.7. 

We will see that the mechanisms probed in each Regime are quite different. The transition from 
one to the other is obviously not abrupt but rather gradual, and the central values defined here provide 
only a priori estimates. For example. Regime II involves an expansion in {'^'^,m'^)/{M^uj) which breaks 
down completely when u: ~ m^/MN < 20 MeV. The discussion in Section 5.2.3 will even show that 
its higher-order terms already have a marked effect for w ~ 50 MeV. Since the available deuteron data 
lie in that intermediate region, 49 MeV < wiab < 95 MeV, we will analyse them in Section 5.3 by a 
formulation which is applicable in both Regimes, introduced in Section 5.2.4. 

5.2.2 Very low energies and the Thomson limit 

The Thomson limit of vanishing photon energy (see Eqs. (1.1) and (2.3)) imposes a stringent constraint 
on the few-nucleon Compton amplitudes in Regime I since the photon cannot resolve the target structure 
or spin: 

limT(u;,^) = -^ e'* ■ e + 0{u) . (5.2) 

Since e^ ~ P^, the amplitude indeed scales as P^. While such a low-energy expansion breaks down 
as u approaches the first intrinsic low-energy scale of the object, it not only dominates the deuteron 
amplitude for u < Pa, but constrains the cross section throughout Regime I; see Section 5.2.4. 



^^Different power-counting notations exist in the literature, including renaming the parameter to Q. In some, the 
deuteron Compton amplitudes start at order e^Q"^ [93, 96, 241, 242], but our choice allows an intuitive translation 
between the one- and two- nucleon sectors; see also the note below Eq. (5.3). 
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The significance of tlie Thomson hmit conies from the fact that it is based on gauge invariance. This, 
in turn, imphes that the total matrix element of all electromagnetic current operators is conserved^^, 
including photon couplings to charged mesons like the pion which provide nuclear binding. Thus, Friar 
used the Thomson constraint in conjunction with the generalised Siegert theorem to avoid specifying 
the exact form of meson-exchange current operators [9, 243]. In later work, Arenhovel and Weyrauch 
constructed explicit expressions for these operators and showed that the theorem enforces cancellations 
between different contributions [244, 245]. These significant simplifications and stringent numerical tests 
form the foundation of the modern theoretical description of deuteron Compton scattering in Regime 
1 [77, 91-93, 246-251]. Ultimately, the practical importance of the Thomson limit is not only the fact 
that it is fulfilled, but also how it is fulfilled. Here, we motivate the results of Refs. [6, 9, 243-245] from 
the EFT perspective, cf. [93, 96, 241, 242]. 

An exact low-energy theorem must be observed at each individual order in the generic expansion 
parameter of any EFT. In the one-nucleon sector of Section 4.2, the Thomson limit is automatically 
fulfilled; but NN rescattering complicates the picture in few-nucleon systems. Figure 5.1 shows all 
leading-order, (9(P^), contributions in Regime I. In diagram (a), the one-nucleon seagull term is con- 



i^-^. 



.^^4 



(a): 



(b): 




I: LO (P2); II: LO (P^) 



I: LO (P2); II: NLO (P^) 



I: LO (P2); II: N^LO (P^) 



Figure 5.1: (Colour online) Deuteron Compton scattering in Regime I at LO, 0{P^): one-nucleon 
seagull term (a); intermediate-nucleon propagation without, (b), and with rescattering (c; hatched 
ellipse: Tnn)- Thick (red) line: nucleon carrying photon energy as in Eq. (5.4); dot: coupling from C)^^ 
via minimal substitution and the magnetic moment. Crossed and permuted diagrams not displayed. 
Also indicated is the order at which each graph contributes in Regimes I and II, respectively. 



voluted with the deuteron wave function, and no momentum is transferred at zero energy: 



Tseaguii(a; ^ 0) = (^c 



Mts 



e'*-e)\^d) 



Mk 



e'* ■ e 



(5.3) 



This is the Thomson limit for an individual nucleon, i.e. twice the deuteron result, and would therefore 
overpredict the deuteron cross section by a factor of 4. Sachs and Austern showed that contributions (b) 
and (c) yield an amplitude +e^ (e'* ■ e)/{2M-^) at zero energy, and so cancel half of the seagull diagram 
for the correct Thomson limit ^'^ [6]. This occurs because, in contrast to (a), the photon energy in (b) 
and (c) flows through the two-nucleon state between emission and absorption, with the (highlighted) 
intermediate-state propagators of the form in Eq. (5.1): 



IMk 



o;^ ± Mncu — 7^ — 5^ 



(5.4) 



The (+) sign applies to the diagrams in Fig. 5.1(b) and (c), while the (— ) sign applies to the "crossed" 
graphs. In Regime I, Mncj ~ 7^ ~ P^, the propagator scales as P~^, and loop momenta q ^ P dominate. 
Intuitively, the initial coherent two-nucleon state is not perturbed very much by photon absorption and 
propagates coherently on a typical length scale 1/q ^ I/'tItt, which is larger than the anomalously large 
NN scattering lengths. Multiple NN interactions are therefore not parametrically suppressed before 



^^Gauge invariance and current conservation are identical in the absence of on-sliell photons in loops, as in all Compton 
calculations to date. 

^■^A 0.1% correction from relativistic effects at higher orders replaces 2Mn with the deuteron mass, M^ = 2A/n — B^- 
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another photon is emitted to produce the final deuteron-photon state, and rescattering must be re- 
summed. Power counting bears this picture out. Diagram (c) contributes at LO, just hke (a) and 

(b) in this Regime, since the scattering amphtude of two nucleons close to their mass shell scales as 
Tnn ~ P~^, as in any EFT with fine-tuned NN scales [93, 96, 164, 165, 241, 242]. 

Clearly, the three contributions are of very different computational complexity: (a) is analytically 
known in the Thomson limit (5.3), irrespective of the deuteron wave function used; (b) contains a 
convolution of the propagator (5.4) with I^E'd), which is usually not known in closed form; (c) involves 
an off-shell rescattering matrix Tnn and thus depends not only on Y^a)., but also directly on all NN 
partial waves and the potential used to generate them. In addition, only the sum of diagrams (b) and 

(c) is independent of the UV regulator. Since the Thomson limit requires that the sum of the last two 
must give — ^ of the seagull, it can be used as a nontrivial check of numerical evaluations. This condition 
must hold irrespective of the particular potential or wave function chosen and does not even depend on 
the particle content of the theory. After recent numerical improvements in the xEFT implementation 
of Hildebrandt et al. [91-93], the cancellation is fulfilled for each helicity amplitude with a relative 
numerical accuracy of < 0.01%. 

All contributions at NLO {P^) in Regime I are listed in Fig. 5.2, providing a projected relative 
accuracy of (P/A^)^ ^ 1/5^ ~ 4%. For them, the constraint is even more intricate. Since the Thomson 
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Figure 5.2: (Colour online) Deuteron Compton scattering in Regime I at NLO, 0{P^): one-nucleon 
structure terms (a); photons coupling to the same pion-exchange currents (b); or to exchange current 
and nucleon without, (c), and with, (d), rescattering. Crossed and permuted diagrams not displayed. 



limit is already fulfilled at LO, all NLO contributions must sum to zero as w — )■ 0. Equation (2.11) 
implies that the one-nucleon structure amplitudes Ai of (a) trivially obey the constraint since they vanish 
as a; ^ 0. In (b), wave functions are convoluted by six-dimensional integrals with an instantaneous 
meson-exchange current to which the photons couple directly. The other graphs involve convolutions of 
intermediate-state nucleon propagators without, (c), and with, (d), rescattering amplitudes Tnn which 
have to be constructed separately. Fortunately, Arenhovel and Weyrauch demonstrated three decades 
ago that, while (b), (c) and (d) are individually large and cutoff-dependent, they indeed cancel at a; = 
in the sum for any potential and wave function [244, 245]. Since rather different techniques are employed 
to calculate the different contributions, this provides another nontrivial cross-check. In the xEFT 
implementation of Hildebrandt et al. [91-93], it holds in each helicity amplitude to a relative accuracy 
of < 0.9%. The Thomson limit therefore provides a stringent check on the numerical implementation and 
the consistent treatment of the potential, the wave function and the NN current operator. Subsequent 
numerical improvements result in the Thomson limit now being restored to better than 0.2%. 
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5.2.3 Beyond the Thomson hmit: reducing rescattering contributions, adding the A(1232) 

At higher energies, u ^ P ^ 140 MeV and a/Mncj ~ Pz^ the Thomson limit is not directly signifi- 
cant [96, 199, 222, 241, 242]. The struck nucleon is far off-sheU, E ~ y/M^ > g^ 7^, u'^, and at LO 
the intermediate-state propagator (5.4) becomes i/u. The nucleon behaves as if static, and the standard 
xEFT counting of the one-nucleon sector prevails. Physically, each nucleon propagates incoherently, 
i.e. as if the other were absent, because its wavelength 1/q ~ l/^/M^uJ ^ 1/(360 MeV) is much shorter 
than the rescattering scale 1/7^ 1/(50 MeV). The struck nucleon has little time to scatter with its 
partner before the second photon is radiated to restore the coherent final state. Thus, rescattering is 
suppressed and can be treated perturbatively. This leads to significant computational simplifications 
and changes the relative importance of many contributions. 

The formal power counting, indicated in Figs. 5.1 and 5.2, comes to the same conclusion. The 
seagull in Fig. 5.1(a) has no two-nucleon intermediate state and thus still scales as P^, constituting 
LO. Since an intermediate-state propagator (Eq. (5.4) and thick lines in Figs. 5.1 and 5.2) now scales 
only as i/u ~ P~^, each of its occurrences moves a diagram to one order higher. A diagram with one 
intermediate NN state (but without rescattering. Fig. 5.1(b)) is demoted from LO to NLO (P^) and the 
diagrams of Fig. 5.2(c) from NLO to N^LO (P^). Since Tnn ~ P^ now, each rescattering in Figs. 5.1(c) 
and 5.2(d) costs one additional power of P (c.f. Regime I). In addition, there is one more power of 
P than in the Regime I counting of this diagram from each intermediate-state NN propagator that is 
present. Thus, these two graphs are each demoted by at least P^, to order P^ (N^LO) and beyond. 
On the other hand, the nucleon-structure contributions Ai of Fig. 5.2(a) are enhanced from order P^ 
(N^LO) in Regime I to order P^ (NLO) since they scale with cu^/m^, see (2.11) and Section 4.2.2. The 
meson-exchange diagrams of Fig. 5.2(b) still count as NLO, P^. Contributions from loop momenta 
g ~ Pa make the counting proceed in half- integer powers of P which however only appear beyond 
NLO. 

This implies that, up to N'^LO corrections, contributions from the off-shell matrices Tnn of NN 
rescattering are absent in Regime II. This simplification was first observed by Beane et al. [222] and is 
crucial to the success of the xEFT deuteron Compton-scattering calculations in Refs. [136, 199, 222]. 
The Compton amplitude for the deuteron can then be computed as 

T^d = (*d|[T^N + r^NN]|*d) . (5.5) 

Both contributions are direct convolutions of irreducible photonuclear kernels with deuteron wave 
functions. T^n is the single-nucleon Compton amplitude, and T^nn the irreducible amplitude for 
7NN — >■ 7NN, with all NN interactions truncated after (before) the departure (arrival) of the out- 
going (incoming) photon. A complete calculation for T^n at N^LO (P^) therefore consists of only the 
diagrams of Figs. 5.1(a) and (b) and Fig. 5.2(a), and for TyNN consists of the "exchange current" dia- 
grams of Figs. 5.2(b) and (c). For other light nuclei, the only change is to replace the wave function in 
Eq. (5.5) with that of the target nucleus. 

The pioneering work of Beane et al. [222] used all NLO (P^) contributions in Regime II, namely 
the isoscalar exchange currents of Fig. 5.2(b) and the one-body mechanisms of Figs. 5.1(a) and (b), 
and 5.2(a) with the single-nucleon T^n amplitude computed in xEFT without explicit A(1232) at the 
matching order, 0{P^). The amplitude is thus complete up to corrections ~ (7^, m'^)/{M-^u), which 
indicates effects in the amplitude of a few percent at a; = 100 MeV, increasing to ~ 10% at the lower end 
of the data range. Indeed, agreement with higher-energy data is good, but the 49 MeV data of Lucas [25] 
are significantly overestimated; see Fig. 5.3. The two-body currents due to Compton scattering from 
one-pion exchange turned out to be sizable at all energies and are crucial to reasonably describe all of 
the data. 

Subsequently, Beane et al. fitted the scalar nucleon polarisabilities from deuteron data [198, 199]. 
They extended the calculation to N^LO by incorporating the P^ mechanisms for single-nucleon scatter- 
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ing discussed in Section 4.2 and the two-body mechanisms of Fig. 5.2(c). At this order, these two-body 
diagrams are indeed all contributions to T^nn with one vertex from C^.^ and one from C^.^. 

Neither implementation included the effects of an explicit A (1232) degree of freedom. As discussed 
in Section 4.2.6, the A contributes substantially in the one-nucleon sector, in particular at backward 
angles for w > 100 MeV due to the rise of Pmi to about four times its static value. It would therefore 
be expected to have a similarly large effect in the deuteron, in particular for Hornidge's SAL data at 
i^iab = 94.5 MeV [26]. Except in the one-nucleon amplitude itself, no new contributions arise with 
a dynamical A. Since the deuteron is isoscalar, there is no AN component in its wave function. 
The dynamical AA component is strongly suppressed below u ~ 2Am- No explicit A appears in 
intermediate- nucleon graphs analogous to Figs. 5.1(b,c) and 5.2 (c,d), or in two-body currents similar 
to Fig. 5.2(b). The two-nucleon part at order P^ (NLO) in the A-less theory is thus identical to the 
one at orders e'^ and 6^5'^ with explicit A. We continue to use P to parameterise contributions of the 
two-nucleon sector for the variants with or without an explicit A, and to indicate the counting of the 
single- nucleon amplitudes embedded in Fig. 5.2(a) separately. Recall that for diagrams with only pions 
and nucleons, the three power countings translate as P" ~ e" ~ 6^(5^""^, see Sections 4.2.6 to 4.2.8. 

To summarise, an explicit A(1232) does not appear in the two-nucleon contributions until high 
orders, and one can therefore insert any of the three versions in Section 4.2 for the single-nucleon 
contributions: without a dynamical A(1232) at O(P^); with explicit A at 0{e^), modified by adding 
the LECs Sa^{, Sl3^^\; and with even more pion-nucleon loops at 0(e^(5^). Since the two-nucleon sector 
must be treated to the same order, the latter case warrants including the two-nucleon contributions of 
Fig. 5.2(c). As discussed in Section 4.2, a single-nucleon calculation at order e^ is actually equivalent to 
one at order e^6^ for u ~ m^, i.e. over the whole energy range of Regimes I and II. The proton fit above 
used the resummed, covariant A propagator, but that is not necessary here. All published deuteron 
data lie below the pion-production threshold, so we can treat the nonzero A width and relativistic 
effects perturbatively and consider corrections for the pion-threshold position to be negligible [92, 136], 
in contrast to the single-nucleon case in Section 4.2. 

Hildebrandt et al. performed such a NLO calculation in Regime II, with a dynamical A(1232) at 
modified order e^ [136]. The total accuracy is e^ relative to LO, i.e. N^LO, as the first order in which not 
all contributions are consistently retained in both the one- and two-nucleon sectors. The formulation 
does not improve the accuracy found by Beane et al. [222] at lower energies, but it markedly improves 
the agreement with Hornidge's SAL data [26]; cf. discussion of Fig. 5.3 below. 

5.2.4 A unified description of Regimes I and II 

The deuteron database covers 49 to 95 MeV, i.e. energies that overlap both Regimes I and II. An 
accuracy better than 10% necessary for a meaningful extraction of nucleon polarisabilities can therefore 
only be achieved if one incorporates all the effects which are LO or NLO in either Regime I or Regime 
II, as well as the energy dependence of the A(1232) at e'^6^ ~ e^(modified) that so prominently affects 
/3mi(i^)- This guarantees both the correct Thomson limit and an accurate description at higher energy. 

This was achieved by Hildebrandt et al. [91-93]. Combining the discussions of Regime I and II above, 
we see that when all diagrams of Figs. 5.1 and 5.2 are included, one actually performs a calculation which 
is complete for the two-nucleon sector up to N^LO (P^) corrections in both Regimes simultaneously and 
can move smoothly between them. In Regime I, the only higher-order diagram is the nucleon-structure 
contribution of Fig. 5.2(a). In Regime II, one includes all corrections at N^LO (Fig. 5.2(c)) as well as 
several higher-order contributions, as in Figs. 5.1(c) and 5.2(d). In each Regime, the accuracy of the 
two-nucleon diagrams is still set by the last order for which all terms are included, namely P^ in Regime 
I and P^ in Regime II. The accuracy of the single-nucleon sector must match and thus can at most 
include all diagrams at orders P^ without explicit A(1232), or, with it, e^ and 6^5^, respectively. 

One may wonder about the predicted decrease in accuracy at lower energies, where an EFT descrip- 
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tion is supposed to be better. However, since the Thomson hmit is already fulfilled at LO, higher-order 
corrections approach zero order-by-order as u is decreased. The relative size of the total NLO correction 
is thus not given by P at low energies, but rather is proportional to P multiplied by powers of cu. As 
the importance of the Thomson limit decreases in Regime II, the relative size approaches P. 

It is this formulation which we will use to extract the isoscalar electric and magnetic dipole polaris- 
abilities from all deuteron data in Section 5.3. Since isovector polarisabilities enter only at higher orders 
in xEFT, agreement with deuteron data should be good when the isoscalar amplitudes obtained in a 
proton fit are inserted into the deuteron calculation, which is then parameter- free. This is confirmed in 
the following discussion. 

First, we assess the transition between Regimes I and II. Figure 5.3 compares two otherwise identical 
calculations in xEFT with a dynamical A(1232) at order €^6^ ~ e^ (both modified) [91, 92] with the 
same static scalar polarisabilities of Eq. 4.13. The first takes the unified approach; the other follows 
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Figure 5.3: (Colour online) Comparison of the deuteron data at 49 and 94.5 MeV with xEFT pre- 
dictions (data key in Table 3.2). Solid (blue): with both dynamical A(1232) and rescattering ("uni- 
fied approach", modified order e^ ~ e'^6^)] dotted (green): with A(1232), but without rescattering 
(strict Regime II counting); dashed (red): without dynamical A(1232), but with rescattering ("unified" 
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Refs. [136, 222] for strict Regime II counting at NLO, using only (5.5) and no rescattering. The results 
indeed converge to each other with increasing u as claimed by Beane et al. [222], but the calculation 
with the correct Thomson limit clearly provides a better description of the data at cuiab = 49 MeV. In the 
strict Regime II calculation, corrections from higher orders in (7^,m^)/(MAra;) (including rescattering) 
increase dramatically as the deuteron breakup at Ucm = -Bd is approached from above, as also seen in 
Fig. 5.4. Comparison with Fig. 4.9 shows that the deuteron and proton results are of similar size for 
> 70 MeV. The angular distribution becomes quickly skewed: forward scattering becomes weaker than 
backward scattering for cjcm ^ 80 MeV, and weaker than even the 90° cross section above 100 MeV. 

Of course, this is exactly the effect of the A(1232) discussed above. In Fig. 5.3, we compare the 
deuteron data to the parameter-free predictions of the unified calculation with two variants of the 
single-nucleon amplitudes which produce the static polarisabilities of Eq. (4.13): without a dynamical 
A at order P^, and with an explicit A at order e'^6^ ~ e^ (modified). The results are depicted in Fig. 5.3 
for the lowest- and highest-energy deuteron data. At low energies, including the A as a dynamical 
degree of freedom has no significant effect, and either variant agrees well with the data, as required by 
the decoupling theorem. With increasing energy, however, effects of the A are particularly important 
to match the angular dependence of Hornidge's SAL backward-angle data at 94.5 MeV [26]. As men- 
tioned in Section 5.2.3, this was first demonstrated in the variant without rescattering contributions 
by Hildebrandt et al. [92, 136]. Indeed, the strong energy dependence from the A excitation helps 
resolve the discrepancy between these data and earlier calculations in both xEFT (Section 5.2.3) or 
models (Section 5.5 below). An alternative remedy is to exclude the two backward-angle points since 
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Figure 5.4: (Colour online) Energy and angle dependence of unpolarised deuteron Compton scattering 
in the centre-of-mass frame in the unified approach of xEFT with explicit A(1232) and the static 
scalar polarisabilities set to the 0{P^) values a^[ = 10/3]^\ = 12.5; see Eq. (4.13). Left: at 6cm = 0° 
(solid red line), 90° (dashed blue), 180° (dotted green). Right: at the Thomson limit (dotted green), 
Wcm = 30 MeV (short-dashed green), 60 MeV (long-dashed green), 90 MeV (dash-dotted blue), 120 MeV 
(solid red). 

they correspond to a momentum transfer beyond xEFT without explicit A(1232) even at 0{P^), as 
described in Section 4.2.8 [198, 199]. 

For the rest of this section, we address some more technical issues of this approach [91-93]. Quite 
different numerical implementations were merged. The Compton amplitudes involving the one- and 
two-nucleon kernels T^-^ and T^nn of Eq. (5.5) are numerically evaluated as three- and six-dimensional 



integrals, respectively, against the deuteron wave function, following Beane et al. [222]. 
nucleon-reducible contributions may be represented in the form: 



The two- 



(^dl 



J,G'(-fid + uj)j] + J\G{-B^ - uj)J,] l^d) e 



(5.6) 



The second term represents the crossed diagrams not shown in Figs. 5.1(b) and (c) and 5.2(c) and 
(d). J denotes the NN current operator including mesonic currents and G is the (interacting) NN 
Green's function, i.e. the off-shell S'-matrix of the NN system. Therefore, Eq. (5.6) includes both free 
propagation of the NN pair and rescattering via Tnn- It is constructed in coordinate space following 
Arenhovel et al. [244, 245, 252]; cf. [253, Chapter 3.7.1] and the particularly clear presentation by 
Karakowski [249]. The fact that the cross section at the Thomson limit is now obtained to better than 
0.2% accuracy [93] when the three contributions are combined provides a nontrivial check. 

Besides providing a unified framework for energies from the Thomson limit to the pion-production 
threshold, a crucial benefit of keeping higher-order rescattering diagrams even in Regime II is that 
they also considerably reduce the dependence on the deuteron wave function; see Fig. 5.5. For the 
counting which is valid only in Regime II, Beane et al. had established that the dependence of the 
cross-section predictions on the NN wave function amounts to about 10% at 94.5 MeV at NLO [222]. 
The effect is not diminished at N^LO (-P^), as they subsequently demonstrated in Ref. [199] using wave 
functions calculated from a variety of NN potentials based on xEFT [254, 255]. Since all deuteron 
wave functions share the same long-distance physics, the variability is associated with differences in 
the short-distance physics of the photon-deuteron interaction. This is largely independent of angle 
and represents an irreducible theoretical uncertainty of calculations that strictly employ the Regime 
II power counting, severely compromising the accuracy with which a^[ could be extracted from the 
deuteron data [199]. However, as discussed in Section 5.2.2, the Thomson limit dictates wave-function 
independence as a; — ?■ 0. Figure 5.5 shows that with the zero-energy point thus fixed, the dependence 
is reduced even at 100 MeV from about ±5% to < ±0.5%, and is therefore virtually eliminated in the 
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Figure 5.5: (Colour online) Unpolarised deuteron Compton cross section and data [26] at high energies 
without (left) and with (right) rescattering effects in NLO xEFT with explicit A(1232), using the static 
scalar values from Eq. (4.13). Deuteron wave functions: xEFT at N^LO (cutoff 650 MeV, identical to 
Fig. 5.3; blue sohd) and NLO (cutoff 600 MeV with Bd = 2.175 MeV, green dotted) [254]; AV18 (red 
dot-dashed) [256]; Nijmegen 93 (red dashed) [257]. 

experimentally relevant energy range [91, 92]. Therefore, a wave function derived from any modern, 
"high-precision" potential can be used. 

The dependence on the choice of NN potential could still be of importance when the partial waves 
which enter rescattering in the NN intermediate state are computed. However, Hildebrandt et al. [91, 92] 
compared rescattering contributions derived from the AV18 potential [256] with those from a rather 
crude LO xEFT potential [258]. They found that nearly perfect agreement at the Thomson limit still 
resulted, and there was only a < 4% deviation at 90 MeV. 

Hildebrandt et al. [91-93] also combined one "high-order chiral" or "high-precision traditional" po- 
tential for rescattering with a deuteron wave function generated by another potential. These apparent 
mismatches did not affect the degree to which the Thomson limit is restored, nor did they show sig- 
nificant differences in a variety of unpolarised and polarised observables at < 100 MeV. This is an 
excellent test for independence from details of different short-distance physics and implies that issues 
of matching electromagnetic currents with both the wave function and the potential can be neglected. 
Indeed, such effects related to gauge invariance only appear two orders higher than considered here. 

All of this suggests that cross sections and polarisabilities in the present approach will be essentially 
unchanged once a fully systematic xEFT is established which is cutoff-independent over a wide range 
of cutoffs and has a strict perturbative expansion. The use of Weinberg's power counting for the NN 
potential and NN operators yields results which should differ from the result in this full theory by only 
a small amount, thanks largely to the stringent constraint imposed by the Thomson limit. 

5.3 A new fit to the deuteron data 

We now use the xEFT variant from Section 5.2.4, which is valid from the Thomson limit to the pion 
production threshold, for a new extraction of the isoscalar electric and magnetic dipole polarisabilities 
from the available deuteron Compton scattering data: Lucas (Illinois) [25], Hornidge (SAL) [26] and 
Lundin (MAX-Lab) [27]. Comparisons to the static values for the proton will then allow us to assess 
whether isovector (scalar dipole) polarisabilities are indeed small, as predicted by xEFT. However, first 
we recall the conclusion of Section 3.4 that the deuteron data are not as refined as the proton data, 
cf. Table 3.2. Limited angle and energy coverage, as well as difficulties in cleanly differentiating elastic 
and inelastic events, gives rise to large statistical and systematic errors. It therefore comes as no surprise 
that the experimental errors are at least as large as the residual theoretical uncertainties in extractions 
of the polarisabilities [91, 92, 136, 198, 199]. This is confirmed by our re-analysis. 
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We present an update of the deuteron results by Hildebrandt et al. [91, 92]. As discussed in Sec- 
tions 5.2.3 and 5.2.4, it actually represents a complete (modified) 0{e'^6^) calculation in both the two- 
and one-nucleon sectors, valid from the Thomson limit up to w ~ rriT^. All two-nucleon contributions 
are included which enter at this order in either Regime I or II, Figs. 5.1 and 5.2. The single-nucleon 
sector diagrams are listed in Figs. 4.2, 4.3 and 4.5. They contain the two short-distance coefficients 
Sa^{ and Sf3^j\ (last graph of Fig. 4.4), which, strictly speaking, only enter at order e^S^. They also 
include the A(1232), but in the deuteron calculation it is treated nonrelativistically and without a 
width, as described in Section 4.2.6. With one exception all parameters are the same as for the proton 
fit, including the A parameters A^/ = 293 MeV, g^f^A = 1.425. The exception is the 7NA coupling 
for which we use the translation to the nonrelativistic value 61 = 5, as explained in Section 4.2.8. This 
treatment of the A(1232) is justified by the excellent agreement with the relativistic, nonzero-width 
approach in the region of the deuteron data, wiab ^ 100 MeV; see the comparison of the energy depen- 
dence of the polarisabilities in the two approaches in Fig. 4.7. The final fits are insensitive to varying 
the value of 61 by as much as 5%. In practice, we use the xEFT deuteron wave function at N^LO 
(cutoff 650 MeV) in the implementation of Epelbaum et al. [254] and the AV18 potential [256] for NN 
rescattering. As discussed in Section 5.2.4, this combination provides an adequate xEFT representation 
of the two-nucleon system and the residual dependence on the deuteron wave function, NN potential 
and numerical implementation is minuscule. These parameters are also used in Figs. 5.3 to 5.5 and 5.7. 

The presentation here differs from that of Hildebrandt et al. [91, 92] in the following: a new parameter 
set (61, (^vrNA; A^/) for the A(1232) from the Breit-Wigner parameters and the proton Compton data, 
and not from the pole position; a more careful analysis of systematic errors and summation of statistical 
and point-to-point systematic errors in quadrature, as discussed in Section 3.4; an implementation of 
correlated systematic errors of the experiments by a floating normalisation as in the proton case (4.19); 
and small technical improvements to implement the Thomson limit and increase numerical accuracy, 
see Section 5.2.2. However, the relatively large experimental errors make the flt rather insensitive to 
the procedure used, and so none of these changes substantially alters the conclusions of Refs. [91, 92]. 

As in the proton case, the fractional error in the extracted polarisabilities is 6"^, namely one order 
past the last one for which all contributions are consistently retained in both the one- and two-nucleon 
sectors. Up to higher-order corrections, the results should also agree with those of the version without 
rescattering but with a dynamical A(1232), fltted in Regime II, u ~ m^. Refs. [92, 136] performed 
flts of a^{ and (3lj\ to the same data without rescattering. However, the assumption tu ~ ?7i^ fails at 
the lower end of the data, as demonstrated in Fig. 5.3, and the residual wave-function dependence of 
the extracted values is rather large, namely ±1, cf. Fig. 5.5. Therefore, we do not use this criterion to 
estimate higher-order corrections, even though we note that our flnal values given below agree with theirs 
within error bars. Ultimately, higher-order effects in the chiral amplitudes at order e^(5^ and beyond 
provide by far the largest residual uncertainty, estimated at ±0.8. As in Section 4.4, this estimate is 
justifled by assuming that the uncertainty is of order 6^ of the LO-to-NLO correction; cf. also [136]. 

With the anticipated accuracy established, the isoscalar spin-independent polarisabilities read: 

aJl = 10.5 ± 2.0(stat) ± 0.8(theory) , /3l^\ = 3.6 ± l.O(stat) ± 0.8(theory) (5.7) 

from a x^ fit, with the corresponding cross sections shown in Fig. 5.6. As before, the statistical error 
is determined from the projection of the Xmin + 1 contours onto the a^[ and P^j\ axes. For 29 data 
points and 5 free parameters, one flnds a total x^ = 24.3, or 1.01 per degree of freedom. The small 
differences compared to the results reported in Refs. [91, 92], a^[ = 11.5 ± 1.4(stat) ± l(theory) and 
^Mi ~ 3-4 ± 1.6(stat) ± l(theory), are well within statistical uncertainties and stem from the minor 
differences discussed above. The collective data set appears consistent since each experiment contributes 
about equally to the overall x^; ^ind the extracted parameters remain nearly unchanged when one 
data set is eliminated. This conflrms the assessment in Section 5.2.4 that the angular dependence of 
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Figure 5.6: (Colour online) Cross sections in the two-parameter (dashed) and one-parameter (solid) 
determinations of the isoscalar spin-independent dipole polarisabilities. Bands: statistical error of the 
one-parameter fit. Data symbols from Table 3.2, statistical errors only. 



Hornidge's SAL data is described well by the theory. The floating normalisation of that set is about 
10% lower than the other two, which in turn both have a floating normalisation close to 1. All floating 
normalisations are thus compatible with the correlated systematic uncertainties (3% for Lucas, 5% for 
Hornidge, 7.5% for Lundin) stated in Section 3.4. Only more accurate experiments, such as those 
ongoing at MAX-Lab [138, 139], can reveal if the lower normalisation of the Hornidge data is more than 
a statistical accident. 

The results (5.7) are in good agreement with the proton values, confirming that isovector effects in 
api and (3mi are indeed small. They also add up, within error bars, to the isoscalar Baldin sum rule. 
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14.5 ± 0.3, Eq. (2.9). One can therefore use the Baldin constraint to reduce statistical 
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uncertainties in a one-parameter fit for a^^[ — (3^^!^ = 7.3 ± 1.8(stat) ± 0.8(theory), thereby obtaining 
very similar results with about half the statistical error: 
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10.9 ± 0.9(stat) ± 0.2(Baldin) ± 0.8(theory) 
3.6 T 0.9(stat) ± 0.2(Baldin) ± 0.8(theory) 



(5. 



The total x^ is unchanged, but drops to 0.97 per degree of freedom since the number of degrees of freedom 
is increased to 25. The cross sections of this fit, shown in Fig. 5.6, also compare well with both the two- 
parameter fit and the data. Extractions including a dynamical A(1232) [91, 92] lead to systematically 
higher values for f3^\ than those without it [198, 199], since the enhancement in the Hornidge (SAL) data 
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at backward angles must be compensated in the latter case by decreasing a^[ — /S^^, cf. Section 5.2.4. 
In addition, agreement with the Baldin sum rule was marginal in Refs. [198, 199]. 

Comparing with the results of Section 4.4, we see that the static polarisabilities of the proton, (4.20) 
and (4.21), and isoscalar nucleon, (5.7) and (5.8), are identical within error bars. They may thus be 
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combined to obtain the neutron polarisabilities in the two-parameter analysis as 

a^^l = 10.5 ± 4.0(stat) ± 0.8(theory) , /3£i = 4.4 ± 2.1(stat) ± 0.8(theory) , (5.9) 

and in the one-parameter fit using the Baldin sum rule as a^l — (3^[ = 7.0 ± 3.6(stat) ± 0.8(theory), 

a^^l = 11.1 ± 1.8(stat) ± 0.4(Baldin) ± 0.8(theory) 

, , (5.10) 

/3W ^ 4j ^ i.8(stat) ± 0.4(Baldin) ± 0.8(theory) 

with statistical errors added in quadrature. These results are in good agreement with those of the quasi- 
free experiments [68], while having somewhat smaller uncertainties. We note again that the proton 
results were obtained in a slightly different xEFT variant, in which A(1232) propagation is treated 
relativistically and its width is included non-perturbatively. However, since the omitted terms are of 
higher order and therefore parametrically small, the two variants must agree within the theoretical 
uncertainties; cf. Fig. 4.7. In turn, comparing proton and isoscalar results confirms that the scalar 
dipole polarisabilities of the proton and neutron are largely isoscalar, i.e. that both share essentially 
the same two-photon response. It is curious to note that the proton-neutron difference of the magnetic 
polarisability appears to have a statistically insignificant tendency to be slightly negative, /3]^^ ^ — 1±2, 
in line with a recent constraint from its contribution to the electromagnetic self-energy of the proton- 
neutron mass difference [18]. 

Finally, smaller theoretical uncertainties are expected from a fully consistent extraction of order e'^6'^ 
which is under way and uses the same variant to the same order for both the proton and deuteron data, 
plus chirally-consistent interactions throughout [137]. Higher-quality data with carefully formulated 
correlated and point-to-point systematic errors would reduce the sizable statistical error. 

5.4 EFTs for very low energies 

5.4.1 Compton scattering without pions 

Since u <^ m^ in Regime I, a photon cannot resolve details of pion-cloud effects and one can formulate a 
more radical EFT by integrating out the pion into contact interactions amongst nucleons and with pho- 
tons. This is the low-energy version of xEFT, called "pion-less" EFT (EFT(/)). Its typical momentum 
scale is 7 ~ 50 MeV; its breakdown scale A^ ~ m^r is set by the mass of the pion as the lightest particle 
not included as a dynamical degree of freedom; and its expansion parameter is Q^ = 7/A^. The lack 
of finite-range forces considerably simplifies calculations, which can be carried out to high orders with 
relative ease, and indeed analytic results are common. This EFT has been used for model-independent 
subtractions of "nuclear effects" and predictions in a variety of very low-energy processes, also with 
electro-weak probes; see e.g. [164, 165, 259, 260] for reviews. In such calculations, the expansion pa- 
rameter is established as being in the range Qf ~ 1/3 to 1/5 for ptyp ~ 7, and accuracy of better that 
1% has been achieved. Since gauge invariance and LETs such as the Thomson limit are automatically 
fulfilled exactly, one can check the numerics of xEFT. Indeed, any consistent theoretical framework, 
whatever its detailed treatment of meson-exchange currents, "off-shell effects" , cutoff dependence, etc., 
must agree with the EFT(^) result within the mutual accuracies at energies where EFT(|f) is applica- 
ble. EFT(^) is thus a theoretically rigorous yet numerically simple tool to study reactions at very low 
energies. 

For NN scattering, the theory is equivalent to Bethe's Effective Range Expansion of the ^Sq and ^Si 
channels [51-55], and at NLO in Q^, the scattering lengths and effective ranges are the only LECs. In the 
latter, momentum-dependent NN couplings also give rise to photon couplings by minimal substitution. 

A particularly compact NLO two-nucleon Lagrangian is given in Refs. [261, 262]. At this order, 
■^Si-^Di mixing is still absent, and the deuteron is a pure S-state. Since contributions with two isovector 
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magnetic-moment interactions are relatively large, they have typically been promoted by one order, 

(fi:(^))2 ^ 5.5 ~ Q-^^ [220-222]. 

In the one-nucleon sector, the NLO terms relevant for Compton scattering are identical to those 
in xEFT with all pion-nucleon couplings set to zero, and lead to the Petrun'kin amplitudes (2.5). 
The photons couple to a point-like nucleon with spin and anomalous magnetic moment, and structure 
effects enter through the photon-nucleon seagull terms of the effective Lagrangian (2.10). As discussed 
there, only the static polarisabilities are relevant for cu <^ A^^ (2.13). At the same order, the relevant 
two-nucleon diagrams are just those in which one or both photons couple to the "effective range" NN 
contact terms. 

EFT(^) is applicable in all of Regime I (w < 20 MeV). In that range, one can again count the 
photon energy as either soft or hard, but now relative to the low scale 7 [220, 221, 263]. When 
w < 3 MeV ~ -Bd = 7^/Mn ~ Q\ is counted as soft, the photon wavelength is larger than the deuteron 
size, and the photon energy may barely be sufficient to disintegrate it. This is the region of the very first 
Compton-scattering calculation on a composite nuclear system, by Bethe and Peierls in 1935 [4]. Here, 
one tests the stiffness of the deuteron as a whole against deformation in the photon field, parameterised 
by the electric and magnetic scalar and tensor polarisabilities of the deuteron. Tensor polarisabilities 
probe the deformation of its quadrupole component. Closed- form expressions exist for each, often to 
high orders. Since it can be shown that the polarisabilities of the individual nucleons only enter at 
the < 0.1% level [263, 264], we do not describe these in more detail here but refer to the original 
literature [164, 220, 265-267] and mention that these parameter-free predictions usually agree very well 
with other theoretical and experimental determinations [268-272]. 

At higher energies, the photon wavelength is comparable to the deuteron size, and cj ~ 7 ~ Q^ is 
counted as hard. However, the expansion is only justified if the hard scale in the intermediate-state 
propagator is also small compared to the breakdown scale, M^io < m^, i.e. cu < 20 MeV as in Regime 
I of xEFT. The nucleon polarisabilities then enter because of their uP' dependence, in principle as 
Q\ ~ 10% corrections, i.e. at N^LO. To this order, this region was studied in EFT(/), with a projected 

accuracy of order Qj; ~ 30% for a^{ and pf}^ [221, 263, 264, 273]. As the energy increases, the relative 
importance of each contribution changes. While partial-wave mixing is still absent at N^LO, two new 
short-distance NN effects enter. The first one couples the ^Sq and '^Si waves by a spin-fiipping magnetic 
photon, with its LEG L\ determined by np — ;■ d7 at thermal energies. In addition, Chen et al. [273] 
pointed out the importance of a seagull term from the relativistic correction to the nucleon magnetic- 
moment interaction, i.e. from the spin-orbit coupling. On top of the scalar amplitudes analogous to 
Al 2, they also noted that vector amplitudes dependent on the deuteron spin and analogous to ^3.6 
appear at N^LO and contribute about 20% at these energies. At leading nonvanishing order, they 
receive contributions from the magnetic moment, spin-orbit coupling and Li, and add about 4 nb/sr to 
the backward cross section at a; = 50 MeV. Since all three are magnetic, spin-fiip effects, they are most 
prominent in backward-angle scattering and thus can mimic a false signal of the combination or^\—PM\- 
In Fig. 5.7 the data at Wiab = 49 and 55 MeV are compared to the cross section of two EFT(/) 
variants which differ only by higher-order contributions and thus map out a band whose width of ^ 5% 
may indicate a typical theoretical uncertainty at N^LO [273]. However, the formal expansion parameter 
~ ^/M^/m^, ~ 1.5 at these energies suggests non-convergence even for the lowest-energy data set at 
49 MeV. In addition, in EFT(^) the deuteron is still a pure S-wave state at this order, the tensor 
components of the Compton amplitude are zero, and rescattering through higher NN partial waves 
is absent. Since all of these effects add to the cross section and are well described by pion-exchange 
contributions, it is no surprise that EFT(-^) misses overall strength compared to the xEFT version. The 
lack of the deuteron tensor component and of rescattering is most prominent at forward angles where 
the xEFT result is up to 20% higher than EFT(^), while both agree within errors at backward angles. 

Nonetheless, this breakdown should be gradual, not catastrophic, so that EFT(^) may serve for 

64 





30 60 90 120 150 180 

ftab [deg] 

Figure 5.7: (Colour online) Deuteron Compton scattering in EFT(^) and xEFT, compared with data 
at 49 and 55 MeV in the lab frame (symbols from Table 3.2, stat. errors only). Solid (blue) line: 
NLO xEFT with dynamical A(1232) and rescattering [91-93]; dashed/dotted (red): N^LO EFT(^) 
with resummed/perturbative effective-range corrections [273]; dash-dotted (green): xEFT in the KSW 
scheme with perturbative pions at NLO [221] (left panel only). All calculations use the static scalar 
polarisabilities (4.13) and are rescaled to the flux factor (2.17). 

qualitative explorations [263]. Indeed, Chen et al. fitted the scalar polarisabilities at N^LO to data 



at 49 and 55 MeV as a^^l = 12.3 ± 1.4, /3j^\ = 5.0 ± 1.6, and ag = 14.2 ± 2.1, PTn = 9.3 ± 2.5 






(statistical errors only) in two variants which should differ by higher-order terms [273]. The difference 
is again an indicator of systematic uncertainties, but the rather large values and disagreement with the 
isoscalar Baldin sum rule (2.9) reflect the fact that the fit must also mock up the physics not captured in 
EFT(^) at this order. Extractions in EFT(^) are thus not reliable for the current data, but the relative 
sensitivity of observables to the polarisabilities may be captured. For example, Chen et al. studied the 
double-polarisation observables S"J^ (2.25) and the beam asymmetry S3 (2.23) on a vector-polarised 
deuteron to N^LO [97]. Parallel results in xEFT [71, 93, 101, 274] will be discussed in Section 6. 

Can nucleon polarisabilities be extracted instead from high-accuracy data inside the formal radius 
of convergence of EFT(/), w < 20 MeV? No data exist, but HI7S can provide the high beam intensity 
for cross sections at < 10% (i.e. N^LO) accuracy needed for polarisabilities with 30% errors [99]. 
Unfortunately, this is not competitive with the accuracy of xEFT at higher energies, where signals from 
nucleon polarisabilities are much stronger; see Section 5.2.3. An attempt to decrease the error by an 
N^LO calculation runs into several unknowns, including contact terms between two photons and two 
nucleons which parameterise the part of the deuteron polarisabilities not determined by long-range NN 
properties. 



5.4.2 Perturbative pions 

Figure 5.7 also contains the result of a calculation at 49 MeV in the KSW variant of xEFT, in which 
pion effects are considered small enough to be included perturbatively [56, 275]. While it still allows 
for analytic expressions, with the Thomson limit manifestly fulfilled, it also adds some aspects of the 
pionic tensor force. However, its radius of convergence appears at best marginally larger than that of 
EFT(^) [276]. Since EFT(^) is recovered by setting the pion-nucleon coupling qa to zero and adjusting 
the parameters, it is not surprising that its results are close to resummed EFT(^). Indeed, Refs. [219- 
221] also provided the first "pion-less" Compton results. In Fig. 5.7, we rescaled the NLO calculation 
of Ref. [220] from the nonrelativistic flux factor to the version in Eq. (2.17) used for the other curves. 
This 5% change allows the dynamical content of the different theories to be isolated, based on different 
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kinematics used^^. Finally, Chen's investigation of unpolarised scattering off a tensor-polarised deuteron 
in the KSW approach [219] led Karakowski to address the same quantity in a potential model [249, 250]. 

5.5 Deuteron Compton scattering: model calculations 

Several calculations of deuteron Compton scattering were performed in the 1950s and 1960s [277-281]. 
The earliest were based on the impulse approximation, which is poor because of the large size of the 
exchange currents depicted in Fig. 5.2(b). In the 1980s, Weyrauch and Arenhovel performed detailed 
calculations of deuteron Compton scattering that included these effects [245]. They also performed an 
energy-dependent multipole decomposition of the deuteron amplitude, akin to that for the nucleon am- 
plitude in Section 2.2. From these dynamical deuteron polarisabilities, all deuteron Compton-scattering 
observables can be reconstructed, with predictions as a function of u available in Ref. [245]. 

Formulating the problem in terms of these quantities has the advantage that their imaginary parts 
can be obtained from the optical theorem and information on 7d — )■ NN from dispersion relations, 
parallel to the discussion in Section 4.1. This accounts for the NN-reducible piece of the Compton 
amplitude, see Eq. (5.6), and Fig. 5.1(b) and (c), as well as Fig. 5.2(c) and (d). The remaining effects 
are part of the amplitude T^d of all diagrams that do not include an NN intermediate state, Eq. (5.5), 
called the "seagull amplitude" in Ref. [245] and elsewhere. As already discussed in Section 5.2.2, 
significant cancellations between the NN-reducible and NN-irreducible diagrams lead to the correct 
Thomson limit. 

Based on this, Weyrauch employed current conservation to simplify the expressions and derived 
an explicit representation of the NN Green's function G, Eq. (5.6), using a separable NN interac- 
tion [246, 247]. This avoids DRs for the deuteron amplitudes. Subsequently, Wilbois [248], and later 
Karakowski and Miller [249, 250], used increasingly sophisticated NN potentials and added contribu- 
tions beyond those mandated by current conservation, such as leading relativistic corrections and some 
meson-exchange currents beyond the Siegert theorem. Both arrived at similar cross sections. 

The most sophisticated and most widely used calculations employing these techniques in a one- 
boson-exchange framework were carried out by Levchuk and L'vov [77, 251]. The NN current operator 
J in Eq. (5.6) and the irreducible diagrams of Tyd (5.5) are derived by minimal substitution from the 
one-boson-exchange Hamiltonian, including form factors. This avoids Siegert-like theorems and makes 
it straightforward to compute important dynamical effects, e.g. nonstatic corrections to the exchange 
propagators and A(1232) effects, both of which Ref. [77] considered. 

Turning to results, Wilbois claimed in Ref. [248] that the model dependence of the differential cross 
section due to the choice of NN potential is ~ 1% or less for photon energies up to 100 MeV. Levchuk 
and L'vov [77] found a variation of ~ 5% for different versions of the Bonn OBEPR potential [282], a 
higher cross section than those of Refs. [248-250], and more angular variation at photon energies > 70 
MeV. Several checks in Ref. [77] were used to ensure that current conservation and its consequences 
had been implemented at a reasonable level of accuracy. For example, the Thomson limit is violated by 
6%. All implementations agreed that rescattering effects in the cross section are ~ 10% at 50 MeV and 
decrease further by 100 MeV, in line with the subsequent xEFT findings of Fig. 5.3 in Section 5.2.4. 

Of all the studies discussed in this subsection, only Refs. [77, 248-251] included terms in the single- 
nucleon Compton amplitude beyond the nucleon Born terms. The first four included the static values of 
a^[ and /3]^\, but no higher functional dependence on to. While good agreement with the forward-angle 
Lucas data [25] was found at low energies, it is not surprising that the rise of the differential cross 
section at backward angles in the higher-energy SAL data [26] could not be reproduced, in view of the 



^^Incidentally, the authors of Ref. [221] meant to include the minuscule contribution from the slope parameters a^^ and 
i^Mi/ ^'^ (2-13) as the leading energy-dependent effects in the scalar polarisabilities, but made the mistake of using instead 
the values of the leading momentum dependence of the generalised polarisabilities for photons of nonzero virtuality. 

66 



importance of the A(1232) established in Section 5.2.3. The study of Levchuk and L'vov [77] imple- 
mented the most sophisticated single-nucleon amplitude of these models. A nonrelativistic reduction 
of the Dirac-Pauli Hamiltonian was performed and terms up to 0{1/M^) were kept in the resulting 
Hamiltonian. They then added the contributions from the static values of a^[, Pmii the four spin 
polarisabilities and four fourth-order polarisabilities, including the slope parameters of Eq. (2.13). The 
eight higher-order polarisabilities were taken from the fixed-t DR calculation of Babusci et al. [88]. As 
discussed in Section 2.2, this incorporates all dynamics of the Compton amplitude, up to this order in u. 
Nevertheless, they could not describe the backward-angle Hornidge (SAL) points, which led Hornidge et 
al. to infer a value of P^j^ ~ 10 — very different from I3l]\ [26]. However, this problem is now known to 
be the result of a coding mistake [283, 284]. Preliminary results from calculations that rectify this error, 
and improve on the treatment of the Thomson limit in Ref. [77], indicate that it is no longer necessary 
to employ a value of /3]|^{ that differs appreciably from I5^^j\ in order to reproduce these data [283, 284]. 
A similar model was used to calculate the inelastic reaction, 7d — )■ 7'np by Levchuk et al. [140, 285, 
286]. The full set of diagrams for this process was not computed, since the focus in these works was 
on making predictions for the inelastic deuteron Compton reaction in neutron quasi-free kinematics. 
In this Regime, the most important diagrams correspond to photon scattering on quasi-free nucleons, 
NN rescattering, and the effect of meson-exchange currents and the A(1232). The OBEPR potential 
was again employed to compute the deuteron wave function and NN amplitude. The nucleon Compton 
amplitudes were taken from the DR calculation of Ref. [287]. Using these ingredients, Wissmann et 
al. [286] showed that the cross section above 200 MeV in the neutron quasi-free peak is large enough 
for an accurate extraction of cc^l- These computations gave strong theoretical motivation to extract 
a^{ from measurements of the inelastic reaction in quasi-free kinematics, as performed in Refs. [66, 68]. 
Different NN models predict triple-differential cross sections in this region that vary by 3% [68]. From 
a theoretical perspective, the absence of a full, current-conserving calculation of 7d — > 7'np is of some 
concern. Also, as emphasised in Section 4.1, significant model assumptions exist in the DRs of Ref. [287] 
that were used to analyse these data, and the unanimity in DR predictions disappears at these higher 
energies. The impact of different 7N — )■ ttN multipole analyses in the DR construction of the single- 
nucleon amplitude was assessed in Refs. [68, 286], but this is not the only source of uncertainty in 
the DR calculation. Overall, the model uncertainty of ±1.1 quoted by Kossert et al. is probably an 
underestimate [288]. Lastly, Levchuk et al. [140] emphasise that, in the kinematics relevant for this 
reaction, the vr'' pole part, Eq. (2.6), has a marked impact on the cross section. As discussed in 
Section 3.5, Kossert et al. exploited this sensitivity in order to perform a simultaneous extraction of 
a^{ — [5^i\ and 77r° , assuming the Baldin sum-rule value (2.7) for a^l + (5^^\ [68]. 

5.6 Compton scattering from ^He 

The first computations of elastic Compton scattering on the ^He nucleus were performed in the past 
five years by Shukla (nee Choudhury) et al., using the xEFT framework without explicit A(1232) 
degrees of freedom at NLO {P^) [69-71]. Since the focus was on photon energies of 60 to 120 MeV, 
one can use the Regime II approximation of instantaneous interaction kernels, Eq. (5.5). A variety 
of potentials were employed in the calculation of the ^He wave function |\E'3hc). Effects from a three- 
nucleon irreducible kernel T^nnn begin only at N^LO, and the xEFT expansion converges well. Given 
an expected accuracy of < 20%, the authors present exploratory studies on the kinds of information 
about neutron polarisabilities that are accessible in such experiments. 

Not surprisingly, ^He has a significantly larger cross section for coherent Compton scattering than 
the deuteron, due to the presence of two protons in the nucleus. Since the proton Thomson terms 
interfere with the neutron polarisability amplitudes in Regime II, neutron polarisabilities have a larger 
absolute effect on cross sections for "^He than for the deuteron; see Figure 5.8. Compared to the proton 
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Figure 5.8: (Colour online) Sensitivity of the differential cross sections for Compton scattering on ^He 
in the centre-of-mass frame at NLO (Regime II, P^) in xEFT without explicit A (1232) on a^l at 60 
(left) and 120 MeV (right). Solid (black) curve: central value a^l = 12.2; long-dashed (blue): a^l 
dot-dashed (red): —2, dotted (magenta): +2, dashed (green): +4. From Ref. [70]. 
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and deuteron results of Figs. 4.9 and 5.4, the process scales roughly with Z^ and is indeed coherent at 
the lower end, Wcm ~ 60 MeV, while cross sections are enhanced over the single-charge results only by 
Z for 120 MeV. 

The sensitivity of the cross section to the scalar polarisabilities of the neutron was tested by includ- 
ing the higher-order terms 5a^l, S(3}^{ for variations about the 0{P^) predictions (4.13). The example 
in Fig. 5.8 shows a significant effect from the neutron electric polarisability. Varying the magnetic po- 
larisability produces effects of a similar magnitude. Experimental prospects and possible improvements 
will be discussed in Section 6, together with results for double-polarisation observables. 

6 The future 

We conclude by outlining the experimental and theoretical progress which can be expected in this area 
over the next few years. 

6.1 The case for investigating spin polarisabilities 

An exciting new frontier in low-energy Compton scattering from light nuclei is the use of polarised beams 
and targets to isolate specific polarisabilities of interest. In particular, the future holds considerable 
promise for gaining access to the dipole spin polarisabilities which have thus far been rather poorly 
constrained (see Section 2.2). The values extracted from experiments on the proton can be compared to 
those of the neutron measured in, e.g., ^He, and to different combinations of isoscalar spin polarisabilities 
accessible on the deuteron. xEFT predictions at order e^ ~ e^5^ (modified) exist for the sensitivity of 
the following observables to both scalar and spin polarisabilities: circularly or linearly polarised beams 
on polarised protons (and, perhaps less useful, neutrons) [90, 92] and circularly or linearly polarised 
beams on unpolarised or vector polarised deuterons [93]. A study of many deuteron observables was 
first conducted in the "A-less" version at NLO (P^) with strict Regime II counting in Refs. [71, 101]. 
For both the deuteron and the proton, the sensitivity to Jeiei seems particularly large for A^'"'^ 
(2.24), the difference of scattering a right- circularly polarised photon on a target polarised in either 
direction perpendicular to the scattering plane. A good signal for 7mimi is also seen in the observable 
Al,™, for which the target is polarised parallel to the beam and the photons are linearly polarised either 
in the scattering plane or perpendicular to it. For all targets, the spin polarisabilities can be reliably 
extracted at photon energies > 100 MeV, after measurements at < 70 MeV provide high-accuracy 
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determinations of the scalar polarisabilities so that they do not contaminate the residual uncertainties 
of the spin polarisabilities [93]. Since full information on sensitivities and impact of constraints like the 
Baldin sum rule cannot adequately be conveyed on paper, Ref. [93] focused only on prominent examples 
and made the complete results available as an interactive Mathematica notebook^^. 

As an example of a double-polarised observable in ^He, Fig. 6.1 shows predictions at 120 MeV for 
A^''''^. In addition to the first ^He Compton-scattering calculation discussed in Section 5.6, Shukla et 
al. [69-71] demonstrated that polarised ^He behaves as an "effective polarised neutron" in Compton 
scattering. For this case, at least up to NLO, the Compton response is dominated by the nuclear 
configuration in which the two protons are paired in a relative S-wave, and the spin of the nucleus is 
carried by the neutron. The two-body currents at NLO were found to be almost completely independent 
of the ^He target spin. These two facts lead to ^He double-polarisation observables that look very 
similar to the neutron asymmetries in Refs. [90, 92, 98]. In each panel, one of the four neutron spin 
polarisabihties in the basis of (2.12) is varied. This shows that a measurement of A^^'''^ at this energy 



should enable an extraction of the combination 7}"^ — (72 + 27^ )cos^, 
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A^:''''^ (2.24) was found to be sensitive to a different combination of neutron 






Figure 6.1: (Colour online) Sensitivity of the double-polarisation observable A"'''^ in ^He Compton 
scattering in the centre-of-mass frame at NLO (Regime II) in xEFT without explicit A(1232) at 120 MeV 
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when varying the neutron spin polarisabilities 7} (top left), 72 (top right), 73 (bottom left) and 
74 (bottom right), around the 0{P^) values (4.13). Solid (black) curve: unperturbed NLO result; 
variations ^7^ of —100% (long-dashed blue), —50% (dot-dashed red), -f-50% (dotted magenta) and 
-1-100% (short-dashed green) of the central value. From Ref. [70]. 



^^Ref. [93] contains two coding errors. All "isoscalar" variations quoted should be interpreted as variations around 
the isoscalar values for the neutron polarisabilities only. Fortunately, this leaves the conclusions qualitatively unchanged, 
albeit with all sensitivities doubled. Finally, the variation of ^eiei around its (correct) central value was implemented 
with wrong signs. 
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6.2 Experiment 

The High Intensity Gamma-Ray Source (HI7S) will soon perform its first Compton scattering exper- 
iment that is dedicated to the extraction of nucleon polarisabilities. At HI7S, the photons produced 
in a free-electron laser (FEL) are backscattered by the electrons in the FEL storage ring to produce 
a high-intensity beam of nearly monochromatic, 100% polarised photons, with either linear or circular 
photon polarisation [99]. The photon energy resolution depends on the coUimation of the incident beam, 
but resolution of 2-3% with photon fluxes of ~ 10^ Hz is readily possible. Photon energies up to 80 
MeV are already available at HI7S, and the region around pion threshold is anticipated to be accessible 
in 2 to 3 years. Because the HI7S photon beam does not involve bremsstrahlung, it is extremely clean, 
with very low background. 

Early Compton scattering at HI7S will focus on exploiting the high intensity of the beam to perform 
precision measurements on the proton and deuteron. Unpolarised measurements on deuterium will 
be done at 65 and 100 MeV using a scintillating active target [289]. For the first time, Compton 
scattering will be measured by detecting the recoil deuteron in coincidence with the scattered photon, 
an advance which is expected to greatly reduce backgrounds from inelastic scattering and other beam- 
induced processes. This new data set will overlap with the existing data of Refs. [25-27], but will have 
significantly better statistics than these previous experiments. 

Single- and double-polarised measurements are anticipated to be a high priority at HI7S [99]. For 
example, another active target will be used for the photon-beam asymmetry S3 (2.23) at 90° and 
82 MeV [290]. This experiment will isolate a^l at an angle of 90° using linearly polarised photons, 
since /3]^{ does not contribute to S3 at that angle, as discussed in Section 2.4. The lowest energy with 
published beam-asymmetry data is presently Ucm = 213.1 MeV [123-125], so this experiment represents 
a substantial increase of the kinematic range. 

Entirely new investigations of double-polarisation observables on the proton are planned at both 
HI7S and Mainz. At HI7S, a circularly polarised beam of energy 100 MeV incident on a newly de- 
signed, scintillating, transversely polarised proton target will enable a measurement of the double-spin 
asymmetry, S^''''^ (2.25), and provide access to 'Jeiei [291]- Meanwhile, experiments with a polarised 
beam and a polarised ^He target will be carried out at a photon energy of 125 MeV [292] to explore 
the neutron spin polarisabilities. HI7S would also be the ideal place to measure the double-polarisation 
observables on the deuteron discussed above. Such a programme requires a polarised deuterium target, 
although as yet there are no specific plans at HI7S for such experiments. However, we reiterate that 
they could provide complementary information on neutron spin polarisabilities to the ^He experiments 
already scheduled. 

A measurement of unpolarised Compton scattering on ^Li at 60 MeV has already been completed at 
HI7S [293], and data taking at 80 MeV is planned. With suitable theoretical interpretation, such cross 
sections could provide an alternative path to the extraction of a^{ and (3m{- The path towards such an 
extraction is further advanced for Z = 2 nuclei, as seen for ^He in Section 5.6. Another good candidate 
for complementary information is ^He, since it is both a scalar and isoscalar target, with cross sections 
that are comparable in magnitude to those on ^He. In these Helium nuclei, nuclear binding is not as 
complex as in ^Li, and so measurements of their Compton cross sections would be very worthwhile. The 
prospects for advances in the theory of all these nuclear Compton reactions are discussed in Section 6.3. 

Another avenue to measure nucleon polarisabilities has been pursued at the S-DALINAC at Darm- 
stadt. A proof-of-principle experiment was performed at low energies using a bremsstrahlung beam 
and a high-pressure ionisation chamber [294]. This chamber, which could be filled with hydrogen or 
deuterium, functioned as both the target and as a particle detector for the recoiling hadron. Even with 
a continuous bremsstrahlung distribution for the incident photon beam, the coincident detection of the 
scattered photon in a 25.4 cm x 35.6 cm Nal detector and the struck-particle recoil in the active tar- 
get volume provided sufficient kinematic over-determination to reduce backgrounds significantly. A test 
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measurement of 7p differential cross sections using bremsstrahlung endpoints of 60 and 79 MeV obtained 
data of comparable statistical quality to that of Federspiel [20] for 30 MeV < cuiab < 50 MeV. Yevetska 
et al. claim that "expected yields from an experiment based on the technique described here are at least 
an order of magnitude larger than in the experiment of [Ref. [24]]" and conclude that high-statistics 
experiments using this setup are feasible. However, no full-scale production run is planned. 

MAX-Lab at Lund is presently performing unpolarised measurements on the deuteron at 145 to 
170 MeV. These will provide the first high- accuracy Compton scattering data on deuterium at higher 
energies [295]. Data from previous runs below 115 MeV are being analysed in parallel [138, 139]. 
For both, the experimental setup at Lund is the same as that described at the end of Section 3.4. 
While proton data exist up to the A(1232) resonance (as discussed in Sections 3 and 4.4), there are no 
analogous measurements on the deuteron, so data on the neutron amplitude in this energy range are 
quite limited. One of the primary motivations for extending the Compton experiments at Lund above 
pion-production threshold is to provide benchmark data for EFT calculations that are envisioned for 
the near future (see below). 

A new synchrotron light facility, MAX-IV, is presently being constructed at Lund [296]. Under a 
proposal to dedicate one of the beamlines to nuclear physics, its 1.5 GeV electron storage ring could open 
up new opportunities. Perhaps as soon as 2014, it will be possible to perform laser backscattering via a 
technique similar to that at LEGS and produce collimated photon beams of > 10^ Hz and polarisations 
of > 70% for 130 MeV < uj\ah < 157 MeV. These could be used for Compton scattering and other 
photonuclear experiments [297]. Even at this early stage, a directed programme of Compton studies is 
being explored. 

MAMI at Mainz has also been gearing up to perform polarised Compton experiments using the 
Crystal Ball detector [298]. These experiments cover the energy range 200-300 MeV and focus on 
the spin polarisabilities of the proton. Polarised photon beams are incident on either an unpolarised 
liquid hydrogen target or a polarised frozen-spin butanol target, with the target located at the centre 
of the Crystal Ball. Higher photon energies are mandatory for this detection scheme because, as 
seen in Section 3.2, Compton scattering cannot be separated from other processes unless the recoil 
proton is detected in coincidence with the scattered photon. MAMI produces linearly polarised beams 
using tagged coherent bremsstrahlung on a diamond radiator, and circularly polarised beams via a 
longitudinally polarised electron beam incident on an amorphous radiator. Three separate polarisation 
observables are measured: E3, as well as S"''^ and T,™'^. These experiments have only recently begun, 
and so far only preliminary results have been obtained. Since MAMI already has a polarised '^He target, 
the possibility to perform similar experiments for 7^He scattering also exists [299]. Such experiments, 
in concert with those at HI7S, will provide access to the proton and neutron spin polarisabilities. 

We close our discussion of the future of experiments by reiterating a central conclusion from the 
proton analysis of Section 4.4. The unpolarised proton database is somewhat noisy, and the data that 
exist between 190 MeV and 250 MeV are contradictory and not of particularly high quality. The 
statement in Section 5.3 about improving the deuteron database also applies to the proton: data of 
higher quality with carefully formulated correlated and point-to-point systematic errors are needed. 
Theorists and experimentalists should work together in order to determine which of the sometimes- 
conflicting data sets are trustworthy and to ensure that part of the future experimental Compton 
programme enhances the set of measurements of unpolarised 7p scattering. 

6.3 Theory 

Indeed, it is notable that the last twenty years of the history of low-energy Compton scattering from 
protons and light nuclei has been marked by significant cooperation between experiment and theory 
(EFTs, DRs,. . . ). We will now provide an outlook for the theoretical advances which can be anticipated, 
bearing in mind that the future path of theoretical developments is somewhat trickier to predict than 
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that of experimental work. 

In Section 4.4, we presented results of an analysis of proton data up to energies around 350 MeV using 
xEFT with a resonant A(1232) pole as well as vrN loops at 0{e'^6^), with a first look at extending to 
(9(e^(5^) (N^LO) by including higher-order vrN loops. It was noted that in order to obtain a reasonable 
fit to the data at 0{e'^6'^), a contact term had to be added to the calculation, thereby allowing the 
value of 7M1M1 to be tuned. This analysis is presently being refined. In particular, the impact of 
other spin polarisabilities on the fit must be examined, and the effect of 7NA vertex dressing must be 
incorporated. The resulting calculation will be complete up to 0{e'^S'^) in the low-energy region and 
up to 0{e'^6^) in the vicinity of the resonance, and it is expected to produce robust constraints on the 
scalar polarisabilities a^l and /3}^{ [214]. 

In order to interpret the ongoing experiment at MAX-Lab described above, a theory capable of 
describing elastic scattering on the deuteron up to Wiab ~ 200 MeV is required, taking into account that 
pion photoproduction channels are open. In particular, as the energy increases towards the A(1232) 
peak, it is expected that resonant pion photoproduction will play an increasingly important role. The 
S expansion [205] mandates that the NN and NA channels should be treated on an equal footing in 
the vicinity of the resonance peak. But experience from the single-nucleon sector, combined with the 
power-counting arguments of Section 4.2, implies that this is only strictly necessary within about 50 
MeV of the peak. Thus, the first step for analysing data up to 200 MeV will be a perturbative inclusion 
of A effects, with pion production. However, it must be noted that a reliable calculation has to include 
the correct d7 — )■ dvr^ and 7d — > NNvr thresholds — ^just as the vrN threshold had to be corrected in the 
proton calculation. The remedy in the NN sector is well-defined, see Refs. [300-305], but more elaborate 
than for the single-nucleon sector. At higher energies, the issues of chirally consistent currents [306- 
309], wave functions and potentials [226, 228, 229, 232-234, 310] must also be revisited. As discussed 
in Section 5.2.4, they are largely absent at lower energies, thanks to the strict constraint imposed by 
the Thomson limit. 

Computing the deuteron breakup reaction 7d — )■ 7'np in xEFT is a high priority for future work. 
The calculation of Levchuk and L'vov was used to extract a^l from data on this reaction in Ref. [68], as 
discussed in Section 5.5. It is important to analyse its model dependence using a theory with a consistent 
description of Compton interactions in the one- and two-nucleon sectors, along with defensible error 
bars. The quasi-free mechanisms are independent of the xEFT cutoff for low missing momenta, but 
corrections from final-state interactions could be sensitive to details of the short-distance physics in 
the NN system [310]. xEFT therefore provides a diagnostic — the cutoff dependence of results — which 
facilitates studying the reliability of the calculation of effects in the NN sector. 

The single-nucleon amplitude described in Section 4.4 and presented in Ref. [214] will also provide an 
alternative to the dispersion-relation methods used by Kossert et al. [68] for the analysis of the quasi-free 
data. In particular, one might be concerned that the asymptotic contributions to the disper son-relation 
integrals introduces uncontrolled model dependence. The xEFT amplitude will be able to demonstrate 
the extent to which different assumptions about the short-distance 7N physics alter the extracted value 
of a El- If xEFT confirms that the kinematics of Ref. [68] are ideal for a reliable extraction of a^l from 
inelastic deuteron Compton data, then a by-product of such an analysis could be a xEFT result for 
'jn , to which these data are also sensitive (see Section 3.5). More broadly, such a calculation could 
be envisioned as a tool to search for other kinematics in order to optimise the extraction of neutron 
polarisabilities. 

Turning from the 2N to the 3N sector, we have already argued that ^He provides a fine opportunity for 
obtaining information on the neutron Compton amplitude (see Sections 5.6 and 6.1). The experiments 
at HI7S and MAMI discussed in Section 6.2 aim to use this opportunity to provide data pertinent to 
the spin polarisabilities of the nucleon, but the calculations reported in Section 5.6 for ^He were clearly 
exploratory in nature. They neither contained a dynamical A(1232) nor assessed the importance of 
rescattering in the 3N intermediate state. In combination, these effects can increase double-polarisation 
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observables by up to a factor of 2 [93] . 

In particular, the results from Compton scattering on deuterium presented in Section 5.2.4 suggest 
that rescattering significantly reduces the wave-function dependence seen in Refs. [69-71]. For the 
deuteron, rescattering increases some double-polarisation observables by up to 40% at 125 MeV [93]. It 
is therefore imperative to include the effect of rescattering for ^He in order to provide a precise xEFT 
extraction of neutron polarisabilities from forthcoming data. A Bochum-Jiilich-Bonn- Washington (DC) 
collaboration is embarking on such a calculation, initially using the "A-less" xEFT NN and NNN 
potential at N^LO. As in the case of the deuteron, consistency between the ^He wave function, the 
nuclear interaction and the nuclear electromagnetic operators is a key ingredient. The correct xEFT 
implementation in the low-energy (Regime I) power counting can again be checked by the Thomson 
and Gell-Mann-Goldberger-Low low-energy theorems for the triton and ^He. 

As for the proton and deuteron, we also expect the A(1232) to have a large effect on observables 
at wiab ^ 90 MeV. It is straightforward to include the single-nucleon amplitudes with an explicit A at 
orders e^ and 6^5^, respectively, into the ^He calculations. In contradistinction to the deuteron case, A 
exchange currents can appear, and so a consistent xEFT calculation requires the vr-exchange diagrams 
of Refs. [69-71] to be complemented by currents with As. In the kinematic domain < a; ~ m^r, the 
ttNA interactions and corrections to the 3N— )-3N Green's function do not have to be resummed but can 
be treated in perturbation theory. In the longer term, a calculation that includes channels with up to 
one explicit A(1232) and a resummation of intermediate-state interactions can be anticipated. It would 
compute the coherent reaction in a consistent manner over the entire kinematic range < a; < 350 
MeV, similar to the ongoing effort for the deuteron. Inelastic Compton scattering on '^He could also be 
used to extract proton or neutron polarisabilities in specific kinematics, as for the deuteron. 

Compton scattering from nuclei such as ^He and ^Li is becoming accessible to ab initio calculation. 
Bampa et al. recently applied the Lorentz Integral Transform (LIT) method to deuteron Compton 
scattering with a xEFT potential [311]. The currents were constructed via Siegert-like theorems. The 
results up to 60 MeV are promising. The same method is presently being applied to ^Li [312], for 
comparison with the recently completed ^Li experiment at 60 MeV at HI7S [293]. This computation 
will then be extended to higher energies by incorporating NN currents beyond the electric dipole and 
the full single-nucleon amplitude. The ability to accurately assess the impact of nucleon polarisabilities 
in Compton scattering from ^Li offers the tantalising prospect that they could be extracted there, where 
the cross section in the coherent region is about twice as large as in ^He — and roughly nine times larger 
than for the proton or deuteron. 

6.4 Conclusion 

An energy-dependent multipole analysis of Compton scattering provides important information on the 
scales, symmetries and mechanisms which govern the interactions amongst the low-energy constituents 
of the nucleon, and with photons. It facilitates exploration of the chiral symmetry and isospin depen- 
dence of the pion cloud, the properties of the A(1232) resonance, and the question of which degrees 
of freedom dominate the response of the nucleon spin to electromagnetic fields. It is in this context 
that Chiral Effective Field Theory has emerged as a systematic and reliable tool to guide, predict and 
analyse experiments. 

A window of opportunity seems to be emerging between the competing demands of theorists and 
experimentalists. With increasing nuclear mass, there is also an increase in the number and viability 
of potential target nuclei. Cross sections grow quadratically with the target charge Z at lower energies 
where the nucleons act coherently, and still linearly with Z at higher energies, u > 100 MeV. This 
makes high-accuracy experiments more feasible for heavier nuclear targets. Breakup experiments on 
such nuclei may permit the isolation of neutron properties in particular, but proton structure effects 
may also be enhanced thanks to the interference effects with more charges in the target. On the other 
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hand, interpreting such data in terms of the properties of the individual nucleons requires a model- 
independent analysis which treats nuclei by ab initio methods. This poses no difficulty for deuteron 
Compton scattering and is within reach for ^He. Given the recent progress, one can speculate that both 
^He and ^Li may be treated in the same xEFT framework. This will also provide important benchmarks 
on the accuracy with which xEFT describes nuclear binding and the charged meson-exchange currents 
which are directly probed by the Compton photons. 

In conclusion, the future is bright for Compton scattering from protons and light nuclei, both 
literally and figuratively. A number of high-luminosity facilities around the world conduct or plan to 
perform experiments that will map out the spin polarisabilities of the proton, with those of the neutron 
not far behind. More sophisticated experiments, with higher fiuxes and improved target techniques, 
combined with ab initio computations based on xEFT Hamiltonians, current operators and photon- 
nucleon amplitudes, could permit the use of nuclei up to at least A = Q. The resulting values for 
the twelve static dipole polarisabilities of the proton and neutron will provide important benchmarks 
for lattice QCD computations and models of the nucleon. The picture that ultimately emerges of the 
manner in which mechanisms driven by chiral symmetry (like pion loops) compete with other shorter- 
distance effects (like the properties of the A(1232) resonance) will provide fascinating and important 
insights into the similarities and differences of the two-photon responses of the proton and the neutron. 
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